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Abstract. Core-collapse supernovae are among the most energetic explosions in the 
universe marking the catastrophic end of massive stars. In spite of rigorous studies for 
several decades, we still don't understand the explosion mechanism completely. Since 
they are related to many astrophysical phenomena such as nucleosynthesis, gamma-ray 
bursts and acceleration of cosmic rays, understanding of their physics has been of wide 
interest to the astrophysical community. 

In this article, we review recent progress in the study of core-collapse supernovae 
focusing on the explosion mechanism, supernova neutrinos, and the gravitational 
waves. As for the explosion mechanism, we present a review paying particular attention 
to the roles of multidimensional aspects, such as convection, rotation, and magnetic 
fields, on the neutrino heating mechanism. Next, we discuss supernova neutrinos, which 
is a powerful tool to probe not only deep inside of the supernovae but also intrinsic 
properties of neutrinos. For this purpose, it is necessary to understand neutrino 
oscillation which has been established recently by a lot of experiments. Gravitational 
astronomy is now also becoming reality. We present an extensive review on the physical 
foundations and the emission mechanism of gravitational waves in detail, and discuss 
the possibility of their detections. 
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1. Overview 

Core-collapse supernovae are among the most energetic explosions in the universe. They 
mark the catastrophic end of stars more massive than 8 solar masses leaving behind 
compact remnants such as neutron stars or stellar mass black holes. Noteworthy, they 
have been thought to be extremely important astrophysical objects and thus have been 
of wide interest to the astrophysical community. The nucleosynthesis in these massive 
stars, and their subsequent explosions, are responsible for most of the heavy element 
enrichment in our galaxy. So naturally, any attempt to address human origins must 
begin with an understanding of core-collapse supernovae. 

At the moment of explosion, most of the binding energy of the core is released as 
neutrinos. These neutrinos, which we call them as supernova neutrinos in the following, 
are temporarily confined in the core and escape to the outer region by diffusion. Thus 
supernova neutrinos will have valuable information of deep inside of the core. In fact, 
the detection of neutrinos from SN1987A paved the way for the Neutrino Astronomy, 
which is an alternative to the conventional astronomy by electromagnetic waves. Even 
though neutrino events from SN1987A were just two dozens, they have been studied 
extensively and allowed us to have a confidence that the basic picture of core-collapse 
supernova is correct. 

Here it is worth mentioning that supernova neutrinos have attracted not only 
astrophysicist but also particle physicist. This is because supernova neutrinos are 
also useful to probe intrinsic properties of neutrinos as well as supernova dynamics. 
Conventionally they are used to set constraints on neutrino mass, lifetime and electric 
charge etc. More recent development involves neutrino oscillation, which have been 
established experimentally in the last decade. Neutrino oscillation on supernova 
neutrinos is important in two ways. First, since neutrino oscillation changes the event 
spectra, we cannot obtain the information on the physical state of the core without a 
consideration of neutrino oscillation. Second, since supernova has a distinct feature as 
a neutrino source compared with other sources such as the sun, atmosphere, accelerator 
and reactor, it also acts as a laboratory for neutrino oscillation. 

Supernova is now about to start even another astronomy, Gravitational- Wave 
Astronomy. In fact, core-collapse supernovae have been supposed to be one of the 
most plausible sources of gravitational waves. Currently a lot of long-baseline laser 
interferometers such as GEO600, LIGO, TAMA300 and VIRGO are running and 
preparing for the first direct observation, by which the prediction by Einstein's theory 
of General Relativity can be confirmed. 

Astrophysicists have been long puzzled by the origins of the gamma-ray bursts 
since their accidental discovery in the late sixties. Some recent observations imply that 
the long-duration gamma-ray bursts are associated with core-collapse supernovae. In a 
theoretical point of view, the gamma-ray bursts are considered to be accompanied by the 
failed core-collapse supernovae, in which not the neutron star but the black hole is left 
behind. It is one of the most exciting issue to understand how the failed core-collapse 
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supernovae can produce the observed properties of the gamma-ray bursts. 

In order to obtain the understanding of these astrophysical phenomena related 
to core-collapse supernovae and the properties of neutrino and gravitational- wave 
emissions, it is indispensable to understand the explosion mechanism of core-collapse 
supernovae. However one still cannot tell it exactly albeit with the elaborate efforts 
during this 40 years. At present, detections of neutrinos and gravitational waves from 
nearby core-collapse supernovae are becoming reality. Since neutrinos and gravitational 
waves can be the only window that enables us to see directly the innermost part of 
core-collapse supernovae, their information is expected to help us to understand the 
explosion mechanism itself. Under the circumstances, the mutual understanding of the 
explosion mechanism, the supernova neutrinos, and the gravitational waves, which we 
will review in turn in this article, will be important. 

The plan of this article is as follows. We begin by a brief description of the standard 
scenario of core-collapse supernovae in section 2. In section 3, we give a tool to discuss 
supernova neutrinos and their observation, neutrino oscillation. Although neutrino 
oscillation is thought to have only a negligible effect on the dynamics of supernova, 
it is necessary when we try to interpret observed neutrinos and extract information 
of supernova from them. Then supernova neutrinos and their neutrino oscillation are 
elaborately reviewed in section 4. With respect to the study of the explosion mechanism, 
good progress in the mult i- dimensional models has been made recently. We review these 
studies in section 5. Finally, gravitational waves in core-collapse supernovae are reviewed 
in section 6, in which we pay a particular attention to the predicted characteristics of 
gravitational waves and their detectability for the currently running and planning laser 
interferometers. So far a number of excellent reviews already exist on various topics in 
this article. This article goes beyond such reviews to cover more the state-of-the-art 
investigations. 
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2. Supernova Theory 

2.1. The Fate of Massive Star 

In a historical view point, supernovae owe their name to astronomers Walter Baarde 
and Fritz Zwicky, who already in the early 1930's realized that these objects show a 
sudden bursts in luminosity that slowly decays, similar to common novae, but much 
more luminous and rare [18]. Their high luminosities comparable to the integrated 
light of their host galaxies and their broad spectral lines led them to conclude that 
supernovae were very energetic explosions produced at the death of the massive star. 
What is amazing is that they suggested that a supernova derive their tremendous energy 
from gravitational collapse, in particular that the inner part of the star collapses to a 
neutron star. Although much observational and theoretical progress have been made 
since then, and many physical principles and important details have been identified, the 
basic picture of the early 1930's still holds nowadays. To begin with, we review the 
current understanding of the fate of the massive stars in the following. 

The fate of a single massive stars, that is to say, whether the remnant formed 
after stellar collapse will be a neutron star or a black hole, is mainly determined by its 
mass at birth and by the history of its mass loss during its evolution. The mass loss is 
expected to be crucially affected by the initial metalicity of the star, because the mass 
loss rate by the stellar winds is sensitive to the photon opacity, which is determined by 
the metalicity. The stars with high initial metalicity have more mass loss, and thus, have 
smaller helium cores and hydrogen envelopes during its evolution. Stellar collapse of 
such stars tends to lead to the formation of a neutron star, while for the lower metalicity 
stars, a black hole [123]. Figure 1 illustrates how the remnants of massive stars will be as 
a function of the initial mass and the metalicity (this figure is taken from [123]). From 
the figure, stellar collapse of the stars with the initial masses above ~ 9M and below 
~ 25M lead to the formation of neutron stars. Above 25M , black holes are expected 
to be formed either by fall-back of matter after the weak explosion (below 4OM ) or 
directly if the stellar core is too massive to produce the outgoing shock wave (above 
4OM ). Given a fixed initial mass above 4OM , the stars with smaller initial metalicity 
tend to form a black hole directly due to the more heavier core as a result of the less 
mass-loss activities during evolution. 

Recently, the fate of massive stars has been paid considerable attention. This 
is mainly due to the accumulating observations that the death of massive stars and 
supernova-like events are associated with the long-duration gamma-ray bursts (GRBs) 
(see, for example, [194]). The fact that accompanying supernovae are in general more 
energetic (they are frequently referred to as "hypernovae" in the literature) than the 
canonical core-collapse supernovae is another reason for this frenzy [225] . According to 
the most widely accepted theoretical models, it is believed that a black hole/an accretion 
disk system supported by the sufficient angular momentum is required [220]. In addition 
to the rapid rotation, the strong magnetic fields, as high as 10 15 ~ 10 16 G in the central 
regions are also pointed out to be helpful for producing the GRBs. In order to determine 
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the progenitor of the gamma-ray bursts, stellar rotation and magnetic fields should be 
taken into consideration. Such investigation has just begun [122]. In addition, the 
astrophysical details of the geometry or environment of the black hole/accretion system 
are currently hidden from us both observationally and computationally. Although 
these are open questions now, this situation may change in the near future with the 
development of gravitational- wave and neutrino observatories and more sophisticated 
astrophysical simulation capabilities (see [262] for a review). 

Very massive stars between 14OM < M < 26OM with the smaller initial 
metalicity are considered to become unstable to the electron-positron pair-instability 
(77 "~ e + e~) during its evolution, which lead to the complete disruption of the 
star. Recently, explosions of metal-poor stars have been paid great attention because 
such stars are related to the first stars (the so-called Population III stars) to form 
in the universe. So far two hyper metal-poor stars, HE0107-5240 [69] and HE1327- 
2326 [96], whose metalicity is smaller than 1/100000 of the sun, have been discovered. 
They provided crucial clues to the star formation history [281] and the synthesis of 
chemical elements [341, 150] in the early universe. Furthermore, neutrino emissions and 
gravitational waves from such stars are one of the most exciting research issues. 

In this review, we focus on the ordinary supernova which lead to the neutron star 
formation (~ 9M < M < 25M with the solar metalicity). As will be explained below, 
the most promising scenario of the explosion mechanism of such stars are the neutrino- 
heating explosion. After we shortly refer to the current status of the presupernova 
models in section 2.2 (for details, see, [364, 123]), we explain the scenario in section 2.3. 
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initial mass (solar masses) 



Figure 1. Remnants of massive single stars as a function of initial metalicity and 
initial mass. In the regions above the thick green line (for the higher initial metalicity), 
the hydrogen envelope is stripped during its evolution due to the active mass loss 
processes. The dashed blue line indicates the border of the regime of direct black hole 
formation. The white strip near the right lower corner indicates the occurrence of the 
pair-instability supernovae. In the white region on the left side at lower mass, the 
stellar cores do not collapse and end their lives as white dwarfs. This figure is taken 
from Heger et al (2003) [123]. 
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2.2. Presupernova Models 



H m 3.B1 % 10 11 cm 
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Figure 2. Structure of the precollapse star (Woosley and Weaver's 15 M© model taken 
from [362]). In the upper panel, the temperature and density profiles are given. L to t, 
e v , and e nuc , represents the total energy loss and the contribution from the neutrino 
emission and from the nuclear-energy generation, respectively. In the lower panel, the 
composition profile is given. 



In Figure 2, an example of the precollapse stellar model by Woosley and Weaver (1995) 
[363] , which has been often employed as an initial condition of core-collapse simulations, 
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is shown. The iron core is surrounded by shells of lighter elements (the bottom panel of 
Figure 2). This is called onion-skin structure. The size of iron core is the order of 10 9 cm 
while the stellar radius is larger than 10 13 cm. At the core and the surrounding shell, the 
density decreases steeply and hence the dynamical timescale of the core (:Td yn ~ (Gp)^ 1 ^ 2 
with G and p being the gravitational constant and the average density) is much shorter 
than that of the envelope (see the upper panel of Figure 2). That is, the dynamics of 
the iron core is not affected by the envelope. Therefore we focus on the core hearafter 
for a while. 

The late evolutional stage of massive stars are strongly affected by weak 
interactions. In fact, it can be seen from the upper panel of Figure 2 that the 
dominant energy loss process in the iron core is the neutrino emissions (see L tot , e u , 
and e nuc in the panel). The generated neutrinos, which are well transparent for densities 
P < 10 11 g cm -3 , escape the star carrying away energy and thus cooling the star. 
Due to the weak interactions, namely electron capture and beta decay, not only the 
core entropy s, but also the electron fraction Y e , which is the electron to baryon ratio, 
changes. Since the mass of the presupernova core can be approximately expressed by 
the effective Chandrasekar mass [66, 334], 



M Ch = 5.83? e 2 



1 + 



( 6 



M , (1) 



with Y e and s e being the average values of electron fraction and electronic entropy per 
baryon in the core, the weak interaction rates play an important role of determining the 
core mass. Putting the typical values of Y e = 0.45, s e = 0.52 in a 15M star into Eq. 
(1), one has an Chandrasekhar mass of 1.34M which is close to the core mass obtained 
by the stellar evolution calculation (see Figure 2). 

So far presupernova models have been constructed by employing the weak 
interaction rates by Fuller, Fowler and Newman (FFN) [104, 105, 106] for electron- 
capture rates with an older set of beta decay rates. As well known, the electron 
capture and its inverse are dominated by Fermi and Gamow-Tellar transitions. A 
correct description of the Gamow-Tellar transitions is difficult because it requires to 
solve the many-body problem in the nuclear structure. Due to the restricted available 
experimental data in the mid 1980's, the tabulations of FFN could not fully describe 
the Gamow- Taylor distributions in nuclei. This has been practicable recently by the 
new-shell model calculation by Langanke and Martinez-Pinedo ([189, 190], see [193] for 
review). According to Heger et al [124, 125], who studied the effect of the shell model 
rates on presupernova models by repeating the calculations of Woosley and Weaver 
(1995) [363] while fixing the other stellar physics, the iron core mass is found to be 
reduced about up to O.2M than the ones in Woosley and Weaver's computations (see 
Figure 3). 
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Figure 3. Comparison of the values of the iron core masses for 11 — 40M© stars 
between the Woosley and Weaver's models (WW) using the FFN interaction rates 
and the ones using the shell model weak interaction rates by Langanke and Martincz- 
Pinedo (LMP) [189, 190]. AMp e represents the difference of the masses between the 
two computations. This figure is taken from [193]. 



2.3. Standard Scenario of Core-collapse Supernova Explosion 

In the following, we shall briefly outline the modern picture of the explosion mechanism 
of core-collapse supernovae (see, also [35, 319] for reviews). 

2.3.1. onset of infall In the late-time iron core of a massive star, the pressure, which 
supports it against the core's own gravitational force, is dominated by a degenerate gas 
of relativistic electrons, 

Pe4(30 1/3 f-)V, (2) 



4 \m u 

where Y c = n e - — e e + is electron fraction per baryon, m u is the atomic mass unit, and 
p is the density. At the typical core densities and temperatures (p ~ 10 10 g cm" 3 and 
T ~ 10 10 K), the electron capture on Fe nuclei occurs via 

56 Fe + e~ -> 56 Mn + u e , (3) 

because the Fermi energy of electrons, 

He = {37r 2 n e )hc (4) 



i 



= ll.lMeV(-— ^ — -V (5) 
1 10 10 g cm- s ' 
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exceeds the mass difference between the nuclei, namely, m Mn — m Fe = 3.7MeV. Decrease 
of the electron fraction results in the reduction of the pressure support and the core 
begins to collapse. Note that neutrinos escape freely from the core before the central 
density for p < 10 11 g cm -3 as will be mentioned in the next subsection. 

The onset of core-collapse can be also understood by the fact that the $ adiabatic 
index: 

<91np N 



dlnp 



(6) 



is lowered below 4/3, which is the instability condition against the radial perturbation 
of a spherical star [66]. From Eq. (2), the adiabatic index becomes 



d\nP P 



dlnp 



dlnP e 



1 + 



s dlnp 
<91nK 



+ 



dlnY e 



S,Ye 



dlnp 



dlnP e 



dlnY P 



dlnS 
s dlnp 



<91nP P 



dlnS 



dlnp 



(7) 
(8) 



where the final term on Eq.(7) is set to zero, because the collapse proceeds almost 



which makes T$ 



adiabatically. Progression of electron capture implies negative 
less than 4/3. 

Furthermore, the endothermic photodissociation of iron nuclei, 

7 + 56 Fe -> 13a + 4n - 124.4MeV (9) 

occurs for the temperature T > 5 x 10 9 K, which leads to the reduction of the thermal 
pressure support. In addition, the internal energy produced by the core contraction is 
exhausted by this reaction. Both of them promote the core collapse. 

Since the degenerate pressure of relativistic electrons in finite temperature can be 
expressed as, 

2" 



127T 2 C 3 fl 3 



2 / S e 
1 + -I- 



3VttF r 



the adiabatic index in Eq. (6) becomes, 



1 + 



d In S P 



l + 2/3a 2 <91np 



(10) 



(11) 



Ye 



where S e = 7r 2 y" e /c B 7 1 /p e is the electron entropy with /c B being the Boltzmann constant 
[30] . The electron entropy decreases with the central density during infall phase because 
the photodissociation proceeds by the loss of the thermal energy of electrons. Hence, 
9 Q^ e in Eq. (11) becomes negative, by which the core is shown to be destabilized 

Ye 

by the reaction. It is noted that the entropy transfer from electron to nucleon occurs 
during the core collapsing phase because the reduction of electron entropy leads to the 
increase of the entropy of nucleons, while conserving the total entropy [30]. 

| Strictly speaking, this adiabatic index is a pressure-averaged adiabatic. See for details, section III - 
(f) in [55]. 
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2.3.2. neutrino trapping After the onset of gravitational collapse, the core proceeds 
to contract under the pull of the self gravitational force, unnoticed by the rest of the 
outer part of the star, on a free-fall time scale, which is of the order of r ~ (Gp)^ 1 ^ 2 ~ 
0.04 sec (p/1 x 10 10 g cm" 3 ) -1 / 2 with the average density of the core of ~ 10 10 g cm~ 3 . 
When the central densities exceed 10 11 — 10 12 g cm -3 , electron neutrinos, which can 
escape freely from the core at first, begins to be trapped in the core, because the 
timescale for electron neutrino diffusion from the core becomes longer than that for the 
timescale of the core-collapse. This is the so-called "neutrino trapping" , which plays 
very important roles in supernova physics [278, 279, 97]. 

During the core collapses, only electron neutrinos (z/ e ) are produced copiously by 
electron captures. Since the wavelength of neutrinos with the typical energy of E u , 

^^KliTSv)" 1 ' (12) 

is longer than the size of the nuclei of I^Fe? 

r nuc ~ 1.2A 1 / 3 fm » 5 (J^J * fm, (13) 

neutrinos are scattered coherently off A nucleons in the nucleus, by which the cross 
section (u e + A — > u e + A) becomes roughly A 2 times of the cross section of each 
scattering of nucleons [y e + n, p — > v e + n, p) . Thus the coherent scattering of neutrinos 
is the dominant opacity source for the neutrinos during the infall phase. 

The mean free path determined by the coherent scattering of the neutrinos on the 
iron nuclei is estimated to be, 

A, = — , (14) 
a A n A 

where ua = p/{Am u ) is the number density of nuclei and a a is the cross section of the 
coherent scattering in the leading order (see the detailed one in [54]), 

= A 1 -h«'** t '- 1 ^' (15) 

where a = AG 2 F (m c c 2 ) 2 / (n(hcY) = 1.705 x 10~ 44 cm 2 is a convenient reference cross 
section of weak interactions, Gp and 6 W is the Fermi coupling constant and the Weinberg 
angle. The mean electron-neutrino energy in the iron (|gFe) core can be estimated as 
follows, 

K « = | ( Stt 2 ^) 1 he » 10.3MeV ( ^) * (-^-) * .(16) 

6 6\ m u J \3 x 10 10 g cm -3 / \ 26/56/ 

Introducing Eqs. (15) and (16) to Eq. (14), the mean free path in Eq. (14) becomes 

K = — - — 
a A n A 

= 8 ' 5 X Cm ( 3 X " (I)' {wjkv) 2 ' (17) 

w io ' cm ( j^wk^ ) (£) (2^5) (is) 
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Since the mean free path becomes smaller than the size of the core, 



*- « (^f ) V3 10?Cm ( 3 x ifej 

as the central density increases (note A„ oc p -5 / 3 , while i? cor c oc p -1 ^ 3 ), neutrinos 
cannot escape freely from the core. This suggests that there is a characteristic surface 
determining the escape or trapping of neutrinos in the core. The radial position of the 
neutrino sphere is usually defined as the surface where the neutrino "optical" depth, 

r°° dr 

r(r,E v )= / -, (20) 

Jr *v 

becomes 2/3. The neutrino sphere is the effective radiating surface for neutrinos, in 
analogy with the "photosphere" of normal light emitting surface. It is noted that its 
position differs from neutrino species to species and is dependent on the neutrino energy. 
The neutrino sphere, which we are now discussing, is of the electron neutrinos defined 
by their mean energy. 

Introducing Eq. (18) to the above equation, one obtains 

poo 

r(r,E u )ocEl p(r)A(r)dr. (21) 

Jr 

Taking the distribution of the density, which can be approximated by 

p(r) = ( with H = 3 x 10 31 ) (22) 

during the collapsing phase [35], and taking the typical values of A = 56 at the central 
density of p = 10 12 g cm -3 , the optical depth becomes 



Thus the typical radial position and the density of the neutrino sphere (r(R u , E v ) = 
2/3) becomes 

fl - 53 10 x 107 cm (islv) • (24) 

p(MEu)) = 3.6 x 10 10 g cm- 3 (j^y) ■ ( 25 ) 

Neutrinos produced at R > R u can freely escape from the core, while neutrinos produced 
inside propagates outwards by a random-walk induced by the coherent scattering. The 
diffusion time for neutrinos to diffuse out from the core of size R, can be estimated as, 

tdiff = -^~2.3x lO^sec ( — -) . (26) 

c\ u \10 i2 gcm~ d / 

Since the dynamical timescale of the core, 

t dyn = 4xlO- 3 sec(^ iol2g P cm _ 3 ^ 2 , (27) 
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Figure 4. Infall velocity and sound velocity versus radius at the central density of 
10 12 g cm -3 of a 15 Mq progenitor model. The region inside and outside the sonic 
point (i? « 200 km, at which the two curves cross) roughly corresponds to the inner 
core and the outer core, respectively. 



is shorter for the core density of 10 11 ~ 10 12 g cm -3 . This also means that neutrinos 
cannot freely escape from the core and trapped in the core. After the neutrino trapping, 
the lepton fraction (Y t = Y e + Y u ), where Y u = n Ue — n Pe is the electron-type neutrino 
fraction per baryon, is kept nearly constant during the collapse stage. Furthermore, 
electron neutrinos also become degenerate like electrons and the (3— equilibrium is 
established between e~ + p — > n + v e and its inverse. After the achievement of the (3— 
equilibrium, the entropy is conserved and the collapse proceeds adiabatically. 

2.3.3. homologous collapse The collapsing core consists of two parts: the 
(homologously collapsing) inner core and the (supersonically inf ailing) outer core. This 
structure is clearly seen in Figure 4. Matter inside the sonic point (the point in 
the star where the sound speed equals the magnitude of the infall velocity) stays in 
communication and collapses homologously (velocity roughly proportional to radius). 
On the other hand, the material outside the sonic point falls in quasi-free fall with 
velocity proportional to the inverse square of the radius. Beautiful analytic studies 
have been done on this phase of collapse by [111, 365], who predict that the size of the 
homologous core is roughly the Chandrasekhar mass (see Eq. 1). For a typical value of 
Y e in the inner core, the mass of the inner core can be estimated M c h ~ 0.5 — O.8M , 
which is in good agreement with the mass of the inner core obtained in a realistic 
numerical simulation [55]. 

2.3.4- core bounce When nuclear densities are reached in the collapsing core (p c ~ 
3 x 10 14 g cm -3 ), repulsive nuclear forces halt the collapse of the inner core driving a shock 
wave into the outer core. As the shock propagates into the outer core with dissociating 
nuclei into free nucleons, the electron capture process e~ + p — > n + u e generates a 
huge amount of electron neutrinos just behind the shock. Before the shock arrives 
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Figure 5. The luminosities and root-mean-square energies of the neutrinos as a 
function of time. The results of the 13 M@ model are drawn with dashed lines and 
the results of the 40 M Q model with solid lines. The time is measure from the core 
bounce. This figure is taken from [199]. 



at the neutrino sphere, these electron neutrinos cannot escape in the hydrodynamical 
scale. Because these regions are opaque to the final state electron neutrinos and they 
are effectively trapped because their diffusion time is much longer than that for the 
shock propagation. As the shock waves move out in outer radius and pass through the 
neutrino sphere, the previously trapped electron neutrinos decouple from the matter 
and propagate ahead of the shock waves. This sudden liberation of electron neutrinos 
is called the neutronization burst (or "breakout" burst) (see the top panel of Figure 5). 
The duration of the neutronization burst is the timescale of the shock propagation and, 
hence, less than ~ 20 msec. While the peak luminosity exceeds 10 53 erg s -1 , the total 
energy emitted in the neutronization burst is only of the order of 10 51 erg due to the short 
duration timescale. This electron-neutrino breakout signal is expected to be detected 
from the Galactic supernova in modern neutrino detectors such as SuperKamiokande 
and the Sudbury Neutrino Observatory (see section 4). 

The breakout of the electron neutrinos is almost simultaneous with the appearance 
of the other neutrino species. In the hot post-bounce region, the electron degeneracy is 
not high and relativistic positrons are also created thermally leading the production of 
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anti-electron neutrinos (P' e s) via reaction e + + n — > p + V e . Mu- and tau- neutrinos are 
also produced in this epoch by the electron-positron annihilation (e + e~ — > ^ M)T P^ )T ), 
nucleon-nucleon bremsstrahlung (NN' — > NN' v^v^t), and neutrino/anti- neutrino 
annihilation (f e P e ~~ *■ ^t^Vt) ( see the bottom panel of Figure 5). Note that each process 
listed above also contributes for i/ e and V e neutrinos, the production of these neutrinos 
are predominantly determined by the charged- current interactions, e~ + p — > n + u e and 
e + + n — > p + z/ e . 

At a radius of 100 km ~ 200 km, the shock generated by core bounce stalls as a 
result of the following two effects. First, as the shock propagates outwards, it dissociates 
the infalling iron-peak nuclei into free nucleons, thus giving up ~ 8.8 MeV per nucleon in 
binding energy. Second, and most importantly, as the shock dissociates nuclei into free 
nucleons, electron capture on the newly-liberated protons to produce electron-neutrinos 
in the reaction e~ + p — > n + v e . 

Only for very special combinations of physical parameters, such as the stellar model 
of the progenitor or the incompressibility of nuclear matter, resulting in an extraordinary 
smaller cores, the so-called prompt explosion might work [130, 13, 131, 30], in which 
the shock wave at core bounce propagates through the outer core to produce explosions 
without the shock-stall (see Figure 6). 

2.3.5. delayed explosion Only several milliseconds after the shock-stall, a quasi- 
hydrostatic equilibrium is maintained between the newly-born protoneutron star (with 
radius of ~ 50 — 80 km) and the stalled accretion shock (with radius of ~ 100 — 200 km). 
The core of the protoneutron star is hot and dense, producing high-energy neutrinos of all 
species. If the energy transfer from the neutrinos to the material near the stalled shock is 
large enough, the stalled shock can be revived to produce the successful explosion. This 
neutrino-heating mechanism was discovered from the numerical simulations by Wilson 
[358] (see Figure 7). It is interesting to note that already in 1960's, Colgate and White 
proposed that neutrino heating was essential for producing the explosions [72]. The 
amount of gravitational binding energy {E zrscv ) released is huge, 

in contrast to the kinetic energy of canonical observed supernovae (.Ekm ~ 10 51 erg), 
where Mns and i?NS are the typical mass and the radius of a neutron star. Therefore 
in order to produce the explosions by the neutrino-heating mechanism, a small fraction 
(~ 1%) of the binding energy should be transfered, via neutrinos, to the mantle above 
the protoneutron star that is ejected as the supernova. 

For the better understanding of the mechanism, we give an order-of-magnitude 
estimation according to [35, 154]. Let us assume the situation that a neutrino sphere 
is formed at radius of R u and from there, the isotropic neutrino with luminosity of L u 
is emitted (see Figure 8). Then the neutrino heating rate of nucleons via the reactions, 
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Figure 6. Success or failure of the prompt explosion mechanism (courtesy of K. 
Sumivoshi and see also [302, 3031). Each panel shows the traicctorics of mass mesh 
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Figure 7. Successful of delayed explosion taken from [358] . The x and y axis represents 
the time in unit of second measured from core bounce and the radius from the stellar 
center in unit of ccnti-meter, respectively. Lines are trajectories of selected mass zones. 
The dashed line represents the shock front. 1.665 M© in the figure shows the mass 
point which is expelled outwards by the second shock due to the neutrino heating. As 
a result, the shock wave once weakened at ~ 500 msec revives and then successfully 
propagates to the surface of the iron core. 



n + v e — > e + p and p + u e ^ e + + n at a radius R (R u < R < R s ) can be estimated as, 
L v a{e v )Y N 
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ro( L » \( V( K Y' Z ( Yn ] \ 1 (29) 

V10 52 erg s- 1 / U5 MeV/ V 150 km/ Vl.O/ Ls • nucleonJ A ' 



10 52 erg s _ 

where L v is a typical neutrino luminosity, e u is the mean energy of neutrinos, a(e u ) is 
the cross section of the above absorption processes. Here we take Y N , the sum of the 
fraction of free nucleon and protons, to be 1 because nuclei are nearly fully dissociated 
into free nucleons after the passage of the shock waves. Outside the stalled shock, on 
the other hand, the above heating rates are suppressed because of the absence of the 
free nucleons. Note that each value assumed in the above estimation is taken from the 
recent ID core-collapse simulation [198]. On the other hand, the gravitational binding 
energy per baryon at a radius of R can be given as follows: 

where we take M^s to be a typical mass scale of 1.4M . Comparing the neutrino heating 
rate (r.h.s. of Eq. (29)) with the binding energy (r.h.s. of Eq. (30)), one can see that the 




Figure 8. Sketch of the stellar core during the shock revival phase. R v is the radius 
of neutrino sphere, from which neutrinos are emitted freely, R ns is the radius of the 
protoneutron star, R g is the radius (see text) and R s is the radius of the stalled shock. 
The shock expansion is impeded by mass infall to the shock front at R s by the mass 
infall at a rate M. This figure is taken from Janka (2001) [154]. 



neutrino heating can give the matter enough energy to be expelled from the core in 0.16 
second. In realistic situations, the cooling of matter occurs simultaneously via the very 
inverse process of the heating reactions, which delays the timescale of the explosion up 
to ~ 1 sec [358]. These timescales are much longer than those of the prompt explosion 
mechanism (O(10 msec)). Thus the neutrino-driven mechanism is sometimes called as 
the delayed explosion mechanism. 

Noteworthy, a characteristic radial position, which is the so-called gain radius, 
in which the neutrino heating and cooling balances and above which the neutrino 
heating dominates over the neutrino cooling, are formed after the shock-stall [33]. In 
the following, we estimate the position of the gain-radius by an order-of-magnitude 
estimation. In addition to the neutrino heating rate (Eq. (29)), the neutrino cooling 
rate of nucleons is given as follows: 

Q- = -a(T)acT\ (31) 

where T is the temperature of the material, cr(T) is the corresponding neutrino 
absorption cross section, a = 7/16 x 1.37- 10 26 erg cm -3 MeV 4 is the radiation density 
constant of neutrinos, and c is the speed of light. Since we assume for simplicity that 
the distribution function of neutrino is the fermi distribution with a vanishing chemical 
potential, then a'T 4 represents the energy density of neutrinos which yields to a black 
body radiation. Here we write L u in equation (29) as follows, 

L v = nRla'cT*, (32) 

assuming again that the neutrinos from the neutrino sphere are a Fermi distribution of 
the temperature T v of the neutrino sphere. Noting that a(T)/a{T v ) = (T/T u ) 2 , the net 
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heating rate can be written, 

Qtot — Qt ~ Qu 



Qt 



(fnr 



Using the simple power law relation, 



T 



T s — , 



(33) 



(34) 



which yields a good approximation in the radiation dominated atmosphere R v < R < i? s 
[154], the position of the gain radius R g becomes 



R„ 



(2R S 



Ry 



(35) 



Taking data obtained from a state-of-the-art numerical simulations [198], namely, 
T v = 4.8MeV,T s = lMeV, R s = 200km, R v = 80km, the gain radius becomes 85 km, 
which is in good agreement with the position numerically obtained by the corresponding 
simulations (see Figure 9). 
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Figure 9. Heating and cooling rates by electron (left) and anti-electron (right) 
neutrinos at 100 msec after core bounce taken from Liebendorefer et al. 2001 [198]. 
At the moment, the neutrino sphere defined by its mean energy and the stalled shock 
is located at R v ~ 80 km and R s ~ 200 km, respectively. In both panels, dashed 
lines with positive values, dashed lines with negative lines, and solid lines represent 
the heating rates, cooling rates, and the net rates, respectively. From the figures, it 
is shown that the location where the net heating changes sign defines the gain radius 
R g - 80 - 100 km. Lines labeled by "NR" or "GR" indicate that the results are 
obtained by the Newtonian or general relativistic simulations. 

The extent of the region of the net neutrino heating and the magnitude of the net 
neutrino energy deposition are responsible for producing the successful explosions and 
dependent crucially on the neutrino energy density and the flux outside the neutrino 
sphere, in which the neutrino semi-transparently couples to the matter. Thus the 
accurate treatment of neutrino transport is an important task in order to address the 
success or failure of the supernova explosions in numerical simulations. 
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If the neutrino heating mechanism works sufficiently to revive the stalled shock 
wave, the shock wave goes into the stellar envelope and finally blows off. This is observed 
as a supernova after the shock breaks out the photosphere. Unlike the case in the iron 
core, the photodissociation and the energy loss due to neutrinos is negligible in the 
stellar envelope and the binding energy is small, the shock wave successfully explodes 
the whole star. The propagation time of the shock wave depends on the stellar radius 
and is in the range from several hours to days. 

Here we shall mention that there is another type of supernovae, driven by a quite 
different physical mechanism. Supernovae Type la characterized by the absence of 
hydrogen lines in their spectra are thought to be caused by a thermonuclear explosion 
of a white dwarf that is completely disrupted in this event (for a review of SN la explosion 
models see [132]). Since the luminosities at the explosions of Type la supernovae are 
almost constant, they are good candidates to determine extragalactic distances and to 
measure the basic cosmological parameters. We will not consider them in this review. 
Supernovae we pay attention to in this thesis are the so-called Type II, Type lb and Ic, 
(for the recent observational classifications of supernovae, see [117] for example.) For 
convenience, we have used the common name "core-collapse supernovae" for supernovae 
of Types II and Ib/c. 
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3. Neutrino Oscillation 

In this section, we give a fundamental tool to discuss supernova neutrinos and their 
observation, neutrino oscillation. These topics have been seldom reviewed systematically 
so far. Starting from the physical foundation, we give an elaborate description of the 
neutrino oscillation, which has been established recently by a lot of experiments. Based 
on this, neutrino oscillations in supernovae are reviewed in the next section. 

3.1. Overview 

So far, we know three types of neutrino, v e , and u T . These are partners of the 
corresponding charged leptons, electron, muon and tauon, respectively, and produced 
via charged current interactions such as /3-decay and decays of muons and tauons. Thus, 
v e , and v T are called flavor (weak) eigenstates, which mean the eigenstates of the 
weak interaction. On the other hand, we can also consider eigenstates of their free 
Hamiltonian. They are called mass eigenstates denoted as v^i = 1,2,3) and have 
definite masses m^. 

These types of eigenstates come from essentially different physical concept so that 
they do not necessarily coincide with each other. In fact, this is the case with the 
quark sector: Flavor eigenstates are linear combinations of mass eigenstates which are 
determined by a unitary matrix called Cabbibo-Kobayashi-Maskawa matrix. Then, like 
the quark sector, it is natural to consider that the leptons are also mixing. 

The lepton mixing means that, for example, is e , which is produced by /9-decay, is 
a linear combination of some mass eigenstates i/j. More generally, neutrinos are always 
produced and detected in flavor eigenstates, which are not eigenstates of the propagation 
Hamiltonian. This mismatch leads to neutrino oscillation. 

Neutrino oscillation can be roughly understood as follows. Expressing the 
wave function of the neutrino by plane wave, each mass eigenstate evolves as 

— * — * 

exp [i(Eit — ki ■ #)], where Ei and ki are energy and momentum of the mass eigenstate 
z/j. Because different masses mj lead to different dispersion relations Ef = kf+mf, phase 
differences between the wave functions would appear as the neutrino evolves. Thus time- 
evolved wave function of a flavor eigenstate is no longer the original linear combination 
of mass eigenstates, which means that there is a probability that the neutrino is detected 
as a different flavor from the original flavor. 

3.2. Vacuum Oscillation 

Let us start with the Klein-Gordon equation neglecting the spin degree of freedom of 
neutrino, which is not important unless neutrino has large magnetic dipole moment. 
The equations of motion of the mass eigenstates in vacuum are 




(36) 
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where \J/( m ) = (^1,^2,^3) is a wave-function vector of the mass eigenstates and M = 
diag(mi, m 2 , m 3 ) is the mass matrix. Let us expand the wave function as 

¥ m \t,x) = e~ lEt ^\x), (37) 

where we assumed all the mass eigenstates have the same energy E. Although this 
assumption is not physically appropriate, the results below are not affected if the 
neutrino is ultra-relativistic. Then the equations of motion become 

If the neutrino is ultra-relativistic (E ~ ki ^> rrii), we have 

-E 2 - V 2 = —{E + iV)(E - iV) « -(£ + iV)2£, (39) 
which leads to 

= (40, 

where we set the direction of motion to z direction. The first term on the r.h.s. of Eq. 
(40) is irrelevant for neutrino oscillation because it just contributes to overall phase and 
will be neglected from now on. 

On the other hand, flavor eigenstates can be written as linear combinations of mass 
eigenstates as, 

H>g> = u¥™\ (41) 

where U is the mixing matrix, which is also referred to as the Maki-Nakagawa-Sakata 
(MNS) matrix [216]. This matrix corresponds to the Kobayashi-Maskawa matrix in the 
quark sector and often parameterized as, 



U 



( C12C13 sv2.cn s 13 e lS \ 

-C23S12 - Ci 2 S23Sl3e J<5 C12C23 - S 12 S2zSize l5 C13S23 
y s 12 s 23 — s 13 el<5 — c 12 s 23 — c 23 s 12 s 13 el<5 C13C23 J 



(42) 



where s^- = sin^- and = cosOy, 9ij(ij = 1,2,3) are mixing angles and 5 is CP 
phase. In terms of the flavor eigenstates, the equations of motion (40) are expressed as, 

= -5b-*SP- («) 

Here it should be noted that the mass matrix for flavor eigenstates, UMW, is not 
diagonal in general. 

As a simple example, let us assume that there are only two neutrino species, v e and 
v^. Then the mixing matrix can be written as, 

( cos6 sine \ 

y - sine cose J y ' 

and the equations of motion reduce to 
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for mass eigenstates and 
d 



-cos 29 


sin 29 \ 




sin 26 


cos 29 ) 


:: 



for flavor eigenstates. Here Am 2 = m\ — m{ and we again neglected a term proportional 
to identity matrix. Then consider a neutrino which is purely electron-type at first. 
Noting that electron-type neutrino can be written in terms of mass eigenstates as, 

\v e ) = cos^i) +sm9\v 2 ), (47) 

the neutrino evolves according to, 

(in 2 \ ( tyi 2 \ 

~^^2E Z J cos @Wi) + ex P ( ~^~2E Z ) sm ^l^)- (48) 

Multiplying (i/ e \, we obtain the probability that this state is electron type, 



P Ve ^ v Sz) = \(veHz)}\ 2 = 1 - sin 2 29 sin 2 (tt-^- 1, (49) 



*-osc 



where 



<- s = 2 48 x 10 cm ( mev) (-ST ) (50) 

is called the oscillation length. It is easy to show that 

P v ^ v Xz) = P v ^{z) (51) 
Pu e ^ e (z) = P n (z) = l- P„.^,„ (z) (52) 

as expected by unitarity. 

The probability P Ve ^ v ( function of propagation distance z is plotted in Fig. 

10. It oscillates with respect to z and the wave length is the oscillation length £ osc . 
Here it will be worth noting that the oscillation length depends on the neutrino energy 
and the mass difference of the two involved mass eigenstates as is seen in (50). The 
amplitude is determined by the mixing angle and is the largest when 9 = ir/4. Thus 
even if a neutrino is produced as an electron-type neutrino, there is non-zero probability 
that it is detected as a muon-type neutrino if there is a mixing between the two neutrino 
flavors. It is this phenomenon which is known as the neutrino oscillation. 

Let us consider a more general case with many neutrino species. A neutrino state 
with a flavor a can be written as a linear combination of the mass eigenstates, 

kHEOO- (53) 

i 

Then the evolution of a neutrino which is originally u a is 

i ^ ' 

and the probability P Va ^ Vf3 {z) is 

(Am 2 \ 
- i -2E~ z ) ( 55 ) 

where Am 2 - = m 2 — m 2 . 
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osc 



Figure 10. Conversion probability P Ve ^ v (z) as a function of propagation distance z 
normalized by the oscillation length £ osc 



averages Finally we consider two averages of the conversion probability concerned with 
the finite size of a source and a finite energy width. 

A neutrino source is, in general, not point-like and has a finite size. For example, 
in the sun, there is a spherical neutrino source with a radius ~ 10 10 cm [20]. In this case, 
the finite size will average the phase of the oscillation of the conversion probability. 
Denoting the source distribution as f(z ), the conversion probability is given by 

2 ( n(z-z o y 



Pu e ^(z) = sin 29 / dz f(z )sm 



(56) 



where z is the distance between the center of the source and the detector. If we assume 
a Gaussian distribution with a width of s, that is, 

1 



/(*b) 



the conversion probability is computed as, 



P,^^)^ sin 2 20 



osc COS 



2nz\ 

^osc / 



(57) 



(58) 



and is shown in Fig. 11. As one can see, if a source has a finite size, the amplitude of 
the probability oscillation become small while the average value remains unchanged. In 
the limit of s — > 0, it reduces to Eq. (49) and the oscillation is completely smoothed for 

S 3> (osc- 

Next, let us consider a finite energy width. If a source has a finite energy width 
we must also average the conversion probability by the energy spectrum because the 
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Figure 11. Conversion probability P Ve ^ v ^(z) averaged by taking the finite size of a 
source into account. 



5=0 5: 


=A)/10 5»A) 

MA: 











Z/lo 



Figure 12. Conversion probability P„ e ^ v (z) averaged by a Gaussian energy 
distribution. 



oscillation length (50) depends on neutrino energy: 
Pu.^(z) = sin 2 26 J dEg(E) sin 2 



2tt 
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(59) 
(60) 



losc(E) 2E 

where g(E) is the energy spectrum of neutrinos. As a simple example, we consider a 
spectrum with Gaussian A(E): 



9(E) 



2S 1 



A = 



2tt 
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(61) 



5y/2n " k{E) 2E ' 

where 5 is the width and E is the central energy. Here it should be noted that this 
distribution is not Gaussian with respect to neutrino energy. Then we have an averaged 
conversion probability, 



- sin 2 26 
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(62) 
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This is plotted in Fig. 12. Although the oscillation is smoothed as in the case of the 
finite-size source, the behavior is different between the two cases. Since difference in 
energy leads to difference in oscillation length, the phase difference of two neutrinos with 
different energies increases as they propagate a long distance. Therefore, the conversion 
probability will cease to oscillate in the end regardless of the magnitude of the energy 
width. 



3.3. Oscillation in Matter 

The behavior of the neutrino oscillation changes if the neutrino propagates in the 
presence of matter, not in vacuum. Due to the interaction with matter, neutrino 
gains effective mass, which modifies the dispersion relation. If the interaction is flavor- 
dependent, like that with electrons, the change of the dispersion relation is also flavor- 
dependent. Remembering that the neutrino oscillation in vacuum is induced by different 
dispersion relations due to different masses, it is easy to imagine that further change in 
dispersion relations will change the behavior of the neutrino oscillation. This effect, the 
MSW effect, was first pointed out by Wolfenstein [353, 354] and discussed in detail 
by Mikheyev and Smirnov [226, 227, 228]. As will be discussed below, if matter 
is homogeneous, the situation is essentially the same as the vacuum oscillation with 
effective mixing angles determined by matter density and the original mixing angles. 

What is interesting and important is a case with varying density. In fact, the MSW 
effect with varying density gave the solution to the long-standing solar neutrino problem 
[314, 298] and will also play a important role in supernovae. 

3.3.1. constant density At low energies only the elastic forward scattering is important 
and it can be described by the refraction index n re f. In terms of the forward scattering 
amplitude f(E), the refraction index is written as 

n ref = l + ^/(£), (63) 
where n is the target density. Then we have the dispersion relation in matter as 

(n rei E) 2 = k 2 + m 2 . (64) 
If we rewrite this dispersion relation as 

(E-V cS ) 2 = k 2 + m 2 , (65) 
we obtain the effective potential V e s as 
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V e s = (1 — n re f)E = — —f(E). (66) 

On the other hand, low-energy effective Hamiltonian for weak interaction between 
a neutrino and a target fermion is 

H int = ^ n ,(C v - C Al ^ s ^{\ - 75 )^ (67) 
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Table 1. Effective weak coupling constant for neutral-current interactions. Here #w 
is the Weinberg angle sin 2 6w ~ 0.23. 



where ipf is the target fermion field, ip u is the neutrino. Here the coupling constant Gp 
is the Fermi constant, 

Gp = 1.2 x l(T 5 GeV~ 2 = 9 x l(T 38 eV cm 3 , (68) 

and Cy and are vector weak charge and axial- vector weak charge, respectively, which 
depend on the species of the target. For the neutral-current interactions, charges are 
shown in Table 1 and for the charged- current interaction, Cy = Ca = 1- 

Using the Hamiltonian and coupling constant, the forward scattering amplitude, 
refraction index and effective potential are computed as, 

f(E) = T C v G F -4r, (69) 
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(70) 
(71) 



for neutrino and anti-neutrino, respectively. Here ub is the baryon density, rif and nj 
are fermion and anti-fermion number density, Yf = (rif — nfj/ns is the fermion number 
fraction per baryon and 
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Assuming charge neutrality (Y e = Y p ), we have 



V eS = ±V2G F n B x 
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With this effective potential, the wave equation (43) is modified as 
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where A is the mass matrix representing the contribution from interactions with matter. 
Neglecting the background neutrino, we have 

{ 3Y e -l \ 

Y e - 1 . (75) 
Y e -1 J 

Here we used Y n = 1 — Y p = 1 — Y e but the contribution from neutrons is not important 
for neutrino oscillation because it is proportional to identity matrix. This reflects the 
fact that the interaction with neutrons is via the neutral-current interaction which occurs 
equally to all flavors. 

Again, let us consider the two-flavor case. The wave equation in matter (74) can 
be rewritten as 

Am 2 




4E 



— cos 29 m 


sin20 m \ 




sin 28 ' m 


cos 28 m J 


:: 



(76) 



up to terms proportional to identity matrix. This is exactly the same form as the 
vacuum case (46) with modified parameters, 

sin20 m = - 77 
v/(£-cos2#) 2 + sin 2 2# 



Am 2 m = Am 2 y / (£-cos2#) 2 + sin 2 2#, (78) 
where £ is the dimensionless density parameter: 
_ 2V2G F n B E 
Am 2 

— 2 (t^)(t^)(^). 

The oscillation length in matter is also defined in the same way, 
_ AttE _ sin2fl £ osc 
° SC ' m "A ffl ^ sin2^ osc cos 2^ + sin 2 29 1 J 

Thus neutrino oscillation occurs with modified mixing angle 6> m and oscillation length 
4>sc, m - Mass eigenvalues can be obtained by diagonalizing (74) as 



2 m\ + m\ Am 2 



m m = h 



1) 



(2y e -l)£± A /(£-cos20) 2 + sin 2 20 

In Fig. 13, behaviors of various parameters in matter as functions of the dimensionless 

density parameter £ are shown. 

When £ = cos 26, the mixing angle in matter, 9 m , becomes maximum (7r/4). This 

is called resonance and the resonance condition can be rewritten as 

2 s n /^0.5\ ZlMeV\ / Am 2 \ . , 

Prcs = 1.3 x 10 2 g cm- 3 cos 29 (-j J (^j . (82) 

At the resonance density, the matter oscillation length and the mass-squared difference 
become maximum and minimum, respectively: 

(/osc,m)rcs = , (Ai^) ra = Am 2 sin 26*. (83) 

sin ^t/ 
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Figure 13. Behaviors of various parameters in matter as functions of the dimensionless 
density parameter £. The mixing angle in vacuum is set as sin 2 29 = 0.1. 

On the other hand, in the case of anti-neutrino, the signature of the matter effect is 
different so that there is no resonance. However, we will discuss possible resonance in 
anti-neutrino sector later. 



3.3.2. varying density The solar neutrinos are produced at the center with the density 
about 150g cm" 3 , and escape outward into vacuum. In this case and many more 
cases in astrophysical systems including supernovae, the neutrinos propagates in an 
inhomogeneous medium and the neutrino oscillation becomes much more complicated. 
As we saw in section 3.3.1, flavor eigenstates and mass eigenstates can be related by the 
effective mixing angles in matter as 



(84) 



In an inhomogeneous matter, the mixing angles 9i >m are functions of z. Due to this 
dependence of the mixing angles on z, the wave equations cannot be solved analytically 
in general. 

To see this, consider a two-flavor case. The wave equations were given in (76): 

Ami, 







f 




\ 




= U{0 i>m ) 








\ v - ) 
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^3,m 


/ 



.d_ 

dz 



4E 



— cos 29 m 


sin20 m \ 




sin 26 m 


cos 29 m J 





(85) 



If the mixing matrix U is constant, we can diagonalize the equations by multiplying U 
to the both sides of (85). However, if the mixing matrix U depends on z, the derivative 
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operator d/dz and U do not commute so that the l.h.s. does not result in a simple form, 
although the r.h.s. is diagonalized: 



d 



dz 



l,m 



Thus the equations are not effectively diagonalized and can be written as, 




m 









2oj 



Ami 




(87) 

9z 2a; 

This shows that even the mass eigenstates are mixing if the density is inhomogeneous 
and the magnitude of the mixing depends on d9 m /dz. 

Let us first consider the mixing of the mass eigenstates qualitatively. If the diagonal 
component is much larger than the off-diagonal component everywhere, that is, 



cm, 



dz 



Ami 
~2uo 



(88) 



the mass eigenstates will propagate without mixing. In other words, the heavier state 
remains heavier and the lighter state remains lighter. There is no energy jump and 
this case can be said to be "adiabatic". Contrastingly, if the condition (88) is not 
satisfied, the heavier state can change to the lighter state and vice versa, which is a 
" non-adiabatic" case. The non-adiabaticity is largest where the change of the mixing 
angle is rapid, which is expected to be around the resonance point as can be expected 
from Fig. 13. 
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Figure 14. Schematic view of two-flavor resonance in an inhomogeneous medium. 
The dashed lines show the masses of the mass eigenstates and the solid lines show 
the effective masses of the flavor eigenstates. Contribution from the neutral-current 
interactions are subtracted. 

Fig. 14 will be helpful to understand the oscillation in an inhomogeneous matter. 
Assume that the vacuum mixing angle is so small that the lighter state is almost v e , that 
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is, v e is effectively " ligher" than in vacuum. In a dense region, on the other hand, v e 
is effectively "heavier" than v^. The two flavors have the same "mass" at the resonance 
point. Then let us consider a case where a v e is produced at a dense region and escape 
into vacuum as the solar neutrinos. If the resonance is adiabatic, the heavier state will 
remain heavier, which means that a v e produced at the dense region will emerge as a 
u^. In contrast, if the resonance is non-adiabatic, a u e will remain a v e . 

Thus, the survival probability of u e depends largely on the adiabaticity of the 
resonance. The importance of the resonance in the solar neutrino problem was first 
pointed out by Mikheyev and Smirnov [226]. It was Bethe who found the essence of 
the MSW effect in an inhomogeneous medium to be the level crossing of the flavor 
eigenstates [34]. 

Let us discuss more quantitatively. The adiabatic condition (88) at the resonance 
point can be rewritten as, 

1 < \ d lnn e \ Sil1 29 taI1 26 ■ ( 89 ) 

E — o — - 

I oz I res 

The adiabaticity parameter 7 is defined as 

7 = m sin 29 tan 29, (90) 

I dz I res 

so that the adiabaticity condition is, 

7>1- (91) 

What we want to know is the conversion probability of the matter eigenstates, Pu 1:m ^u 2 , m - 
For a general profile of matter density, this cannot be obtained analytically. But for 
some special cases, analytic expression is known [182]. If we write the probability as, 



exp [-ffi-e xp [-- '-- 
l-exp[- 2snr 



p _ ^' 2ik??J /qo\ 



the factor F is given by, 

{1 (n e oc z) 

feSSf (n e ^z--) . (93) 

l-tan 2 # (n e oce- 2 ^ ) 

Note that P Ul . in ->u2, m ~ 1 when 7 ^> 1 and P V1 -n^ ~ when 7 <^ 1 in any cases, as 
expected. 

If we obtain the conversion probability P Ul . m ^u2, m m some way, we can compute, 
for example, the survival probability of u e , P Ue ^ Ue . When a v e is produced at the center 
of the sun, the probabilities that it is v\ )Xa and z/ 2im are cos 2 9 C and sin 2 9 C , respectively, 
where 9 C is the mixing angle at the center. First, if there is no conversion between v\ )Xa 
and i/ 2 ,m, that is, if P vl ^ V2<m = 0, 

Pv e -n; e = cos 2 9 C cos 2 9 + sin 2 9 C sin 2 9 
1 + cos 29 c cos 29 
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On the other hand, when P v , _>,,„ is non-zero, 



1 -P) cos 2 ^H-Psm 2 9 C 



cos 2 9 + 



Pcos 2 9 C + (1-P) sin 9 C 



cos 29 c cos 29 



Pcos 29 c cos 20, 



sin 2 ^ 
(95) 



where P = P, lim _^ 2 , m . 

In most cases, the survival probability is determined by the adiabaticity parameter 
7. Because 7 depends on the neutrino parameters, mixing angles and mass differences, 
and density profile of matter, neutrino observation from various systems will allow us 
to investigate them. 

If the matter density is much larger than the solar case, the two-flavor analysis 
is invalid and we have to take three flavors into account. This is exactly what we do 
later to consider neutrino oscillation in supernova. For a three flavor case, there are two 
resonance points. Although the situation will become more complicated, the essence is 
the same as the two-flavor case, the adiabaticity of the resonance. 



3.4- Experiment of Neutrino Oscillation 

Because the neutrino oscillation is a phenomenon beyond the standard model of particle 
physics, many experiments have been conducted to verify it in various systems. One 
of the attractive features of neutrino oscillation experiment is that it does not need 
high-energy accelerator. 

In this section, we review neutrino oscillation experiments starting from general 
remarks about the experiments. 



3.4-1- general remarks The basic of neutrino oscillation experiment is to observe 
neutrinos from a known source. Although all flavors except sterile neutrino can be 
ideally detected, v e and v e are easier to detect than other flavors so that they are often 
used as signal. In this respect, neutrino oscillation experiment can be classified into 
two types. One is called "appearance experiment", in which, for example, we detect 
neutrinos from a source. If we observe even a single event of v e , this is an evidence 
of flavor conversion. Another type is called "disappearance experiment", in which we 
observe u e s (P e s) from a v e (P e ) source with known flux. If we observe less number of 
is e s, this can also be an evidence of flavor conversion. 

Each type has its own advantage and disadvantage. By the appearance experiment, 
we cannot reject neutrino oscillation phenomenon even if we did not observe the signal. 
This is because might have changed into v T , not u e . Contrastingly, the disappearance 
experiment can tell whether neutrino oscillation occurred or not, if only z/ e s were 
converted into any type of neutrino. However, we cannot know the oscillation channel, 
that is, which flavor z/ e s were converted to. In this respect, the appearance experiment 
can probe a selected channel of neutrino oscillation. 

Many kinds of experiments have been done so far and each has different parameter 
region it can probe. Here we discuss how the neutrino parameters can be probed. More 
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Table 2. Probe of squared mass difference by various experiments. 



source energy E(MeV) baseline L(m) Am 2 (eV 2 ) 



accelerator 


10 3 ~ 10 5 


10 2 


~ 10 6 


io- 3 ~ 10 2 


reactor 


1 ~ 2 


10 


~ 10 5 


10~ 5 ~ 10" 1 


atmosphere 


~ 10 3 


10 5 


~ 10 7 


io~ 2 ~ 10~ 4 


Sun 


~ 1 




10 11 


~ io- 11 



concretely, let us consider what we could know if there was no signal in an appearance 
experiment. For simplicity, we consider just two-flavor oscillation in vacuum. As we saw 
in section 3.2, conversion probability P u ^u e as a function of the propagation distance 
is, 

P^ e (^) = sin 2 2£sin 2 (f^) (96) 

\ 'osc / 

If we detect no v e signal, it means that the conversion probability is smaller than a 
certain value 5 which is determined by the noise level of the experiment. When the 
baseline L is much smaller than the oscillation length £ osc , it is written as, 



nL\ 2 



P^ e (L) » sin 2 29 — <5. (97) 

\ ^osc / 

Substituting the definition of £ osc (50), this reduces to 

Am 2 sin2#<^. (98) 

On the other hand, when L ^> £ osc , finite energy width of the neutrino beam will average 
the oscillation of the conversion probability so that the no signal means, 

P^ e (L) = ^sin 2 20<5, (99) 

which reduces to, 

sin2^<v / 27. (100) 

Note that we cannot obtain information about Am 2 in this case. 

Thus if we want to probe small Am 2 , experiments with small E/L are advantageous. 
Various systems with characteristic neutrino energy, baseline and possible Am 2 which 
can be probed are shown in Table 2. Analysis of an disappearance experiment can be 
done essentially in the same way. 

3.4-2. accelerator experiment Accelerator experiment is the most popular experiment 
of neutrino oscillation. One of the advantage of accelerator experiment is that we can 
control the neutrino source while Sun, atmosphere and supernovae are uncontrolled and 
rather unknown sources. Although there have been a lot of accelerator experiments so 
far, the basic concept is similar as we will review below. 
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Table 3. Accelerator experiments 



experiment 


baseline 


neutrino energy 


detector 


status 


reference 


CCFR 


~ 1km 


30 - 500GeV 


690 ton target calorimeter 


completed 


[63, 64] 


NuTeV 


1.4km 


~ lOOGeV 


690 ton target calorimeter 


completed 


[63, 255] 


KARMEN 


17.5m 


~ 50MeV 


56ton liquid scintillator 


completed 


[167, 168] 


LSND 


30m 


~ 50MeV 


167ton liquid scintillator 


completed 


[207, 208, 209] 


NOMAD 


625m 


~ 50GeV 


2.7 ton drift chambers 


completed 


[253, 254] 


MiniBooNE 


490m 


0.1 ~ 2GeV 


445ton mineral oil 


ongoing 


[230, 231] 


K2K 


250km 


~ lGeV 


22500ton water (SK) 


ongoing 


[158, 159, 160] 
and [161] 



First, protons are accelerated and collided with target nuclei to produce tt ± : 

P + N ^7T + , (101) 

Then 7r + s or n~s are absorbed and the others decay to produce neutrinos. For example, 

7T+^/X + + ^, (102) 

H + -> e+ + u e + P„ (103) 

In this way, if 7r~s are absorbed completely, we have a neutrino beam which consists of 
v e , and v^. Thus, if we detect u e s in this beam, we can confirm neutrino oscillation 
— > u e . In fact, 7r~s can not be absorbed completely and decay like, 

7i- + Vp, (104) 

IT -> + v e + iy (105) 

Consequently, some u e s will be produced and they become one of main noises in this 
kind of experiment. 

In this experiment, we can control the amount of the produced neutrinos. If there 
are v e events we can compute the conversion probability Pp M _»p e . Or, we can obtain an 
upper limit for the conversion probability if no signal was obtained. Anyway, we can 
extract information about neutrino parameters as we saw in section 3.4.1. 

The current and past accelerator experiments are classified into long-baseline 
experiments and short-baseline experiments. Of course, the latter is technically easier so 
that early experiments have rather short baselines. However, implication from the recent 
solar and atmospheric neutrino observations has made it necessary to probe very small 
Am 2 by long-baseline experiments. Features of some of the current and past accelerators 
are shown in Table 3. Among them, CCFR, NuTeV, KARMEN and NOMAD had no 
signal for neutrino oscillation and obtained constraints on oscillation parameters which 
are shown in Fig. 15. 

The LSND experiment performed at Los Alamos observed 87.9 ± 22.4 ± 6.0 excess 
events in the — > V e appearance channel [209]. This signal corresponds to a transition 



39 




3V - 



Figure 15. Parame- 
ter regions allowed by 
the LSND observation 
(purple, blue). Exclu- 
sion curves from vari- 
ous other experiments 
are shown as well, with 
the region on the right 
side excluded [254]. 
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Figure 16. Ex- 
pected excluded re- 
gions in case of a non- 
observation of signal in 
MiniBooNE [231]. 




Figure 17. The recon- 
structed energy spectrum for 
the ^-like sample. Points 
with error bars are data. The 
solid line is the best fit spec- 
trum and the dashed line is 
the expected spectrum with- 
out oscillation [161]. 
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Figure 18. Allowed regions 
of oscillation parameters for 
— ► v r . Dashed, solid and 
dot-dashed lines are 68.4 %, 
90 % and 99 % CL. contours, 
respectively [161]. 



probability of P = (0.264 ± 0.067 ± 0.045)%, which is ~ 3.3a away from zero. If we 
interpret it in terms of two-flavor oscillation, parameter regions shown in Fig. 15 are 
allowed. Although most of the allowed regions are excluded by other experiments, there 
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Table 4. Reactor experiments 

experiment baseline status reference 

Bugey (France) 15,40,95m completed [1] 

CHOOZ (France) 1km completed [67, 68] 

Palo Verde (USA) 750m completed [39, 40] 

KamLAND (Japan) 100km ~ running [164, 165, 166] 



are still some surviving regions with Amlsnd ~ 0.1 — leV 2 . This remaining regions are 
expected to be confirmed or denied by the MiniBooNE experiment [230, 231] in the near 
future. Fig. 16 shows expected excluded regions in case of a non-observation of signal 
in MiniBooNE. 

The K2K experiment [158, 159, 160, 161], the KEK to Kamioka long-baseline 
neutrino oscillation experiment, is an accelerator based project with 250 km baseline 
which is much longer than those of the past experiments. This long baseline make 
it possible to explore neutrino oscillation in the same Am 2 region as atmospheric 
neutrinos. In [161], five-year data with 57 candidates was reported and their energy 
distribution is shown in Fig. 17. The distortion of energy spectrum, which signals 
neutrino oscillation, is clearly seen in Fig. 17 and the probability that the result would 
be observed without neutrino oscillation is 0.0050%(4.0<r). Fig. 18 shows a two-flavor 
neutrino oscillation analysis with disappearance. The best fit point is, 

sin 2 29 = 1.0, Am 2 = 2.8 x 10~ 3 eV 2 . (106) 

3.4-3. reactor experiment In nuclear reactors, u e s are isotropically emitted by /3-decay 
of neutron-rich nuclei. Reactor experiment of neutrino oscillation is detecting these z/ e s 
and seeing if there is a deficit compared with the expected flux. The flux and spectrum 
of v e are determined by the power of the reactor and abundance of 235 U, 238 U, 239 Pu 
and 241 Pu. Because reactor neutrinos have relatively low energy, they are well suited 
in exploring the region of small Am 2 at modest baselines. For example, to explore 
the parameter Am 2 down to 10~ 3 eV 2 a reactor experiment with energy around 5 MeV 
requires a baseline of L = 1 km, while an accelerator experiment with E = 5 GeV would 
require L — 1, 000 km. 

The current and past reactor experiments are shown in Table 4. Among these, 
CHOOZ [68] gives the strongest constraint on the mixing angle of u e — > v x for 
Am 2 > 10~ 3 eV 2 (Fig. 19). 

Inspired by the recent development of solar neutrino observation, a long-base line 
experiment called KamLAND [164, 165, 166] was constructed where there was once the 
Kamiokande detector. KamLAND consists of lOOOton liquid scintillator and its primary 
purpose is to confirm the solution to the solar neutrino problem, which will be discussed 
later. There are 16 reactors with distances about 100 ~ 1000km from KamLAND so 
that the v e flux at KamLAND is sufficiently large to probe neutrino oscillation Because 
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Figure 19. Con- 
straint to conversion 
v e -► P x from CHOOZ 
experiment [68] and 
that of conversion 
— > v e from K2K 
experiment [160]. 
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Figure 20. Prompt 
event energy spectrum 
of v e candidate events 
with associated back- 
ground spectra. The 
shaded band indicates 
the systematic error 
in the best-fit reac- 
tor spectrum above 2.6 
MeV [166]. 
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Figure 21. (a) Neutrino oscillation parameter allowed region from KamLAND anti- 
neutrino data (shaded regions) and solar neutrino experiments (lines) which we discuss 
later [166]. (b) Result of a combined two-neutrino oscillation analysis of KamLAND 
and the observed solar neutrino fluxes. 



of the long baseline and detectability of low-energy neutrinos, KamLAND can probe 
much smaller Am 2 compared with the past experiments. 

Fig. 20 shows the prompt event energy spectrum of v e candidate events with 
associated background spectra [166]. They observed 258 v e candidate events with 
energies above 3.4 MeV compared to 365.2 events expected in the absence of neutrino 
oscillation. Accounting for 17.8 expected background events, the statistical significance 
for reactor v e disappearance is 99.998%. Also the observed energy spectrum disagrees 
with the expected spectral shape in the absence of neutrino oscillation at 99.6% 
significance and prefers the distortion expected from v e oscillation effects rather than 
those from neutrino decay and decoherence. A two-neutrino oscillation analysis of the 



42 

KamLAND data gives 

Am 2 = 7.9+°;^ x l(T 5 eV 2 , (107) 

as shown in Fig. 21. A global analysis of data from KamLAND and solar neutrino 
experiments, which will be discussed later, yields 

Am 2 = 7.9t°oi x 10" 5 eV 2 , tan 2 9 = 0.40t°;$. (108) 
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Figure 22. The zenith angle distribution for fully-contained 1-ring events [315]. The 
points show the data, box histograms show the non-oscillated Monte Carlo events and 
the lines show the best- fit expectations for <-> v T oscillations with sin 2 29 = 1.00 
and Am 2 = 2.1 x 10~ 3 cV 2 . The height of the boxes shows the statistical error of the 
Monte Carlo. 

3.4-4- atmospheric neutrino When cosmic rays enter the atmosphere, they collide with 
the atmospheric nuclei to produce a lot of mesons (mostly pions): 

iVcR + A^^T^Tr ,---. (109) 

The pions decay and the decay product, such muons, further decay to produce neutrinos: 
vr ± -> fj^ + Vpipp), -> e ± + v e {y e ) + z^(z^) (110) 
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When cosmic rays enter the atmosphere, they collide with the atmospheric nuclei 
to produce a lot of mesons (mostly pions): 

iVcR + iVair-^TfV--. (HI) 

The pions decay and the decay product, such as muons, further decay to produce 
neutrinos: 

7r ± ->^ ± + i/ M (P /i ), ^ -> e ± + v e (V e ) + v^Vp) (112) 

These neutrinos are called the atmospheric neutrinos and they provided the first strong 
indication for neutrino oscillation. Their energy range and path length varies from 0.1 
to 10 GeV and from 10 to 10,000 km, respectively, which indicates that atmospheric 
neutrinos can provide an opportunity for oscillation studies over a wide range of energies 
and distances. From (112), we simply expect that the flux ratio (z/ M + P M ) / [y e + 9 e ) = 2. 
This is roughly correct, even though the ratio depends on a lot of factors such as 
neutrino energy and zenith angle of the incoming neutrino if detailed decay processes 
and geometric effects are taken into account [29, 139, 140]. However, the results from 
many experiments showed that this ratio was about unity. This was once called the 
atmospheric neutrino problem, which is now interpreted successfully in terms of neutrino 
oscillation. In fact, the SK was the first to prove neutrino oscillation phenomenon by 
its observation of the atmospheric neutrino with a high accuracy. 

It was reliable identification of that allowed SK to prove neutrino oscillation. 
Because the SK detector is huge (diameter ~ 39m and height ~ 42m), z/^s with energy 
~ 1 GeV can be identified as fully contained events, for which all neutrino-induced 
interactions occur in the detector and neutrino energies and directions can be accurately 
obtained. The results from SK are shown in Fig. 22. As can be seen, v e flux is consistent 
with the theoretical calculation while there is a deficit in flux. It indicates that 
<-> v T oscillation would be solution to the atmospheric neutrino problem. In Fig. 23, 
allowed oscillation parameters for <-> v T oscillations are shown. The best fit values 
are 

Am^ tm = 2.1 x 10" 3 eV 2 , sin 2 # atm = 1.00. (113) 

These results are confirmed by other experiments such as Soudan 2 [316] and MACRO 
[214, 215] (see also [109]). 

As we saw in Eq. (49), the conversion and survival probabilities depends on L/E, 
where L and E are neutrino path length and energy, and oscillates with respect to L/E. 
This behavior was confirmed by the atmospheric neutrino observation at SK [313]. This 
means that it rejected other possibilities such as neutrino decay and decoherence which 
had different dependence on L/E and proved that solution to the atmospheric neutrino 
problem was truly neutrino oscillation. 

3.4-5. solar neutrino In the central region of the sun, z/ e s are continuously produced 
by the hydrogen fusion, 

4p -> 4 He + 2e + + 2u e + 26.731MeV. (114) 
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Figure 23. Allowed oscillation parameters for <-» v T oscillations from the 
atmospheric neutrino observation at SK [315]. Three contours correspond to the 68% 
(dotted line), 90% (solid line) and 99% (dashed line) C.L. allowed regions. 



These neutrinos are called solar neutrinos. The first experiment of the solar neutrino 
observation was the chlorine experiment at Homestake by R. Davis and his collaborators 
in the 1960's [74, 71]. On the other hand, the theoretical study of the solar neutrino was 
pioneered by J. N. Bahcall. The neutrino production rate in the sun has been calculated 
based on the standard solar model, which reproduces the current state of the sun by 
following the evolution of a main-sequence star with solar metalicity and mass. The 
solar neutrino flux calculated by the current standard solar model [22, 24] is shown in 
Fig. 24. 

However, the fluxes observed by several detectors have been substantially smaller 
than the predicted flux (Fig. 25). This is so called solar neutrino problem and has 
been studied for several decades (for reviews, see [114, 21], and a text book by Bahcall 
[20]). Historically, the solar neutrino problem was attributed to incompleteness of the 
solar model and/or unknown neutrino property including neutrino oscillation. With the 
improvement of both the solar model and observation, especially at SuperKamiokande 
[307, 310, 311, 314], it is now common to think that the uncertainties of the solar model 
cannot solely explain the gap between the observations and prediction. 

The critical observation has been conducted by Sudbury Neutrino Observatory 
(SNO) [295, 296, 297, 298], which clearly showed that electron neutrinos are converted 
to the other flavors. The SNO was designed primarily to search for a clear indication of 
neutrino flavor conversion for solar neutrinos without relying on solar model calculations. 
Its significant feature is the use of 1000 tons of heavy water which allows the distinction 
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Figure 24. Solar neutrino flux from various nuclear reactions [241. 




Figure 25. Observed and predicted fluxes of solar neutrino [19]. 



between the following three signal: 

u e + d — > p + p + e", (115) 
v + d^p + n + v, (116) 
z/ + e~^z/ + e~. (117) 

Here the reaction (115) occurs through the charged current interaction and is relevant 
only to u e , while the other two reactions are through the neutral current interaction 
and sensitive to all flavors. It should be noted that SK can identify only the electron 
scattering event (117) for energies of the solar neutrinos (< lOMeV), although its volume 
is much larger than that of SNO. 

Fig. 26 shows the fluxes of fi + r neutrinos and electron neutrinos obtained from 
SNO [298] and SK [312]. Combining the signals from the three channels, it clearly shows 
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that there is non-zero flux of v ^ and v T , which is a strong evidence of neutrino flavor 
conversion. If we interpret these data by 2-flavor neutrino oscillation v e <-> u x , we obtain 
constraint on the mixing angle and mass-squared difference as in Fig. 27 with the best 
fit: 

AmLar = 6.5jtl x 10" 5 eV 2 , tan 2 # solar = 0.45+°;^. (118) 

Also it should be noted that, as can be seen in Fig. 26, the prediction of the 
standard solar model is, at least concerning the high-energy neutrinos, confirmed by 
the observation. 




$ e (x lCTcm 2 s 1 ) 

Figure 26. Flux of fx + 
t neutrinos versus flux of 
electron neutrinos obtained 
from SNO [298] and SK 
[312]. The total 8 Be solar 
neutrino flux predicted by the 
Standard Solar Model [24] is 
shown as dashed lines. The 
point represents <f> e from the 
CC flux and (j)^ T from the 
NC-CC difference with 68%, 
95%, and 99% C.L. contours 
included. 
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Figure 27. Global 
neutrino oscillation 
analysis using solar 
neutrino data [298]. 
The solar neutrino 
data includes SNO 
data, the rate mea- 
surements from the CI, 
SAGE, Gallex/GNO, 
and SK-I zenith 
spectra. 



3.4-6. current status of neutrino parameters Here we summarize the information of 
neutrino oscillation parameters obtained so far. Combined analysis based on all the 
data from various neutrino oscillation experiments has been done by many authors 
[92, 113, 217, 218, 300]. Basically three-flavor neutrino oscillation can explain all 
the data reasonably except the LSND data discussed in section 3.4.2. Maltoni et 
al. [218] performed a general fit to the global data in the five-dimensional parameter 
space (9i2 , 2 3 , 0i 3 , Aml 2 and Am 2 3 ), and showed projections onto various one- or two- 
dimensional subspaces. They gave, 

• sin 2 #i2 ~ 0.30, mostly from solar neutrino data 
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Table 5. Best-fit values, 2a, 3a, and 4a intervals (1 d.o.f.) for the three-flavor neutrino 
oscillation parameters from global data [218]. 



parameter 


best fit 


2(7 


3(7 


4(7 


A77& [10" 5 eV 2 ] 


8.1 


7.5-8.7 


7.2-9.1 


7.0-9.4 


Am 2 31 [10~ 3 eV 2 ] 


2.2 


1.7-2.9 


1.4-3.3 


1.1-3.7 


sin 2 0i2 


0.30 


0.25-0.34 


0.23-0.38 


0.21-0.41 


sin 2 6> 23 


0.50 


0.38-0.64 


0.34-0.68 


0.30-0.72 


sin 2 6*i3 


0.000 


< 0.028 


< 0.047 


< 0.068 



• sin 2 6*23 ~ 0.50, mostly from atmospheric neutrino data 

• sin 2 #13 < 0.30, mostly from atmospheric neutrino and CHOOZ data 

• Am 2 2 ps 8.1 x 10~ 5 eV" 2 , mostly from KamLAND data 

• |Am 2 3 | ps 2.2 x 10~ 3 eV~ 2 , mostly from atmospheric neutrino data 

and allowed regions of the parameters are summarized in Table 5 and Fig. 28. We can 
see that most of the parameters are known with high accuracies. 

Contrastingly, we have rather poor information on some of key neutrino parameters. 
First, only a loose upper bound has been obtained for sin 2 6>i 3 . As we will discuss in 
section ??, this parameter acts an important role in supernova neutrino oscillation. 

The next is the signature of Am 2 3 . There are two mass schemes according to the 
signature (Fig. 29). One is called normal hierarchy with m 2 ps m 2 , <C m\ and another 
is called inverted hierarchy with m\ -C m\ ps m|. The mass scheme is also crucial when 
we consider neutrino oscillation in supernova. 

The value of CP violation parameter (see Eq. (42)) is not known either. Although 
it would have small impact on supernova neutrino oscillation so that we have neglected 
it, it is important in considering the structure and origin of neutrino masses. Also it will 
be important particle-theoretically whether #23 is maximal (7r/4) or not and whether #13 
is exactly zero or just small. 

As we saw in section 3.4.2, the LSND experiment gave us an implication of — > v e 
oscillation with Am1 SND ps leV 2 . However, noting that Am 2 SND ^> Am 2 tm ^> Am 2 olar , 
it is easy to see that three-flavor oscillation scheme discussed above cannot explain the 
LSND data. This is because we have only two independent mass-squared differences with 
three flavors and they are completely determined by the solar, atmospheric, reactor and 
accelerator experiments. If the LSND results are confirmed by another experiment like 
MiniBooNE, we will need some new physics beyond the standard three-flavor neutrino 
oscillation. One possibility is to add an extra neutrino which mixes with standard 
neutrinos. It must not have a charge of weak interaction because the LEP experiments 
imply that there are no very light degrees of freedom which couple to Z-boson [196]. 
Thus the extra neutrino must be sterile. We will not discuss sterile neutrino further 
in this review. For further study about sterile neutrino, see [70, 299] and references 
therein. 
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Figure 28. Projections of the allowed regions from the global oscillation data at 90%, 
95%, 99%, and 3cr C.L. for 2 d.o.f. for various parameter combinations. Also shown is 
Ax 2 as a function of the oscillation parameters sin 2 612, sin 2 623, sin 2 #13, Am| 1; Ato 2 1; 
minimized with respect to all undisplayed parameters [218]. 
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Figure 29. Possible configurations of neutrino mass states as suggested by oscillations: 
the normal (left) and inverted (right) hierarchy. The flavor composition is shown as 
well [89, 90]. 
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4. Neutrino Oscillation in Supernova 

4-1. Overview 

As we saw in section 2, core-collapse supernovae are powerful sources of neutrinos with 
total energies about 10 53 erg. Since neutrinos are considered to dominate the dynamics 
of supernova, they reflect the physical state of deep inside of the supernova, which cannot 
be seen by electromagnetic waves. Neutrinos are emitted by the core and pass through 
the mantle and envelope of the progenitor star. Since the interactions between matter 
and neutrinos are extremely weak, one may expect that neutrinos bring no information 
about the mantle and envelope. In fact, they do bring the information through neutrino 
oscillation because resonant oscillation discussed in section 3.3.2 depends on the density 
profile around the resonance point. Thus neutrinos are also a useful tool to probe the 
outer structure of supernova, including propagation of shock waves. 

On the other hand, supernova has been attracting attention of particle physicist, 
too, because it has some striking features as a neutrino source. As we discussed in 
section 3.4, there have been a lot of neutrino oscillation experiment, which allowed us to 
know many important parameters such as mixing angles and mass-squared differences. 
However, there are still some unknown parameters and physical structure of neutrinos 
which are difficult to probe by the conventional approaches. In this situation, supernova 
has been expected to give us information on fundamental properties of neutrinos which 
cannot be obtained from other sources. 

4-2. Supernova Neutrino 

Here we review the basic properties of neutrinos emitted during various phases from the 
onset of the gravitational collapse to the explosion. Supernova is roughly a blackbody 
source for neutrinos of all flavors with a temperature of several MeV. What is important, 
in the context of neutrino oscillation, is that each flavor has a different temperature, flux 
and its time evolution. The differences are significant especially among v e , v e and the 
other flavors denoted v x . Although the quantitative understanding of the differences are 
not fully established in the current numerical simulation, we can still have qualitative 
predictions and some quantitative predictions. 

4-2.1. neutrino emission during various phases Here let us follow again the supernova 
processes discussed in section 2 focusing on neutrinos. First of all, the core collapse is 
induced by electron capture, 

e~ + A -> v e + A', (119) 

which produces v e s. They can escape freely from the core because the core is optically 
thin for the neutrinos in the early stage of the collapse. However, the luminosity is 
negligible compared with the later phases. 

As the core density increases, the mean free path of u e , X u , becomes smaller due to 
the coherent scattering with nuclei, v e A — > v e . Neutrinosphere is formed when the mean 
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free path A„ becomes smaller than the core size, R. Further, if the diffusion timescale 
of u e , 

3R 2 



is larger than the dynamical timescale of the core, 

1 

+ . ^ ■ 

''dyn r 



(120) 



(121) 



u e s cannot escape from the core during the collapse, that is, u e s are trapped. When 
neutrinos are trapped and become degenerate, the average neutrino energy increases and 
the core become optically-thicker because the cross section of the coherent scattering 
increases as (x coh oc E%. Since low-energy neutrinos can escape easily from the core, most 
of neutrinos emitted during the collapse phase have relatively low energy (< 30MeV). 
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Figure 30. Time evolution 
of neutrino luminosities cal- 
culated by a numerical model 
of a 25M Q progenitor star 
with a 2.05M© core [219]. 
The neutronization peak is 
drawn to be half the actual 
value. Here v, D and fi denote 
v e ,v e and v x , respectively. 




Figure 31. Electron neu- 
trino luminosity and motion 
of mass shells based on the 
same model as the one in Fig. 
30 [219]. The neutrino lumi- 
nosities increase when outer 
core matter accretes onto the 
protoneutron star. 



The collapse of the inner core stops when the central density exceeds the nucleus 
density and a shock wave stands between the inner core and the outer core falling with 
a super-sonic velocity. It should be noted that the shock wave stands in a region with a 
high density (~ 10 12 ~ 10 14 g/cc), which is much deeper than where the neutrinosphere 
is formed (w 10 10 ~ 10 12 g/cc). 
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In the shocked region, nuclei are decomposed into free nucleons. Because the 
cross section of the coherent scattering is proportional to the square of mass number 
(cx co h oc A 2 ), v e can freely propagate in the shocked region. Besides, the cross section 
of electron capture is much larger for free proton than nuclei. Thus a lot of u e s are 
emitted like a burst while the shock wave propagate through the core. This process, 
neutronization burst, works for about 10 msec and the emitted neutrino energy is 
estimated as, 

peak luminosity ~ 10 53 erg s _1 , (122) 

total energy ~ 10 51 erg. (123) 

Some fraction of the shocked outer core accretes onto the protoneutron star, where 
the gravitational energy is converted into thermal energy. Through thermal processes 
like 7 + 7 — > e + + e~, positrons are produced and through processes like, 

e + + n^i? e +p, e + + e~^v + v, (124) 

in addition to the electron capture, neutrinos of all flavors are produced. This accretion 
phase continue for O(10) msec for the prompt explosion and 0(1) sec for the delayed 
explosion. 

Finally, the protoneutron star cools and deleptonizes to form a neutron star. In this 
process, thermal neutrinos of all flavors are emitted with a timescale 0(10) sec, which 
is the timescale of the neutrino diffusion. Dominant production process of the neutrinos 
depends on the temperature: pair annihilation of electrons and positrons e~ + e + — > u+u 
for relatively high temperatures and nucleon bremsstrahlung N + N'—^N + N' + u + u 
for low energies. 

The total energy of neutrinos emitted in the cooling phase of the protoneutron star 
is roughly the same as the binding energy of the neutron star, 



GI^ „ /M NS \ 2 /10km\ 

E m ~-^~3xl0 53 erg ( ^ ) ( J . (125) 



■^« 3 x 10 53 t_ c . , 

i?NS V M Q y 

About 99% of the energy is emitted as neutrinos. 

Fig. 30 shows time evolution of neutrino luminosities calculated by a simulation of 
the Livermore group [219]. This is based on a numerical progenitor model with mass 
25M and a 2.O5M core. The time evolution of v e luminosity is superimposed on the 
motion of mass shells in Fig. 31. The first v e peak is the neutronization burst, whose 
amplitude is drawn to be half the actual value in the figure. The next 0(1) sec is the 
matter accretion phase. The luminosities of v e and u e are greater than that of u x because 
there are additional contributions to u e and u e luminosities from the charged current 
interaction of the pair-annihilation of e + e~. As can be seen in Fig. 31, the neutrino 
luminosities increase when outer core matter accretes onto the protoneutron star. The 
final is the cooling phase of the protoneutron star, during which neutrino luminosities 
decrease exponentially with the neutrino diffusion timescale, O(10) sec. 

To summarize, neutrino emission from a supernova can divided into three phases. 
Because different mechanisms work during the three phases, they have different 
timescales and neutrino luminosities. 
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Figure 32. Time evolution of average neutrino energies [219]. Here v, v and /x denote 
v e , v e and v x , respectively, and the subscripts S and C stand for the calculation method. 
As neutronization of the protoneutron star proceeds, the difference between v e and v x 
energies decreases. 



4-2.2. average energy Average energy of emitted neutrinos reflects the temperature 
of matter around the neutrinosphere. Interactions between neutrinos and matter are 
sufficiently strong inside the neutrinosphere so that thermal equilibrium is realized 
there. Since the temperature is lower in the outer region, neutrino average energy 
becomes lower as the radius of the neutrinosphere is larger. Then the problem is what 
determines the radius of the neutrinosphere. Basically, it is determined by the strength 
of interactions between neutrinos and matter. 

Interactions between neutrinos and matter are, 

u e + n <-> e~ + p, (126) 

v e +p^e + + n, (127) 

u + e ± ^u + e ± , (128) 

v + iV <-> v + N. (129) 

Here it should be noted that the reactions (126) and (127) are relevant only to v e and P e , 
respectively. Furthermore, although all flavors interact with matter through the reaction 

(128) , interactions for u e and V e are contributed from both the neutral and charged 
current, while that for u x is contributed only from the neutral current. Interaction 

(129) occurs equally to all flavors. Therefore, interactions of v e and v e are stronger than 
those of v x . Because there are more neutrons than protons in the protoneutron star, 
u e couples stronger to matter than V e . Thus, it is expected that average energies of 
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neutrinos have the following inequality: 

(E VB ) < (E Pe ) < (E Ux ). (130) 
Although this hierarchy would be a robust prediction of the current supernova theory, 
it is highly difficult to estimate the differences of the average energies without detailed 
numerical simulations. 

Fig. 32 shows time evolution of average neutrino energies obtained from the 
Livermore simulation [219]. We can see the hierarchy of neutrino average energies, 
Eq. (130). The difference between V e and v x energies decreases in time because number 
of protons decreases as neutronization of the protoneutron star proceeds. 

The differences of average energies are important particularly for neutrino 
oscillation. As an extreme case, neutrino oscillation does not affect neutrino spectra at 
all if all flavors have exactly the same energy spectrum. However, prediction of neutrino 
spectra by numerical simulation is highly sensitive and model-dependent although the 
qualitative feature, Eq. (130), is confirmed by a lot of simulations. Simulations by the 
Livermore group [219, 352] predict relatively large differences of average energies, while 
simulations of protoneutron star cooling by Suzuki predict much smaller differences 
[318, 301]. 

4-2.3. energy spectrum The position of the neutrinosphere is determined by the 
strength of the interactions between neutrinos and matter. However, since the cross 
sections of the interactions depend on neutrino energy, the neutrinosphere has a finite 
width even for one flavor. Therefore, the energy spectra of neutrinos are not simple 
blackbodies. Because neutrinos with lower energies interact relatively weakly with 
matter, their neutrinospheres have smaller radii compared to those of high-energy 
neutrinos. As a result, the energy spectrum has a pinched shape compared to the 
Fermi-Dirac distribution. 

Fig. 33 shows energy spectra of V e at different times and the Fermi-Dirac 
distribution with the same average energies [340]. We can see the pinched Fermi-Dirac 
distribution at each time. Time-integrated energy spectra are shown in Fig. ??. 

In the literature, the neutrino spectrum is sometimes parameterized in several 
forms. One popular way to parametrize is called a "pinched" Fermi-Dirac spectrum 
(e.g. [210]), 

L E 2 

F ° a{E) = F( V °)T*eE/T.-*. + V ^ 
where L a and T a are the luminosity and effective temperature of u a , respectively, and 
r) a is a dimensionless pinching parameter. The normalization factor F(r) a ) is 

roc 3 

F M = / ^r— T' ( 132 ) 

J e x + 1 

where F(0) = 77r 4 /120 ~ 5.68. Their typical values obtained from numerical simulations 

are, 

(£? g ) = (14-22)MeV, ^ = (1.1 - 1.6), ^ = (0.5-0.8), (133) 

[Eg) {E s ) 
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Figure 33. Energy spec- 
tra of D e at different times 
and the Fermi-Dirac distribu- 
tion with the same average 
energies [340]. The chemical 
potentials of the Fermi-Dirac 
distributions are set to zero. 



Figure 34. Time-integrated 
energy spectra [340]. 



^ = (0.5-2), ^£ = (0.5-2), (134) 

77e = (0-3), 77 g = (0-3), Tfe = (0 - 2). (135) 

Note that the average energy depends on both T a and r) a , and for i] a = we have 
(E a ) ~ 3.15T Q . 

On the other hand, Keil et al. suggested the following form [170, 171], 

*>,™ i« * ( E V"~' - In E 



^msmrAm ex H-W (136) 

where L a and E a denote the flux normalization and average energy, respectively, and 
P a is a dimensionless parameter that relates to the width of the neutrino spectrum and 
typically takes on values 3.5 — 6. It should be noted that these quantities are dependent 
on both the flavor and time. Spectra obtained from numerical simulations are well fitted 
by 

(E e ) = 12MeV, (E s ) = 15MeV, (E x ) = 24MeV, (137) 

h- = 2.0, y~ = 1.6, (138) 

for the ones by the Livermore group and 

(E e ) = 12MeV, (E^) = 15MeV, (E x ) = (15 - 18)MeV, (139) 

h- = (0.5 - 0.8), ^ = (0.5-0.8), (140) 

for the ones by the Garching group which will be mentioned in the next section. 
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Figure 35. (a) Position of the accretion front as a function of time, (b) Neutrino 
luminosities and rms energies as functions of time. In both figures, thick and thin lines 
show result from AGILE-BOLTZTRAN and VERTEX, respectively. The differences 
in the neutrino results are mainly indirect consequences of the approximate treatment 
of general relativity in the VERTEX simulation. [201] 

4-2.4- recent developments Around neutrinosphere and shock front, neutrinos strongly 
couple the dynamics of different layers on short propagation time scales so that neither 
diffusion nor free streaming is a good approximation in this region. An accurate 
treatment of the neutrino transport and neutrino-matter interactions therefore are 
important not only for obtaining reliable neutrino spectra but also for following the 
dynamics of supernova correctly. It requires to solve energy- and angle-dependent 
Boltzmann transport equation, which is an extremely tough job. However, recent 
growing computer capability has made it possible to solve the Boltzmann equation in 
consistent with hydrodynamics [256, 197, 198, 265, 338, 201]. The result of simulations 
of several groups agree that spherically symmetric models with standard microphysical 
input fail to explode by the neutrino-driven mechanism. 

In [201], Liebendorfer et al. gave a direct and detailed comparison between two 
independent codes of neutrino radiation-hydrodynamics, AGILE-BOLTZTRAN of the 
Oak Ridge-Basel group and VERTEX of the Garching group. Fig. 35 shows (a) the 
position of the accretion front as a function of time, and (b) the neutrino luminosities 
and rms energies as functions of time. The two codes are reasonably consistent as to 
the position of the accretion front. The differences in the neutrino results are mainly 
indirect consequences of the approximate treatment of general relativity in the VERTEX 
simulation. It can be seen that differences of rms neutrino energies are, 

(E Pe ) - (E Ue ) ~ 2 - 3MeV, (E Ux ) - (E Pe ) ~ 3 - 4MeV. (141) 
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However, since recent sophisticated simulations have not succeeded in explosion, it is 
conservative to consider that we do not have definite quantitative predictions about 
neutrino spectra. 

Despite of the advent of the sophisticated simulations, traditional ones by the 
Livermore group are still useful for neutrino oscillation study of supernova neutrinos. 
This is because the latter has a great advantage that it covers the full evolution from the 
collapse over the explosion to the Kelvin-Helmholtz cooling phase of the newly formed 
neutron star, while the former covers typically at most 1 sec. 

4-3. Neutrino Oscillation 

4-3.1. overview Observation of supernova neutrinos can give us information on deep 
inside of supernova which cannot be seen by electromagnetic waves. However, in general, 
neutrinos do not reach the earth as they were produced at the core due to neutrino 
oscillation. As in the case of the solar neutrino, resonant neutrino oscillation occurs 
in the star. In the current case, however, there are two resonance points involving 
three generations of neutrino because neutrinos of all three flavors are produced in the 
supernova and the core density is sufficiently high. 

A key point in the three-generation resonance is the mass hierarchy of neutrino. 
Because there are three generation, there are two mass differences, which are obtained 
by several experiments as 

Am 2 12 «8x 10" 5 eV 2 , 

Am 2 13 ~ Am 2 23 « 2 x 10" 3 eV 2 . (142) 

As was discussed in section 3.4.6, there are two ways to order these two mass difference: 
one is normal hierarchy, for which m\ m m\ <C m|, and another is inverted hierarchy, 
for which m\ <C m\ ~ m|. 

The two resonance points are called H-resonance which occurs at denser region 
and L-resonance which occurs at less dense region. Because the resonance density is 
proportional to the mass difference of the two involved mass eigenstates (see Eq. (82)), 
H-resonance and L-resonance correspond to Am^ and Am 2 2 , respectively. Likewise, 
the involved mixing angles are #13 and #12 for H- and L-resonance, respectively. Thus, 
adiabaticities at the two resonance points depend on the following parameters: 

H — resonance at higher density — > Am] 3 , #13, (143) 
L — resonance at lower density — > Am| 2 , 9 12 . (144) 

The situation is quite different for the two mass hierarchies. If the mass hierarchy is 
normal, there are two resonances in the neutrino sector, while there is one resonance 
in both the neutrino and anti-neutrino sector for the inverted mass hierarchy. The 
schematic views of the two situations are shown in Fig. 36. The vertical line corresponds 
to vacuum and (anti-)neutrino sector is right (left) hand side of it. Here it should be 
noted that v e and v e become effectively heavier and lighter, respectively, in matter if 
neutral-current contributions to the effective mass are subtracted. As can be seen in 
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Fig. 36, L-resonance always occurs in the neutrino sector, while H-resonance occurs 
in the neutrino sector for the normal hierarchy and in the anti-neutrino sector for the 
inverted hierarchy. 




SN core vacuum SN core SN core vacuum SN core 



Figure 36. Resonance schemes for normal mass hierarchy (left) and inverted hierarchy 
(right). 

Neutrinos produced at the core have energy spectra which reflect the physical state 
of the core. Then neutrino oscillation changes the spectra according to the adiabaticities 
of the resonances. Since the adiabaticities depend on the density profile of the progenitor 
star and neutrino oscillation parameters, neutrino spectra observed at the earth are 
determined by: 

• spectra at the core, 

• density profile of the progenitor star, 

• neutrino oscillation parameters (mixing angles, mass differences and mass 
hierarchy) . 

If we observe supernova neutrinos, a mixture of these information will be obtained. 
The former two of the three directly reflect the explosion mechanism of supernova, 
of which we have the basic picture as we discussed in section 2. However, since we 
have not succeeded in reproduce an explosion with a numerical simulation and we have 
not observed supernova neutrinos besides the small number of events (19 events) from 
SN1987A, it would be conservative to think that our understanding of the explosion 
mechanism, especially quantitative understanding, is incomplete. Thus, if we want to 
extract information on neutrino oscillation parameters, we have to conduct an analysis 
based on qualitative features of supernova neutrinos such as, 

• v e emission by the neutronization burst, 

• differences in the average energies and fluxed between flavors. 

However, as we discuss later, number of the neutronization burst events would not be 
statistically sufficient unless supernova occurs several persec from the earth. Anyway, 
supernova is a unique source to probe neutrino oscillation parameters as is implied by 
the comparison to other neutrino sources (Table 6). Supernova neutrinos have quite 
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Table 6. Comparison of some neutrino sources. 



neutrino source 


solar 


atmospheric 


supernova 


production site 


solar center 


atmosphere 


supernova core 


flavor 








energy 


~ 10 MeV 


~ 1 GeV 


10 ~ 70 MeV 


resonance 


once 


none 


twice 


distance 


1AU 


10km 
~ 10 4 km 


lOkpc 
(Galactic center) 



different features from those of other sources so that useful information which cannot 
be obtained from other experiments is expected to be obtained. 

On the other hand, if we want information on supernova itself, we have to subtract 
the effect of neutrino oscillation. Some of the neutrino oscillation parameter, such 
as 2 mass differences, Q\2 and 6*23, have been obtained with excellent accuracies by 
recent experiments as we saw in section 3.4.6. However, only a loose constraint has 
been obtained on 9 13 and we have no idea on the mass hierarchy. Therefore it is not 
straightforward to obtain the intrinsic property of supernova neutrinos. 

4-3.2. resonance points As we saw in Eq. (82), the resonance density is, 



„ /0.5\ /10MeV\ {Am 2 \ 
p res = 1.3 x 10 6 g cm" 3 cos 29 i^—j [^—) [j^2 J > ( 145 ) 

from which, knowing the approximate values of the mass differences Eq. (142), the 
densities of H- and L-resonance regions are obtained as, 

p H ~ (1 - 10) x 10 3 g cm" 3 , (146) 
p L « (20 - 200) g cm" 3 . (147) 

We show a density profile of a progenitor star with mass 15M just before the core 
collapse in Fig. 37, which is based on a numerical model by Woosley and Weaver 
[356, 317]. Electron fraction is also shown in Fig. 38 for reference. As the figure shows, 
the resonance regions are far from the core: 

r H ~ (0.05 - 0.1)^0, (148) 
r L « (0.1 - O.2)i2 . (149) 

Therefore, it is conventional to assume that the dynamics of supernova is not affected 
by neutrino oscillation. Also it has often been assumed that the shock wave does not 
affect the structure of the region where the resonances occur so that static models of 
progenitor star have often been used to analyze neutrino oscillation. In fact, as we 
discuss later, shock wave reaches the resonance region in 0(1) sec and can change the 
adiabaticity of the resonances. 
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Figure 37. Density pro- 
file of a progenitor star with 
15M just before core col- 
lapse [356, 317]. Also shown 
corresponding layers. 



Figure 38. Electron- 
fraction profile of a progeni- 
tor star with 15M Q just be- 
fore core collapse [356, 317]. 



4-3.3. conversion in supernova Here we give a general discussion on flavor conversion 
in supernova following [78]. Our purpose here is to express neutrino fluxes at the surface 
of the star in terms of the original fluxes. Let us denote the original flux of flavor a 
at the core as F®. Because v T , and v T interact with matter almost equally in 
supernova, their fluxes can also be considered to be equal and we define 

F° = F; = F T = F° = F° (150) 

as the non-electron neutrino flux. For possible difference in and v T fluxes, see [3]. 

Although three-generation resonance is complicated in general, the hierarchy of the 
two mass-squared differences and the smallness of U e3 simplify it substantially so that the 
two resonances can be considered to be two two-generation resonances [181, 182, 229]: 

P^l = P„ H ^ x . (151) 

neutrino sector For extremely high densities as in supernova core (p pu,Ph), all 
mixings are suppressed so that the flavor eigenstates coincide with the eigenstates in 
the medium: 

^3,m = Ve, ^2,m = *V, ^l,m = V" ( 152 ) 

Correspondingly, the original fluxes of the eigenstates in medium are written as, 

-^3,m-^ e) ^2,m ~ * r' ~ *x > l,m //,' x ' { Lb6 ) 

We now want to calculate the fluxes of mass eigenstates at the surface of the star 
assuming normal hierarchy. They will be written in terms of the original fluxes 
F®,(a = e,fi,r) and flip probabilities at the H- and L-resonances Ph and P^. First 
let us consider the fate of F e °. At the H-resonance, a fraction P H of F e ° flip to lighter 
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state i>2,m and (1 — Ph) remain to be z/ 3 m . Among the z/ 2m with flux PhF®, a fraction 
Pl flip to the lightest state u ljm and (1 — P L ) remain to be i/ 2>m . As a consequence, 

PlPhP" 

(l-P L )P H P e ° • (154) 
(1 " Pn)F? 

In the same way, Pj are contributed from and iv as, 

( F 1 : (1 - P L )P°, f Pi : P L (1 - P H )P T °, 

P0^<^ P 2 :PlP; , F»^\ P 2 :(1-P L )(1-P H )P T °, -(155) 

[ P 3 : [ P 3 : P H P T , 

Summing all the contributions, we obtain 

F 1 = P L P H P e ° + (1 - P L Ph)P°, (156) 
P 2 = (1 - PL)P H P e ° + (1 - P H + PMFl (157) 
P3 = (l-P H )P e + P H P°, (158) 

which are rewritten as, 

P = a,P e ° + (1 - a;)P°, (159) 

where 

ai = P L P R , a 2 = (l-P L )P H , a 3 = l-P H . (160) 

Because the mass eigenstates are the eigenstates of the Hamiltonian in vacuum, they 
propagate independently to the earth. Further, their coherence is lost on the way to the 
earth so that the neutrinos arrive at the surface of the earth as incoherent fluxes of the 
mass eigenstates. 

Taking into account the neutrino mixing, the net flux of v e is, upto the geometrical 
factor 1/4ttL 2 , 

F e ^ ^ | U e i | F{ 



= pP e ° + (l-p)P°, (161) 

where we used the unitarity condition J2i \ U e i\ 2 = 1 and the "total survival probability 
of u e " , p, is defined as, 

— \TT I 2 

p = \U ei \ di 

= We^P^Pu + WaHl-P^PK + W^il-Pn). (162) 
Since the total flux is conserved, that is, 

P e ° + 2F° = P e + 2F X , (163) 
we obtain the net flux of v x as 

F, + F T = 2F X = (1 - p)P e ° + (1 + p)P°. (164) 
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Thus the fluxes at the surface of the earth can be expressed by the original fluxes and 
the survival probability p. Here it should be noted that not only the fluxes but also the 
survival probability depend on the neutrino energy. 

So far, we have assumed the normal hierarchy. In fact, the case with the inverted 
hierarchy reduces the same expressions with Ph = 1. 



anti-neutrino sector Next we consider the antineutrino sector. Again, the flavor 
eigenstates coincide with the eigenstates in the medium at the core with sufficiently 
high densities: 

(165) 



Vr. 



so that the original fluxes of the eigenstates in the medium are given by, 

Because the small mixing angle 6 e3 is further suppressed in the medium, the V e <-> P 3 
transitions are negligible, Also, the state z/ 3 is so far from the level crossing that it 
propagates adiabatically, that is, v' T — > v%. The other two states can flip through El- 
resonance for the inverted hierarchy. (For a possible L-resonance of anti-neutrino sector, 
see [292].) In the same way as the neutrino sector, we obtain, 

F- e = pF» + {l-p)F x) (167) 
2F x = {l-p)F° e +{l+p)Fl. (168) 



+ Ff 



where 



P=\U el \ 2 P a + \U e3 \ 2 (l-P a ), (169) 
is the effective survival probability of v e . For the normal hierarchy, Ph = 1- 



summary Summarizing the considerations above, the neutrino fluxes can be written 

as, 
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with survival probabilities, 



p = \U el \ 2 P L P H + \U e2 \ 2 (1 - P L )P H + \U e3 f (1 - P H ), 

|2 



(171) 

p=\U el \ z P u + \U e3 \ z (l-P R ), (172) 
which depend on neutrino energy, neutrino oscillation parameters and mass scheme. 
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4-3.4- survival probabilities In the previous section, we obtained an general expression 
of neutrino fluxes at the surface of the earth in terms of original fluxes and survival 
probabilities, which can be expressed by the conversion probabilities at the two 
resonances and oscillation parameters. 

The survival probabilities p and p can be calculated from the following Schrodinger- 
like equation, 
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where the effective Hamiltonian is given by 
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Am( 3 /2E J 
and matter function A(t) is given by 
A(t) = V2G F n e (t). 

When the two resonances are perfectly adiabatic or non-adiabatic, the situation is 
simple. Let us first estimate the adiabaticies of the two resonances. If we assume 
the density profile of the progenitor star as 

p(r) = p r~ n , 

the adiabaticity parameter 7, defined in Eq. (90), is written as, 
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Adopting an approximate value n = — 3 in the mantle, we have 



7 = 9 x W 
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(cos 20 13 ) 4 /3 V2 x 10" 3 eV 2 
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(for H — resonance) 
(178) 



10 d 



sin 2 2#i2 



Amj 2 



2/3 



^lOMeV^ 



2/3 



(cos20 12 ) 4 /3 ^8 x io-5 e v 2 . 
(for L — resonance) 

Substituting the mixing angle obtained from solar neutrino observation into Eq. (179), 
it is seen that the L-resonance is perfectly adiabatic for an energy range of supernova 
neutrino (E = 1 — lOOeV). On the other hand, for the current constraint on 13 , the 
H-resonance can be either perfectly adiabatic or perfectly non-adiabatic. 

Takahashi et al. followed the evolution of survival probabilities, P Ue ^ Ue and Pp e ^p e , 
in [325, 321]. Their calculations were performed by solving the Schrodinger equation 
(173) numerically along the density profile of the progenitor star given in Fig. 37. 
Fig. 39 shows the evolution of the survival probabilities with sufficiently large (L, 
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Figure 39. Survival probabilities of v e (left) and D e (right) in the progenitor star [321]. 
Survival probabilities depend drastically on the mass scheme (normal or inverted) and 
the value of 6*13 (Large or Small). When the resonances are sufficiently adiabatic or 
non-adiabatic, difference in neutrino energy leads to just difference of the position of 
resonance. 



sin 2 2#i3 > lCr 3 ) and small (S, sin 2 29 13 < lCT 5 ) values of 9 13 [321]. Here "normal-L" 
means a model with sufficiently large 6>i 3 and normal mass hierarchy. As can be seen 
the survival probabilities depend drastically on the mass scheme and the value of 9 13 . 
It is also seen that difference in neutrino energy leads to just difference of the position 
of resonance because the resonances are sufficiently adiabatic or non-adiabatic. Here 
it should be noted that mass hierarchy is not important when #13 is sufficiently small 
because, in this case, the H-resonance is perfectly non-adiabatic, which is equivalent to 
the absence of the H-resonance. 




Energy(MeV) time(sec) 



Figure 40. Original neu- 
trino spectra [340] used in 
[321]. 



Figure 41. Original evo- 
lution of neutrino flux [340] 
used in [321]. 



4-3.5. neutrino spectra With survival probabilities obtained in section 4.3.4, 
information of detector, which will be discussed in section 4.5, and original neutrino 
fluxes, we can calculate number of events at the detector. In [321], original neutrino 
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Table 7. Number of events at SK 



hierarchy 


normal 




inverted 






$13 


large 


small 


large 


small 


no osc 




9459 


9427 


12269 


9441 


8036 


v e e~ 


186 


171 


171 


171 


132 


v e e~ 


46 


46 


56 


46 


42 


v x e~ 


98 


98 


77 


98 




0v e 


297 


214 


297 


214 


31 


0P e 


160 


158 


296 


159 


92 


total 


10245 


10114 


13084 


10129 


8441 


burst 


15.7 


16.7 


20.1 


16.7 


12.4 



Table 8. Number of events (CC) at SNO 



hierarchy 


normal 




inverted 






$13 


large 


small 


large 


small 


no osc 


u e d{CC) 


237 


185 


185 


185 


68 


u e d(CC) 


118 


117 


190 


117 


82 


total 


355 


302 


375 


302 


150 


burst 


0.6 


1.1 


1.1 


1.1 


2.1 



fluxes calculated by a realistic model of a collapse-driven supernova by the Lawrence 
Livermore group [352, 340] was used and is shown in Figs. 40 and 41. Figs. 42, 43 and 
44 show time-integrated energy spectra and time evolution of the number of neutrino 
events at SK and SNO. Event number of each interaction is also shown in Tables 7 and 8. 
The distance between the earth and the supernova was set to 10 kpc, which corresponds 
to the galactic center. As in section 4.3.4, "normal" and "inverted" represent the mass 
hierarchy and "L" and "S" mean that #13 is enough large and small for the H- resonance 
to be perfectly adiabatic and non-adiabatic, respectively. 

With a supernova at 10 kpc, SK will have an enormous number of events, which is 
mostly u e p events, and make a statistical study possible. SNO will have much smaller 
events than SK but it can count the number of u e s. Thus it is expected that we can 
obtain many useful information by combining data from the two detectors. However, at 
both detectors, event number of neutronization burst will be rather small. 

In general, neutrino oscillation makes the u e and V e spectra harder, since the 
original average energies of u e and V e are smaller than that of v x . In other words, 
neutrino oscillation produces high energy v e and v e from v x . As a result, the number 
of high-energy events increases and that of low-energy events decreases. The boundary 
between high energy and low energy is around 20 MeV. Note that the amounts of these 
increase and decrease depend on the adiabaticity parameters, and therefore the neutrino 



65 



oscillation parameters, as can be seen in Figs. 42, 43 and 44. This feature can be used 
as a criterion for the magnitude of the neutrino oscillation effects. In [321], the following 
simple ratios were introduced criterion: 

_ number of events at 20MeV < E v < 70MeV 

~~ number of events at 5MeV < E v < 20MeV ' 
The ratios at SK and SNO are plotted in the left of Fig. 45. In this figure, only v e d CC 
events are taken into account for -Rsno and error bars represent the statistical errors 
only. The ratios -Rsk and -Rsno can be considered as estimators of neutrino conversion 
at neutrino sector and anti-neutrino sector, respectively. As can be seen, the ratios will 
give a reasonable implication for the mass hierarchy and the value of 6*13, especially SK 
is expected to give valuable information on 9 13 with its large event number in case of 
inverted hierarchy. However the mass hierarchy cannot be distinguished if #13 is very 
small. Note that the v e flux and the v e flux contain essentially different information 
about the neutrino oscillation parameters. For example, inverted-L and inverted-S are 
distinguishable from V e events, but they are not distinguishable from v e events. 

To this point, we have considered only two extreme cases with perfectly adiabatic 
and non-adiabatic H-resonances. It is also interesting to investigate the intermediate 
cases. In the right of Fig. 45 the 6> 13 dependences of R SK and -Rsno are plotted. Here 
it should be note that -Rsk and -Rsno vary only in the cases of inverted and normal 
hierarchies, respectively, as expected from Figs. 42, 43 and 44. In the case of the 
normal hierarchy, it would be difficult to determine the value of #13, due to the small 
event number at SNO and large statistical errors, but it will give a hint whether it is 
very large or very small. In the case of the inverted hierarchy, the overlap of the error 
bars is small even in the intermediate cases. If #13 is rather large (sin 2 29 i3 > 1CT 3 ), 
these rations will give useful information of the mass hierarchy. 

4-3.6. Earth effects We have discussed neutrino oscillation in supernova and its effects 
on neutrino signal at detectors. Here we will consider matter effect inside the earth. 
Actually, neutrinos can pass through the earth before reaching the detector depending 
on the configuration of the supernova, the earth and the detector (Fig. 46). Fig. 47 
shows path lengths in the earth of neutrinos from the galactic center to various detectors 
as functions of the time of the day. The distance in the earth is normalized by the earth 
diameter. As one can see, at any time of the day, at least one detector among SK 
(Japan), SNO (Canada) and LVD (Italy) observes neutrinos which have experienced 
the earth matter. 

The earth matter can have extra effects on the supernova neutrino spectra and 
observation/non-observation of them can give us further information on neutrino 
oscillation parameters and original neutrino fluxes. First we give a general discussion 
on the earth matter effect following again [78]. 

First, we consider the neutrino sector assuming the normal hierarchy. But it is 
equivalent for the inverted hierarchy if we set Ph = and also for the anti-neutrino 
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sector with the normal and inverted hierarchy if we set Ph — Pl — and Pl — 0, 
respectively. Flux of v e at the surface of the earth, P e , is written as Eq. (161), 

Fe = J2 \ U *\ 2F * = P F e + C 1 - P) F l ( 181 ) 

i 

where Pj is the flux of i-th mass eigenstate and p is the survival probability of v e 
defined in Eq. (162). Now let P« e be the probability that a z/j entering the earth reaches 
the detector as a u e . Then the flux of u e at the detector, P e D , is, 

P e D = ^P e p. (182) 

j 

Rewriting Ps by the original fluxes P° using Eq. (159), we obtain, 

F? = Fe E fl * P - + ^ ^ - E 

= p D P e ° + (l-p D )P°, (183) 
where we have used the unitarity condition J2i Fie = 1 in the first line and, 

P D = E>^e, (184) 

i 

is the effective survival probability of u e at the detector. The earth effect can be 
expressed by the difference between P e and P e D : 

Fe D -F e = (p»-p)(F°-F°). (185) 

The difference of the survival probabilities are computed as, 

P D -P= J2< P ie-\Uie\ 2 ) 

i 

= (1 - 2P L )P H (P 2e - \U 2e \ 2 ) + (1 - P H - P L PH)(P 3e - |t/ 3e | 2 Xl86) 

where ^ Pj e = 1 was used. Because the earth density is rather low and the mixing 
angle 6> 13 is known to be small, we can neglect the second term. Thus we have, 

P e D - P e = (1 - 2P L )P H (P 2e - \U 2e \ 2 ) (P e ° - P°) . (187) 

Now let us evaluate the factor (P 2 e — |^2e| 2 )- For simplicity, we assume a constant density 
and two-flavor oscillation. Remembering that the flavor eigenstates can be written by 
both vacuum mass eigenstates and eigenstates in matter, that is, 





v \*)= v -\Zy (188) 

where U and U m are mixing matrix in vacuum and in matter, respectively, we have 

u- l u n 

COS (0i2,m — #12) Sin (^12,m — ^12; \ I "l,m \ (10,0) 

-sin(0 12 , m -0 12 ) cos(0 12 , m -0 12 ) 11 1 
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Now consider a neutrino which is initially a u 2 and enter the earth. The wave function 
evolves as, 



m 



V2{z) = - sin (0i2, m - On) exp ( ~ i ^^ z ) ^1,: 



+ cos (6>i2, m - On) exp ( -i-^r z ) ^ 
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Then P 2e , the probability that a z/ 2 entering the earth is observed as a u e is, 

^ 9 „ £ sin 2 20 12 9 ( \ / \ 

P 2e = sin 2 12 - NO sin 2 , 191 

(£-cos20 12 ) 2 + sin 2 20 12 VoscmJ 1 ' 

where £ is the dimensionless density parameter defined in Eq. (79), 
2V2G F n B E 



Am 2 



and £ sc,m is the oscillation length in matter defined in Eq. (80), 



£ 

_ *-osc 

*-osc,m 



cos 20) 2 + sin 2 20 

3.1 x m 2 km — 1 ( — - — ] x 

v/(£-cos20) 2 + sin 2 20 VlOMeV; 

8 x l(T 5 eV 



Am 2 

Eventually, we have the final expression: 



(193) 



F? ~F e = (l- 2P L )P U - esin ^ 12 2 x 

1 ' (£-cos20 12 ) 2 + sin 2 20 12 



sin 2 



71Z 



{F° x - F e °) . (194) 



We can easily understand this expression. First of all, there must be a difference between 
the fluxes of v e and v x in order for neutrino oscillation to have any effects. Since flux of v e 
is larger and smaller than that of v x at lower and higher energies, respectively, the earth 
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effect is expected to change sign at some critical energy where F e °(_E crit ) = F®(E crit ). 
This critical energy is about 20 MeV in the simulation of the Livermore group. 

Then there must also be a difference between the fluxes of v\ and v 2 at the surface 
of the earth because we have neglected the contribution from the third generation. This 
fact is reflected in the factor (1 — 2Pl)-Ph~: if -Ph = 0, both v\ and z/ 2 have the original 
spectrum of v x and if Pu = 1/2, both v x and u 2 have the same mixture of the original 
spectra of v e and v x . Therefore it is expected that if #13 is so large (sin 2 2# 13 > 1CT 3 ) 
that the H-resonance is perfectly adiabatic, there is no earth matter effect on u e and u e 
spectra for the normal and inverted hierarchy, respectively. Therefore, existence of the 
earth effect is determined by #13 and the mass hierarchy: 

• small # 13 — > u e and u e 

• large #13 and normal hierarchy — > u e 

• large #13 and inverted hierarchy — > u e 

Thus observation/non-observation of the earth effects will give further information on 
neutrino parameters and mass hierarchy. 

Finally, the remaining factor represents the magnitude of the earth matter effect 
which oscillates with an amplitude |£ sin 2 20 12 /[(£ — cos26> 12 ) 2 + sin 2 26> 12 ]| and with an 
oscillation length £ osc ,m- The amplitude is plotted in Fig. 48 as a function of £. Anti- 
neutrino sector is shown the negative energy region. As one can see, the earth matter 
effect will be larger for higher energies. Comparing the neutrino and anti-neutrino 
sector, it is expected that the earth effect in the neutrino sector is larger than that in 
the anti-neutrino sector because both the above amplitude and the difference of neutrino 
spectrum compared with v x are larger for v e than v e . 

In Fig. 49, fluxes of v e with and without the earth effect are shown. Here it is 
assumed that u e s propagate 4,000 km in matter with the density of the earth core. At 
low and high energies, u e flux increases and decreases, respectively, due to the earth 
effect. As is expected, the earth effect vanishes at the critical energy. Note also that 
there is a modulation in the spectrum, which is a unique feature of the earth effect 
and comes from the factor sin 2 (nz/£ OSCjm ) in Eq. (194). This feature will be further 
discussed in the next section. 

4-3.7. detection and implication from earth effect Here we discuss the earth effects 
more in detail. Because the oscillation length is almost the same order of the earth 
radius, the earth effects depend highly on the path length of neutrinos inside the earth. 
The path length is determined by the position of detector on the earth, the direction of 
supernova and the time of the day. As we saw in Fig. 47, at least one detector among 
SK, SNO and LVD will observe neutrinos with the earth effect and at least one detector 
will observe neutrinos without the earth effect at anytime if a supernova occurs at the 
galactic center. 

If we can see a supernova optically, its direction can be determined accurately. 
However, it is likely that the light from supernova is obscured by interstellar dust if 
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the supernova occurs around the galactic center. Methods of direction determination 
by neutrinos are suggested and studied by several authors. One strategy is to use 
electron scattering events [?]. As we will see in section 4.5, scattered electron in electron 
scattering event reflects the direction of the incident neutrino. In [9], the accuracy of 
the direction determination was estimated to be about 9 degree using expected event 
data at SuperKamiokande. Much better accuracy about 0.6° is expected if a megaton 
water Cherenkov detector is available [336]. Another strategy is triangulation using 
time delay of neutrino signals at different neutrino detectors [47]. However, given the 
expected statistics, this method was shown to be not so effective [32]. 

As we discussed in the previous section, #13 and the mass hierarchy determines the 
existence of the earth effects. Further, since the oscillation length depends on Am^, 
the spectral modulation due to the earth effects will be highly dependent on the value 
of Am\ 2 - In [320], the dependence of the earth effects on Am^ and the path length of 
neutrinos inside the earth was studied in detail. The computational method is essentially 
the same as that without the earth effects: the Schrodinger equation (173) needs to be 
solved along the density profile of the earth. The standard density profile of the earth 
is shown in Fig. 50 [85]. 

Fig. 51 shows event spectra at SK and SNO with earth effect with some values 
of zenith angle and Am^ 2 [320]. In the left row, AmJ 2 is fixed to 2 x 10 _5 eV 2 and, in 
the right row, zenith angle is fixed to zero. It is seen that the spectral shape varies 
with the zenith angle and Am^ 2 . Since the oscillation length is shorter for larger Am^ 2 , 
the frequency of spectral modulation with respect to neutrino energy is larger for larger 
Amf 2 - Thus, two detectors located at separate cites observe neutrinos with different 
energy spectra because their path lengths in the earth are different in general. Also, 
one can see that the earth effect at SK is smaller than that at SNO. This is because the 
dominant events at SK are v e while both v e and u e contribute to the signal at SNO. 

Thus the detection of the earth effect will give us information on #13, the mass 
hierarchy and Amf 2 . There are basically two strategies for detecting the earth effect. 
One is to compare the neutrino fluxes at two or more sites and another is to identify the 
spectral modulation mentioned above. Because the earth effect is a rather small effect, 
the former strategy needs large detectors. On the other hand, the latter strategy needs 
only one detector although it must have a sufficient energy resolution to identify the 
spectral modulation. 

Comparison of neutrino flux at two detector was considered in [76]. This method 
is simple but the detectors must be enough large to accomplish statistically significant 
detection of the earth effect. They considered IceCube and SuperKamiokande, both 
of which are Cherenkov detectors. Because the dominant events come from v e at both 
the detectors, the earth effects are relatively small. Fig. 52 shows the variation of the 
expected IceCube signal with neutrino earth-crossing length. The signal is normalized 
to unity when no earth effect is present (L — 0). The solid and dashed lines correspond 
to the accretion phase and cooling phase, respectively. Here they used neutrino spectra 
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Table 9. Regions in Fig. 53 for the earth effect in IccCubc and Supcr-Kamiokandc. 



Region 


Sky fraction 


Neutrinos come from 
IceCube Super-K 


N SK /N IC 


A 


0.35 


below 


above 


1.070 


B 


0.35 


above 


below 


0.935 


C 


0.15 


below 


below 


Fluctuations around 1 


D 


0.15 


above 


above 


1 



with, 

(El) = 15MeV, 
for the accretion phase and, 

(E° Pe ) = 15MeV, 

for the cooling phase, and neutrino mixing parameters they used Am 2 2 = 6 x 10~ 5 eV 2 
and sin 2 26 12 = 0.9. 

Denoting the number of Cherenkov photon at IceCube as Niq and the equivalent 
IceCube signal measured by SK as Ask, the contour of the ratio Agx/Aic is shown in 
Fig. 53. The sky can be divided into four regions according to the direction of neutrinos 
at the two detectors (Table 9). As one can see, deviation of the ratio Ask /Arc from 
unity is significant when either IceCube or SK observes neutrinos directly and the other 
observes neutrinos from below and the deviation is typically 0.07 in such cases. This 
deviation is rather large because, if a supernova explodes at 10 kpc from the earth, 
the statistical precision for the total neutrino energy deposition is about 0.2% for the 
IceCube and 1% for the SK. Further, such ideal cases is realized for about 70% of the 
sky. Thus the geographical position of IceCube with respect to SK at a latitude of 36.4° 
is well-suited for the detection of the earth effect through a combination of the signals. 

The second possibility, making use of the spectral modulation, was pursued by 
Dighe et al. [77, 75]. Their basic idea is to Fourier-transform the "inverse-energy" 
spectrum of the signal. As we saw in the previous section, the spectral modulation due 
to the earth effect comes from the following factor in Eq. (194), 

Ff-F e ocsm 2 (^|^). (197) 

Therefore the earth effect will have a clear peak in the power spectrum if it is Fourier- 
transformed with respect to E' 1 . The position of the peak is determined by a factor 
Am 2 2m z, that is, Am 2 2 , the density of the region which neutrinos propagate and the 
path length. It is important to note that the peak position is not affected by the primary 
neutrino spectra so that the value of Am 2 2 can be determined accurately from the peak 
position independently of supernova model if we know the neutrino path length inside 
the earth. 



«) = 17MeV, |f = l-5, (195) 



(E°J = 18MeV, S = 0.8, (196) 



71 



To illustrate the effectiveness of their method, they consider a 32 kton scintillator 
detector and a megaton water Cherenkov detector. The major difference between the 
scintillation and the water Cherenkov detector is that the energy resolution of the 
scintillation detector is roughly six times better than that of the Cherenkov detector. 
Therefore even the volume of the SuperKamiokande will not be sufficient to identify the 
earth effect and much larger detector, HyperKamiokande (HK), is needed. 

Figs. 54 and 55 show expected power spectra at the scintillator detector and HK, 
respectively, for different SN models, Garching (G) and Livermore (L), and distances 
traveled through the earth. Supernova is assumed to be 10 kpc from the earth and the 
power spectra are averaged over 1000 SN simulations. The left rows show the power 
spectra of neutrinos which propagate only in the earth mantle while the right rows 
correspond to those which propagate the earth core as well as the mantle. 

A clear peak is seen in each figure in the left rows, corresponding to the assumed 
values of Am^ 2 , the mantle density, the average neutrino energy and the path length. 
The peak width reflects the finite energy resolution of the detectors and weak energy 
dependence of Eq. (194) besides the factor Eq. (197). The large power at small k is the 
dominant contribution from the energy dependence of the neutrino spectra without the 
earth effect. On the other hand, if neutrinos propagate in both the mantle and core, 
multiple peaks appear as in the right rows of the figures corresponding to their densities 
and path lengths in them. 

The Livermore model was used as a representative model of conventional 
simulations which predict relatively large average-energy differences between flavors. 
On the other hand, the Garching model, which treats neutrinos more carefully, predicts 
rather small average-energy differences between flavors. As we saw in the previous 
section, the earth effects as well as the effects of neutrino oscillation in supernova are 
larger for a case with larger average-energy differences. However, as seen in the figure, 
the peak positions are independent of the supernova models. 
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Figure 42. Time-integrated energy spectra (left) and the time evolution of the number 
of neutrino events (right) at SK. Only v e p CC interaction, which is the dominant event 
at SK, is taken into account [321]. 




Figure 43. Time-integrated energy spectra (left) and the time evolution of the number 
oiv e d events (right) at SNO. Neutronization burst is suppressed for normal-S, inverted- 
L and inverted-S, and absent for normal-L [321]. 




Energy (MeV) time (sec) 



Figure 44. Time- integrated energy spectra (left) and the time evolution of the number 
of D e d events (right) at SNO [321]. 
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Figure 45. (left) Ratios of high-energy to low-energy events at SK and SNO. (right) 
Dependence of the ratios on sin 2 6*13. In both figures, only v e d events are taken into 
account for R(SNO) and error bars represent the statistical errors only [321]. 
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Figure 46. Neutrinos pass- 
ing through the earth before 
reaching a detector. 
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Figure 47. Path lengths in the earth 
of neutrinos from the galactic center to 
various detectors as functions of the time 
of the day. The distance in the earth is 
normalized by the earth diameter. 
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Figure 48. Amplitude 
of the earth matter effect 
£sm 2 20 12 /[(£ - cos26>i 2 ) 2 + 
sin 2 20i2] as a function of £. 
Here we set sin 2 0i2 = 0.3. 
Anti-neutrino sector is shown 
in the negative energy region. 
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Figure 49. Fluxes of v e with 
the earth effect (dashed line) 
and without the earth effect 
(solid). Here it is assumed 
that v e s propagate 4,000 km 
in matter with the density of 
the earth core. 
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Figure 50. Density profile of the earth [85]. 
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Figure 51. Event spectra at SK and SNO with earth effect with some values of zenith 
angle and Amf 2 [320]. In the left row, Amf 2 is fixed to 2 x 10~ 5 cV 2 and, in the right 
row, zenith angle is fixed to zero. Only charged-current events are considered in the 
SNO spectra. The mixing angle 9i 3 is taken sufficiently small that the H-resonance is 
perfectly non-adiabatic. 
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Figure 52. Variation of 
the expected IceCube signal 
with neutrino earth-crossing 
length [76]. The signal is 
normalized to unity when no 
earth effect is present (L = 
0). The solid and dashed 
lines correspond to the accre- 
tion phase and cooling phase, 
respectively. 
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Figure 53. Contours of 
Nsk/Njc °n the map of the 
sky projected on the earth 
[76] 
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Figure 54. Averaged power spectra in the case of a large scintillator detector for 
AM^ r ,„i- cm ri„„„u; — , (n\ anc j Livermore (L), and distances traveled through 
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Figure 55. Same as the Fig. 54 but for the case of HK [75]. 
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Figure 56. Left: Realistic spectrum from a single simulation. Right: Area distribution 
of the background (black) and the signal (red) obtained for a 32 kton scintillator 
detector and Garching model for 77 = 60 [75]. 




Figure 57. Left: Comparison of P95 and pgg for the Garching (G) and Livermore (L) 
SN models in a 32 kton scintillator detector. Right: Comparison of P95 in this large 
scintillator detector (SC) and in the case of a megaton water Cherenkov (HK), for the 
Garching model. From [75] . 

To quantify the effectiveness of their method, they introduced an algorithm to 
identify the peaks in the real world with the presence of background fluctuations. It is 
based on the integration of the area around the expected position of the peak. Once we 
know the neutrino path, we can calculate the position of the peaks. It is more robust 
to consider the area around the expected position of the peak than to look for the 
maximum in the height of the power spectrum, as shown in the left of Fig. 56. They 
considered, for a single peak case, the interval of integration /c pea k ± AA; with A A; = 30, 
which is roughly the expected width of the peak. For a multiple peak case, they measure 
the area k = 40 - 160. 

Then we must consider the statistical significance of the result obtained. For this 
purpose, they compared the value of the measured area with the distribution of the area 
in the case of no earth effect. The right of Fig. 56 shows the area distribution of Monte 
Carlo simulations with and without the earth effect. Here A 95 ~ 100 denote the area 
corresponding to 95% C.L. detection of the earth effect. The problem is the probability 
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that we have larger area than A 95 , which is denoted by P95. This probability depends on 
the distance traveled by the neutrinos through the earth, which is in turn determined 
by the location of the supernova in the sky. In Fig. 57, the probability as a function 
of the nadir angle r\ of supernova is plotted. The passage through the core corresponds 
to rj < 33°. The left is comparison of P95 and pg$ for the Garching and Livermore SN 
models in a 32 kton scintillator detector. The right is comparison of P95 in this large 
scintillator detector (SC) and in the case of a megaton water Cherenkov (HK), for the 
Garching model. As can be seen, this method is very effective if the average-energy 
differences of the original neutrino spectra are as large as predicted by the Livermore 
simulation, while the effectiveness decreases by half with the Garching simulation. 

4-3.8. shock wave and neutrino oscillation So far we have assumed the structure of 
supernova envelope to be static during ~ lOsec of neutrino emission. Actually, as 
pointed out by Schirato and Fuller [282] and studied further in [323, 211, 110, 337, 93], 
the shock wave produced at the bounce reaches the resonance regions in several seconds 
(left of Fig. 58) and can affect the adiabaticities of the resonances. This fact implies 
the possibility to probe the propagation of the shock wave by neutrino observation. 




Figure 58. Left: Propagation of forward and reverse shocks with densities of H- and 
L-resonance regions [337]. The width of the density bands reflects the energy range 
of supernova neutrinos. Right: Schematic density profile in the presence of a forward 
and reverse shock wave with a contact discontinuity between them. Four key densities 
are denoted as pib, P2t>, P2a and p\ a from dense to sparse region, which correspond to 
the edges of the forward and reverse shocks. 



Before we discuss neutrino oscillation, let us summarize basic features of shock 
wave. While the shock sweeps the infalling matter, the gain region behind the shock is 
accelerated by the neutrinos from the neutron star to form the neutrino-driven wind. 
At the interface between this neutrino-driven wind and the shock-accerelated ejecta is 
formed a contact discontinuity, which is characterized by a density jump but continuous 
velocity and pressure. Even farther behind the forward shock, a reverse shock can appear 
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due to the collision between the neutrino-driven wind and more slowly moving material. 
The forward and reverse shocks are sharp discontinuities where density, pressure and 
velocity change on the microscopic (sub-millimeter) scale of the ion mean free path. 
The right of Fig. 58 shows a schematic density profile in the presence of a forward 
and reverse shock wave with a contact discontinuity between them. The existence of 
the contact discontinuity and reverse shock seems to be a generic feature in supernova 
explosion according to recent numerical simulation. However, because the structure 
and time evolution of the shock wave depend on the detailed dynamics during the 
early stages of the supernova explosion, it is difficult to make an exact prediction. 
The situation becomes even more complicated if we consider non-spherical supernova 
explosion where violent convective instabilities and large anisotropies make the density 
structure chaotic, although the some generic features of the one-dimensional situation 
are retained. Here we discuss effects of shock propagation on neutrino oscillation in 
one-dimensional supernova following [337]. 



P 




Figure 59. Left: Survival probability as a function of energy with cos 2 612 = 0.7. The 
dotted line corresponds to the probability based on the simple argument in the text 
while the solid line is based on a more realistic analytic evaluation in [337]. Right: 
Survival probability p(E, t) as a function of energy at different times averaging in 
energies with the energy resolution of Super-Kamiokande: for a profile with only a 
forward shock (left) and a profile with forward and reverse shock (right). At t = 5 sec, 
p(E, t) including Earth matter effects for a zenith angle of 62° is also shown with a 
black line. Both are from [337]. 

First, let us consider the survival probability of neutrinos which experience the H- 
resonance, that is, z/ e s and v e for the normal and inverted hierarchy, respectively. In this 
section, we assume the inverted hierarchy and consider v e . We can obtain a qualitative 
feature of the survival probability as a function of neutrino energy from the schematic 
density profile of the shock shown in the right of Fig. 58. In the figure, four key densities 
are denoted as pit,, P2b, P2a and pi a from dense to sparse region. They correspond to the 
edges of the forward and reverse shock: the density jumps from p 2a to p 2 b at the reverse 
shock and from pib to pi a at the forward shock. As we saw in Eq. (82), the resonance 
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Figure 60. Upper and lower left: The average energy of Dp —* ne + events at a 
megaton water Cherenkov detector binned in time for simulations by the Livermore 
group (L) and the Garching group (Gl and G2). The average energy is assumed to be 
static and the error bars represent la errors in any bin. Lower right: Time dependence 
of (E e ) for a profile with a forward and reverse shock for several values of tan 2 813 and 
for model L. Figures from [337]. 



density is proportional to the neutrino energy, 

1 1 „ /0.5\ /10MeV\ / Am 2 \ , , 

so that there are four energies whose resonance densities are the above four densities, 
E\hi Eib, E2& and E\ a from lower to higher energies. Here it should be noted that the 
four densities are time-dependent and so are the four energies. 

The essential point in the effect of shock propagation on neutrino oscillation is that 
the density gradient at the shocks is so steep that the resonance there tends to be non- 
adiabatic even if the mixing angle involved is rather large. Let us assume first that the 
resonance is completely non-adiabatic at the shocks and completely adiabatic elsewhere. 
The latter means that #13 is sufficiently large. Then, a v e with energy En, < E < E 2 b 
experiences the H-resonance three times, that is, behind the reverse shock, between the 
reverse and forward shock, and at the forward shock, which are adiabatic, adiabatic and 
non-adiabatic, respectively. The three resonances are equivalent to one non-adiabatic 
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resonance. In the same way, we obtain the adiabaticities of resonances of various 
neutrino energies, 

• E < E lh =>- adiabatic 

• Ei b < E < E 2h =>■ adiabatic, adiabatic, non-adiabatic effectively non-adiabatic 

• E 2 b < E < E 2a =>■ adiabatic, non-adiabatic, non-adiabatic =>- effectively adiabatic 

• E 2a < E < E la =>- non-adiabatic 

• E > E la =>- adiabatic 

and the survival probability as a function of neutrino energy become like the left of Fig. 
59. As one can see, there are two peaks in the survival probability, which result, as we 
will see later, in characteristic features in the time evolution of observed neutrino average 
energy and event number. Note that there will be only one peak without the reverse 
shock and that the behavior of the survival probability will become more complicated 
if the contact discontinuity is so sharp that resonance there is non-adiabatic. Anyway, 
since we assumed #13 to be sufficiently large so that the survival probability is zero in 
absence of shock propagation, the effect of the shock can be thought to be to make the 
resonance more non-adiabatic effectively. 

Because the shock propagate into low-density regions, the key densities decrease in 
time and so the key energies increase in time. The right of Fig. 59 shows the survival 
probability p(E, t) as a function of energy at different times averaging in energies with 
the energy resolution of Super-Kamiokande: for a profile with only a forward shock (left) 
and a profile with forward and reverse shock (right). At t — 5 sec, p(E,t) including 
earth matter effects for a zenith angle of 62° is also shown with a black line. It can be 
seen than there is only one peak without the reverse shock and are two peaks with both 
the reverse and forward shocks and the peaks shift to high-energy regions as the shock 
propagates, as expected. 

Upper and lower left of Fig. 60 show the average energy of u e p — > ne + events 
binned in time for simulations by the Livermore group (L) and the Garching group (Gl 
and G2). Magenta, red and blue lines correspond to a static density profile, a profile 
with only a forward shock and with forward and reverse shock, respectively. The error 
bars are la statistical errors assuming a megaton water Cherenkov detector. As one can 
see, there are some dips in the time evolution of the average energy: one dip with only 
forward shock and double dips with both forward and reverse shock. This characteristic 
feature is the direct result of the dips in the survival probability in Fig. 59 and can be 
seen in every SN model. Lower right of Fig. 60 is the time dependence of the average 
energy for a profile with a forward and reverse shock for several values of tan 2 9 13 and 
for model L. It is seen that the dips are smaller for smaller values of tan 2 # 13 . This 
can be understood by remembering that the propagation of the shock waves through 
the resonance region tends to make the resonance more non-adiabatic effectively. When 
tan 2 9 1 3 is small, the H-resonance is non-adiabatic even without the shock and then the 
shock has no effect on neutrino oscillation. 
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Figure 61. Left: Total number of events as a function of time for a static density 
profile and a profile with only forward shock and with both forward and reverse shock. 
Right: Same as the left but the energy range is restricted to E — 40 ± 5MeV. The 
error bars are statistical errors expected in a megaton water Cherenkov detector. Both 
are from [337]. 



The imprint of the shock waves can also seen in the time evolution of the event 
number. The left of Fig. 61 shows the total number of events as a function of time for a 
static density profile and a profile with only forward shock and with both forward and 
reverse shock. As one can see, a dip can be seen in the presence of both the forward and 
reverse shocks. This is even striking if the energy range is restricted to E = 40 ± 5MeV 
as can be seen in the right of Fig. 61. 

Thus, if 6*13 is sufficiently large (tan 2 6*13 > 10~ 5 ), dips appear in the time evolution 
the neutrino average energy and event number, which would be relatively model- 
independent. A megaton water Cherenkov detector, like the HyperKamiokande, can 
probe the shock wave propagation efficiently if the mass hierarchy is inverted one so 
that the H-resonance occurs at the anti-neutrino sector. 
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Figure 62. Density profiles (left) and scale height n e /\dn e /dr\ of stars just before 
supernova explosion with initial mass 15M©, 20M©, 4OM and 75M [322]. 
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4-3.9. toward model-independent predictions While type la supernovae have rather 
universal features as far as we observe them, core-collapse supernovae seem to have a 
wide variety in luminosity and spectrum. Since this will reflect the diversity of the 
presupernova structure, neutrino emission will also be dependent on the physical state 
of the progenitor star. Some important physical quantities about the progenitor star 
will include its mass, metalicity, magnetic field and rotation, whose effects on neutrino 
emission are still unclear. If we want to extract information about neutrino oscillation 
parameters, the effects of these supernova parameters on the neutrino emission and 
dynamics of neutrino oscillation must be studied and some model-independent analyses 
are required. In [322], Takahashi et al. studied the effect of initial mass of a progenitor 
star on the neutrino emission and neutrino oscillation. Then some input-physics 
dependences of neutrino emission were studied in [162] concentrating on neutronization- 
burst signal. Here we review [322] and summarize [162] briefly. 

The mass of the progenitor star affects the neutrino oscillation signature of 
supernova neutrinos through differences in both the mantle and core structures. The 
density profile of the progenitor star, especially of the mantle, is important because it is 
related to the dynamics of neutrino flavor conversion. On the other hand, the structure 
of the iron core at the collapse determines the characteristics of the neutrino burst, e.g., 
the average energy and luminosity for each flavor. 

Let us first discuss the density profile of a progenitor star. In numerical simulations 
of supernova, the initial condition is a star just before the collapse whose structure 
is often given by a numerical presupernova model, which is obtained by following the 
evolution of a massive star. The evolution of a massive star is significantly affected by 
mass loss due to a stellar wind, which we still do not have a definite understanding. 
Indeed, the mass loss can become so strong for a star with initial mass more than about 
35 M Q and solar metalicity that the entire hydrogen envelope can be lost prior to the 
explosion of the star. It is suggested in [355] that the maximum in the final mass is 
about 20 Mq. Thus, a massive star just before the supernova explosion may have a 
universal structure independent of the initial mass of the star. The final density profiles 
of stars with various initial masses just before the collapse are shown in the left of Fig. 
62 [317]. As is expected, they are similar to each other. 

The density profile of the star comes into the adiabaticity parameter as the scale 
height n e /\dn e /dr\. A smaller scale height, that is, a steeper density profile results in less 
adiabatic resonance. The scale heights of stars with various initial masses, calculated 
from the density profiles shown in the left of Fig. 62, are given in the right of Fig. 62 
as functions of the density. The scale heights vary significantly, but, at the densities 
relevant for the resonances, differences between different initial masses are factors of 2 
or 3. 

In the left of Fig. 63, the evolution of the survival probability P[y e — > v e ) is shown. 
The neutrino energies on the plot are 5 MeV and 40 MeV and we set sin 2 2# 13 = 10 -4 . 
The H-resonance radii and the final probabilities can be quite different for different 
progenitors and neutrino energies. The observationally important quantity is the final 
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Figure 63. Left: Evolution of survival probabilities P(y e — > v e ) as functions of 
radius for different progenitor models. Upper and lower curves correspond to neutrino 
energies of 40 MeV and 5 MeV, respectively. Right: Energy dependence of survival 
probability P{y e — > v e ) for different progenitor models. Also shown are the perfectly 
adiabatic and non-adiabatic cases for the H-resonance. 



probability and this is shown in the right of Fig. 63 as a function of the neutrino energy. 
The differences are not so small, about O(10)% at all energies. But if sin 2 2# 13 is so 
large or very small that H-resonance is perfectly adiabatic or non-adiabatic, respectively, 
difference in scale height will not affect the neutrino conversion probabilities. 

Next we consider the difference in core structure. As was stated in the above, 
massive stars experience significant mass loss. For current empirical mass loss rates, all 
solar-metalicity stars initially more massive than about 35 M are thought to become 
hydrogen-free objects of roughly 5M at the end of their thermonuclear evolution. The 
corresponding upper limit to the mass of the final iron core is about 2M [355]. 

As mentioned in subsection 2.2, the mass of the iron core is determined roughly 
by the Chandrasekhar mass. For a zero-temperature and constant Y e , its Newtonian 
structure is given by, 

M C ho = 5.83F e 2 M . (199) 

However, there are numerous corrections, some of which are large [339]. To the first 
approximation, the non-zero entropy of the core is important and 

2" 

(200) 



M Ch ~ M Ch0 



1 + 



where 



s ==°- 5 °(io4^) 1/3 (o^r (ji^v) (2oi » 

is the electronic entropy per baryon. More massive stars have higher entropy and contain 
larger iron cores on average. However, this general tendency is moderated by the loss 
and redistribution of entropy that occurs during the late burning stages. Thus, the 
mass of the iron core as a function of the initial mass will be somewhat uncertain in 
that a small change in the initial mass results in a large difference in the iron core mass. 
According to [355], the mass of the iron core is 1.2(1.4) — 1.6M when the initial mass 
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is between 1O(2O)M and 4OM . This weak dependence of the iron core mass on the 
ZAMS progenitor mass leads to a somewhat universal neutrino burst. 

Fig. 64 shows the evolution of the average neutrino energy and number luminosity 
in the early phase up to 200 milliseconds after bounce. The calculation is based on 
dynamical models of core-collapse supernovae in one spatial dimension, employing a 
Boltzmann neutrino radiation transport algorithm, coupled to Newtonian Lagrangean 
hydrodynamics and a consistent high-density nuclear equation of state. Details of these 
simulations are described in [338]. As can be seen, the major features of the early 
neutrino burst are almost independent of the initial mass. 

Combined with the discussion on the scale height above, we conclude that the 
mantle structure and the features of the neutrino burst depend little on the initial mass 
of the progenitor star if sin 2 2^13 < 10~ 5 or sin 2 2^13 > 10~ 3 . On one hand, this means 
that we can not easily obtain information about the initial mass from observations 
of neutrinos during the first 200 milliseconds after bounce. On the other hand, this 
situation is desirable for extracting information about the neutrino parameters. 

In [162], Kachelriess studied dependence of neutronization burst on input neutrino 
physics in numerical simulation, as well as progenitor mass. They found neutronization 
burst to be relatively independent of the progenitor mass, electron capture rate 
and equation of state at high densities. Especially, uncertainties in the number of 
neutronization-burst events due to the input physics was estimated to be less than 10%. 
This feature can be used not only to probe neutrino parameters such as #13 and the mass 
hierarchy but also to determine the distance to the supernova. The latter is important 
because it is likely that the supernova is optically obscured by dust if it occurs around 
the galactic center. Their estimation is that distance to the next galactic supernova is 
determined with a 6% error if we have a megaton water Cherenkov detector. 

As we have discussed above, effects of some of neutrino parameters and input 
physics on supernova neutrinos and their oscillation have been studied so far. 
Fortunately to particle physics, the progenitor mass, electron capture rate and equation 
of state have relatively small effects on neutrinos and we will not suffer from uncertainties 
due to them. However, analyses so far are rather restricted to one-dimensional 
supernova, that is, we do not understand multi-dimensional effects such as magnetic field 
and rotation. This is because multi-dimensional simulation with sufficiently detailed 
treatment of the neutrino physics is still challenging now and much progress in this field 
is indispensable. 
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Figure 64. Left: Evolution of the number flux at the Earth of neutrinos from a 
supernova by a progenitor with initial masses 15M Q , 20M©, 40M Q and 75M© at a 
distance of 10 kiloparsecs. Right: Evolution of neutrino average energies. In all figures, 
time at the bounce is set to zero. 
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Table 10. Arrival times of neutrinos, energies of prompt electrons and angles between 
momenta of neutrinos and corresponding prompt electrons at Kamiokande II and 1MB. 



KII 


1MB 


event 


time 


energy 


angle 


event 


time 


energy 


angle 




(sec) 


(MeV) 


(cleg) 




(sec) 


(MeV) 


(cleg) 


1 


0.000 


20.0 ± 2.9 


18 ± 18 


1 


0.000 


38 ± 7 


80 ± 10 


2 


0.107 


13.5 ± 3.2 


40 ± 27 


2 


0.412 


37 ± 7 


44 ± 15 


3 


0.303 


7.5 ± 2.0 


108 ± 32 


3 


0.650 


28 ± 6 


56 ± 20 


4 


0.324 


9.2 ± 2.7 


70 ± 30 


4 


1.141 


39 ± 7 


65 ± 20 


5 


0.507 


12.8 ± 2.9 


135 ± 23 


5 


1.562 


36 ± 9 


33 ± 15 


6 


1.541 


35.4 ± 8.0 


32 ± 16 


6 


2.684 


36 ± 6 


52 ± 10 


7 


1.728 


21.0 ± 4.2 


30 ± 18 


7 


5.010 


19 ± 5 


42 ± 20 


8 


1.915 


19.8 ± 3.2 


38 ± 22 


8 


5.582 


22 ± 5 


104 ± 20 


9 


9.219 


8.6 ± 2.7 


122 ± 30 










10 


10.433 


13.0 ± 2.6 


49 ± 26 










11 


12.439 


8.9 ± 1.9 


91 ± 39 











44. Neutrinos from SN1987A 

On February 23 in 1987, a supernova was found in the Large Magellanic Cloud. This 
is the closest supernova from us since a galactic supernova in 17th century. From 
various analyses, this supernova, SN1987A, was identified as a type II supernova whose 
progenitor star is a blue giant with mass M ~ 1OM at 50 kpc. Here we review the 
basic facts of neutrinos from SN1987A and their implication. 
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Figure 65. Arrival times and es- 
timated energies of neutrinos ob- 
served at Kamiokande II and 1MB 
[134]. 
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Figure 66. Detection ef- 
ficiencies at Kamiokande 
II [134] and 1MB [41]. 



4-4-1- observational facts Just after the discovery of SN1987A by optical observations, 
it was expected that neutrinos from such a close supernova must have been detected. 
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Actually Kamiokande II [133, 134] and Irvine-Michigan-Brookhaven detector (1MB) 
[36, 41] observed 11 and 8 events, respectively. (For possible detection of neutrinos at 
Baksan, see [5].) 

Table 10 shows the arrival times of neutrinos, energies of prompt electrons 
and angles between momenta of neutrinos and corresponding prompt electrons at 
Kamiokande II and 1MB. Arrival times and estimated energies of neutrinos observed 
at the Kamiokande II and 1MB are plotted in Fig. 65. From this figure, one might 
think that neutrinos detected at the 1MB have higher energies than those detected at 
the Kamiokande II. But it is not obvious because we have to consider the difference in 
detection efficiencies showed in Fig. 66. As one can see, the Kamiokande II had a high 
efficiency at more than 20 MeV while the 1MB was not effective to detect low-energy 
neutrinos. We summarize basic information on the observed neutrinos below. 

duration The duration of neutrino events is 12.4 sec at the Kamiokande II and 5.6 sec 
at the 1MB. This is consistent with the diffusion timescale of neutrino, about 10 sec, 
discussed in 4.2. Thus it is confirmed that neutrinos are confined in the protoneutron 
star and escape by diffusion. 

angular distribution of events Since both the Kamiokande II and 1MB are water 
Cherenkov detectors, main events come from u e p — > e + n which has about hundred 
times larger cross section than that of electron scattering ve~ — > ve~ . A positron is 
emitted isotropically in u e p — > e + n while an electron scattered by a neutrino has a 
forward peak. Angular distribution in Table 10 confirms these arguments. 

Let us focus on the first event at the Kamiokande II. In this event the primary 
electron is emitted forward. The probability that an electron, which has an isotropic 
distribution, is emitted forward inside 20° is 3%, which leads to the expectation value 
of 0.6 for total event number of 19. On the other hand, as we saw in Table 7, event 
number of electron scattering is about 5% of that of u e p — > e + n event. Thus, it is not 
obvious which reaction the first event came from. If the electron scattering is the case, 
u e is the most likely for the event because v e has the largest cross section in electron 
scattering. 

neutronization burst? If the first event at the Kamiokande II was v e event, it is possible 
that the v e is emitted during the neutronization burst. However, the event number from 
the neutronization burst was estimated to be about 0.01 in [280]. Thus it is unlikely 
that the first event at the Kamiokande II was from the neutronization burst. 

neutrino temperature and luminosity Neutrinos from supernova have roughly a thermal 
distribution though they do not exactly. The effective temperature and luminosity of 
the observed P e were estimated by several authors [11, 152]. They are roughly consistent 
with each other and give, 

T Pe = (3 - 4)MeV, L Pe = (3 - 6) x 10 52 erg. (202) 
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If we assume that all flavors have the same luminosity, the total neutrino luminosity 
is about 3 x 10 53 erg. This is about the same as the binding energy of a neutron star, 
which indicates that the current supernova theory is roughly correct. 

4 -4 constraints on neutrino properties Observational feature of the neutrinos from 
SN1987A can be summarized as follows: 

• Duration of the neutrino events is about 12 sec 

• Temperature of v e is about (3 — 4) MeV 

• Neutrino total energy is roughly the same as the binding energy of a neutron star 

We can put constraints on any new physics, unknown processes and exotic particles 
which lead to any contradictions with the above observational facts. Here we will discuss 
topics related to neutrinos. 

mass If neutrinos have mass, neutrinos with different energies have different velocities, 
vu « 1 - (203) 

so that the propagation time from the SN1987A to the earth are also different. Denoting 
the departure and arrival times as t d and t a , respectively, the propagation time can be 
written as 

*--*' = ir D ( 1 + 5S)' (204) 

where D ~ 50kpc is the distance between SN1987A and the earth. Therefore the 
difference in the propagation times of two neutrinos is, 

|Af a - Af d | = \ D ^%2El^ ( 205 ) 

where At a and At d are the differences in the departure and arrival times, respectively 
of the two neutrinos. As is stated above, At a < 12sec for any two of the 19 events. To 
obtain a constraint on neutrino mass, we need a statistical analysis. Let us first consider 
a simple case with the event 3 and 10 at Kamiokande II. In this case At a = lO.lsec, and 
then we have 

|10.1sec- Ai e | = 0.06sec . (206) 

If we assume At d <C lOsec, we obtain m v < 13eV. Statistical studies give slightly weak 
constraints, m v < (19 - 30)eV [14, 23, 173, 280] 

lifetime Since the estimated total energies of the emitted neutrinos is roughly the same 
as the binding energy of a neutron star, it can be said that most of the neutrinos did 
not decay before reaching the earth. Thus, denoting neutrino lifetime as t v , 

—r u > — ^ 5 x 10 12 sec, (207) 
m u c 
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then we have 

Note that neutrinos cannot decay without having mass. 
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Figure 67. Deflection of path of a charged particle by magnetic field. 



electric charge If neutrinos have electric charge, even if it is extremely small, neutrino 
trajectory is deflected by the galactic magnetic field. Then the trajectory of a low-energy 
neutrino is longer than that of a high-energy neutrino. Larmor radius of a neutrino with 
charge Q v in a magnetic field _B ga i is 

R (209) 
Qi/eBga\ 

From Fig. 67, the extra distance due to charge is 
Then time delay of two neutrinos is, 



^=^ + ^w^(M-m)- (211) 

Taking At d < 20sec and B gal = 10~ 6 Gauss, we obtain Q u < 10~ 18 . A statistical analysis 
considering the inhomogeneity of the galactic magnetic field gives Q u < 10~ 17 [25]. 

weak equivalence principle The difference between arrival times of photons and 
neutrinos was several hours. This means that the gravitational constants for them 
are not so different. Specifically, 
G,, — G^ 



G v + Gj 
is obtained in [205, 180]. 



< 10" d , (212) 
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4-4-3. neutrino oscillation Here we discuss neutrino oscillation of neutrinos from 
SN1987A. Because most of the observed neutrinos were u e , we concentrate on the 
anti-neutrino sector. Then the key point of the neutrino oscillation dynamics is the 
adiabaticity of the H-resonance if it exists. If the mass hierarchy is inverted and #13 is 
so large that the H-resonance is perfectly adiabatic, u e flux can be written as (see Eqs. 
(170) and (172)), 

F- e = \U e3 \ 2 F° + (l-\U e3 \ 2 )F° 

« Fl (213) 

where we used \U e3 \ 2 <C 1. Thus the u e flux observed at the earth reflects directly the 
original flux of v x . 

On the other hand, if the mass hierarchy is normal or #13 is very small, the situation 
becomes more complicated, which was studied extensively by Lunardini and Smirnov 
[212]. To interpret the data in terms of neutrino oscillation, first we have to calculate the 
survival probability of u e ,p. It is important to note that neutrinos detected Kamiokande 
II and 1MB have different survival probabilities due to the different positions of the two 
detectors on the earth. Fig. 68 shows the permutation factors (1 — p) as a function 
of the neutrino energy at Kamiokande II, 1MB and Baksan. Here the H-resonance 
was assumed to be perfectly non-adiabatic and neutrino oscillation parameters were set 
as Am| 2 = 7.1 x l(T 5 eV 2 and sin 612 = 0.28. Although the average behavior is the 
same for all detectors the phase of oscillation is different. Further, due to the larger 
distance crossed by neutrinos inside the earth, for the 1MB detector the frequency of the 
oscillatory curve in the energy scale is twice as large as the frequency for Kamiokande 
II. 
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Figure 68. Permutation factor (1 — p) as a function of the neutrino energy at 
Kamiokande II, 1MB and Baksan [212]. 

In [212], they took the spectrum parameterization in Eq. (136). Then, given the 
parameters E 0e , E 0x , L e , L x , j3 e and (3 X , we can calculate the observed average energies 
and event numbers expected at Kamiokande II and 1MB. Fig. 69 shows the average 
energy of positrons (left) and event numbers (right) at the Kamiokande II and 1MB as 
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Figure 69. Left: Average energy of positrons, e 1 , in the detectors Kamiokande II and 
1MB as a function of the average energy of the original v e , E og , and different values 
of Eox/Eo e and L x /L e . Right: The predicted numbers of events at Kamiokande II 
and 1MB as a function of Eq £ . Here L e = 5.3 x 10 52 erg was used. In both figures, 
the horizontal lines represent the experimental results with the la error and the solid, 
dashed and dotted-dashed lines correspond to L x /L e = 1, 1.5, 0.667 respectively. Both 
figures are from [212]. 



functions of the average energy of the original v e , E QS , and different values of E 0x /E 0e 
and L x /L e . Neutrino oscillation parameters are set to the same value as in Fig. 68. 

Let us consider first the no-oscillation case. From the left of Fig. 69, we obtain the 
v e average energy: from Kamiokande II data, 

= 8.7±0.9MeV, (214) 

and from 1MB data, 

E™ B = 14.7 ± 1.9MeV. (215) 

Thus, the 1MB result is more than 3a above the Kamiokande II result. According to 
the right of Fig. 69, the expected event numbers at Kamiokande II and 1MB for the 
energies E™ and E™ B and L e = 5.3 x 10 52 erg are Nf 2 = 7.8 ± 1.5 and iV e IMB = 10±|, 
respectively. Therefore, to reproduce the observed event numbers at Kamiokande II and 
1MB, 

Lf = 8.2 x 10 52 erg, L : e MB = 4.2 x 10 52 erg, (216) 

are required, respectively. Thus, the 1MB signal implies about 2 times higher average 
energy and 2 times smaller luminosity in comparison with Kamiokande II. 

Neutrino oscillation improves the agreement between the Kamiokande II and 1MB. 
For Eo x /Eo e = 1.6 and L x /L e = 1, we obtain, 

E™ = 6.8 ± 0.8 MeV, E™ B = 10.3 ± 1.7 MeV, (217) 
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which now agree in 2a level. Actually, a "concordance" model with 

E 0e = 8MeV, E 0x = 12.8MeV, 
L e = L x = 8 x 10 52 erg, 

Pe = Px = 4, (218) 

gives the best fit to the all available data with x 2 = H-0 while the best-fit no-oscillation 
model with 

E 0e = HMeV, L e = 5.3 x 10 52 erg, (3 e = 4, (219) 

gives x 2 — 16.2. Thus neutrino oscillation leads to a certain improvement of the global 
fit of the data. The improvement requires lower average energy of the original v e 
spectrum and larger v e luminosity. The combination of smaller average energy and 
larger luminosity corresponds to a larger radius of the neutrinosphere: R 11S oc E^L l J 2 . 
It follows that in the concordance model R ns is about 2.4 times larger than in the 
no-oscillation model, which gives R ns = (20 — 30)km [206]. 

The concordance model gives rather small average energies of u e and v x compared 
to those predicted by numerical simulations. The situation becomes even worse 
with adiabatic H-resonance, where most of the observed neutrinos were originally 
v x s. Because u x s are expected to have larger average energy than v e , E 0x must be 
substantially smaller than that predicted by numerical simulations. Anyway, since the 
event numbers of SN1987A neutrinos would be too small to make a definitive conclusion. 
We are looking forward to seeing the next galactic supernova. 

4-5. Neutrino detectors 

It took 25 years to prove the existence of neutrinos experimentally since Pauli predicted 
theoretically in 1931. It took another 30 years to use neutrinos as tools in particle 
physics and astrophysics. This is because neutrinos have very weak interaction with 
other particles. Since the cross section of neutrino reaction is typically 10~ 40 cm 2 , mean 
free path of neutrinos in water is about 10 16 cm. This is why the current neutrino 
detectors are huge in volume. In this section we explain various neutrino detectors and 
their detection principles. 

4-5.1. water, heavy water and ice Neutrino interactions in water produces a relativistic 
charged particle which is mostly an electron or positron. If the velocity of the charged 
particle is larger than velocity of light in water, the charged particle emits Cherenkov 
light. This Cherenkov light is the signal of water Cherenkov detector. There are many 
neutrino detectors which utilize this method, starting from Kamiokande and 1MB to the 
current experiment SuperKamiokande (SK), Sudbury Neutrino Observatory (SNO) and 
IceCube (AMANDA). They are based on the same detection strategy but have different 
features which come from different target media: SuperKamiokande, SNO and IceCube 
use pure water, heavy water and antarctic ice, respectively. 
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neutrino reactions in water Cross sections of neutrino reactions in water at low 
energies, s < rri^, where s is the center-of-mass energy and m w is the W boson mass, 
are given by, 

( Ey 



a(v e e -> v e e) 
a{p e e -> V e e) 
a(v x e -> i/ x e) 
a(u x e -> P x e) 



rr(z/ e 16 ^e" 16 N) 
and shown in Fig. 70. 
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Figure 70. Cross sections of 
neutrino interactions in wa- 
ter. The dominant interac- 
tion is D e p — ► e+n at every 
energy. 
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Figure 71. Angular distri- 
bution of scatter electron in 
v e e~ — > ^ e e~ with £?„ e = 5 
and 10 MeV. 



Among these, the inverse beta decay (225) has the largest cross section at all 
energies so that most of the events at water Cherenkov come from this reaction. Because 
the emitted positron has an almost isotropic distribution for low-energy neutrinos like 
supernova neutrinos, it is difficult to know the direction of the incident neutrino. 

On the other hand, electron scattering event gives us information of the direction 
of the incident neutrino, although the cross section is much smaller than that of the 
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inverse beta decay. The angle between the scattered electron and the incident neutrino 
is given by, 



where m e and E e are the electron mass and electron energy, respectively. Angular 



should be noted that electron scattering occurs for all the flavors and they cannot be 
distinguished at low energies (s <C m^), although the cross sections are different. The 
u e e~ cross section is larger than those of the other flavors because electron scattering 
of v e is contributed from the charged current interaction as well as the neutral current 
interaction. 

SuperKamiokande The SuperKamiokande [306] detector is a cylindrical 50,000 ton 
water Cherenkov detector located at the Kamioka mine in Japan. It lies 1,000 m 
underneath the top of Mt. Ikenoyama, (i.e. 2,700 m water equivalent underground), 
resulting in a cosmic ray muon rate of 2.2 Hz, a reduction of 10.5 compared to the 
rate at the surface. As a supernova neutrino detector, SK has a fiducial volume of 
32,000 ton. The detector is optically separated into two regions, the inner and outer 
detectors. The inner detector of the SuperKamiokande-I detector, which operated 
from April 1996 to July 2001, was instrumented with 11,146 50-cm diameter inward 
facing photomultiplier tubes (PMTs) which provide 40% photocathode coverage. This 
photocathode coverage made it possible to detect low energy electrons down to ~ 5 MeV. 
Also, SK uses the anti-counter which surrounds the inner detector detect and remove 
events due to cosmic-ray muons. The direction of a charged particle is reconstructed 
using the directionality of the Cherenkov light. Angular resolution is about 25 degree 
for a 10 MeV electron. Energy resolution a for low-energy neutrinos is given by, 



As discussed in section 3.4, SK has been playing a central role in the field of neutrino 
oscillation experiment by observing neutrinos from the sun, atmosphere and accelerator. 
It is also expected to give tremendous information on supernova if it occurs in the future. 
For a detailed description of SK detector, see [309]. 

Due to an unfortunate accident in 2001, 60% of the PMTs were destroyed and 
the observation was interrupted for a while. However, the observation was restarted 
by redistributing the survived PMTs and SK is expected to resume normal observation 
with the original number of PMTs in 2006. 

IceCube IceCube [147] is an extended experiment of AMANDA [6, 7] and consists of 
an array of 4800 optical modules on 80 strings, regularly spaced by 125 m in antarctic 
km 3 ice. It covers an area of approximately 1km 2 , with the optical modules at depths of 
1.4 to 2.4 km below surface. Each string carries 60 optical modules, vertically spaced by 
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v e e is plotted in Fig. 71. Here it 
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17 m. IceCube is primarily designed to observe high-energy neutrinos (E > ITeV) from 
astrophysical sources. In order to reach the large volume needed to detect the expected 
small fluxes at high energies, the density of optical modules must be too sparse to 
measure low-energy neutrinos such as solar neutrinos. However, it is expected that 
IceCube can detect a supernova neutrino burst because the Cherenkov glow of the ice 
can be identified as time-correlated noise among all phototubes. The observed quantity 
is the number of Cherenkov photons caused by supernova neutrinos as a function of time. 
Thus, although IceCube cannot identify individual neutrino, it can measure neutrino 
luminosity of supernova. For a detailed description of IceCube, see [148]. Supernova 
neutrino detection at IceCube was studied in detail in [116, 76]. 
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Figure 72. Cross sections of neutrino interactions with deuteron. Solid lines are 
neutrinos and dashed lines are antineutrinos. 

Sudbury Neutrino Observatory Sudbury Neutrino Observatory (SNO) [293] is located 
in a large cavity excavated at the 2,039 m level (6,000 m water equivalent) in the 
Creighton mine near Sudbury. The 1,000 tons of heavy water are contained in an 
acrylic vessel surrounded by a light water shield. 

SNO's uniqueness is the use of heavy water as its neutrino detection medium. 
Neutrinos interact in heavy water in two additional ways not possible in ordinary water. 
One is through the charged current interaction: 

v e + d^ e~ +p + p (E th = 1.44MeV), (230) 

v e + d^ e + + n + n (E th = 4.03MeV), (231) 

whose cross sections are plotted in Fig. 72. These neutrino absorption reactions can 
only happen if the neutrino is an electron neutrino or electron anti-neutrino at low- 
energies (s < mV). Thus the neutrino absorption reaction exclusively counts electron 
neutrinos and electron anti-neutrinos. 
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Another is the deuteron breakup reaction: 

v + d^ v + n + p (E th = 2.22MeV). (232) 

This reaction occurs through the neutral current interaction and it occurs with equal 
probability for all neutrino flavors. In this reaction, no new charged particle is created 
and the free neutron cannot by itself create Cerenkov light. However, after scattering off 
of the nuclei in the heavy water it is eventually captured by another deuteron, creating 
a tritium nucleus and releasing a high energy 7 ray. This 7 ray then scatters an electron 
in the heavy water and it is this secondary electron which creates the Cerenkov light. 

SNO played a critical role in solving the solar neutrino problem with its ability to 
identify v e and neutral current events, as discussed in section 3.4.5. This feature will also 
have a great impact on observation of supernova neutrinos. For a detailed description 
of SNO, see [294]. 

4-5.2. scintillator As a high energy particle propagates in medium, it loses energy 
exciting electrons in the medium. A part of the deposited energy is emitted as 
scintillation photons. Scintillator detector uses the scintillation photons gnal 
to detect high energy particles. Here we summarize basic feature of scintillator detector 
as a neutrino detector. 

As to neutrino detector, liquid hydrocarbon is often used as a medium. In this 
medium, neutrinos interact with electrons, protons and carbon nuclei. Dominant 
contribution to the events comes from inverse beta decay, 

v e + p -> e + + n, (233) 

where the positron has almost the same energy as that of the incident neutrino. Then 
the neutron is absorbed into proton in about 170/zsec, 

n + v -> d + 7 2 . 2 McV (E th = 1 .80MeV) , (234) 
and the emitted photon scatter an electron, which then emits Cherenkov photons. Thus 
the positron and the delayed 2.2 MeV photon are the signal of inverse beta decay. 

Since carbon nuclei are abundant in scintillator, the following reactions also 
contribute to the events, 

(i) neutral current interaction 

12 C + v -> 12 C* + v {E th = 15.11MeV), (235) 
i2 C * 12 C + 7i5 nMeVj (236) 

(ii) v e capture by charged current interaction 

12 C + v e -> 12 N + e~ (E th = 17.34MeV), (237) 
12 N -> 12 C + e+ + v e (n /2 = 11.00msec), (238) 

(iii) v e capture by charged current interaction 

12 C + v e -> 12 B + e + (E th = 14.39MeV), (239) 
12 B -> 12 C + e+ + v e (n /2 = 20.20msec). (240) 
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The cross section of these reactions are plotted in Fig. 73. Because each reaction has its 
unique signal, we can identify them including inverse beta decay. This is an advantage 
of scintillator detector. 




energy (MeV) 



Figure 73. Cross sections of neutrino interactions with 12 C. Solid lines are neutrinos 
and dashed lines are antineutrinos. 

KamLAND KamLAND (Kamioka Liquid scintillator Anti-Neutrino Detector) is 
located underneath Mt. Ikenoyama in Gifu prefecture in central Japan, where 
Kamiokande was once located. It is a 1,000 ton liquid scintillator which is composed 
of 80% dodecane and 20% pseudocumene, whose typical composition is C n H 2ri . The 
primary purpose of the KamLAND is to probe the LMA solution of the solar neutrino 
problem by observing reactor neutrinos from the entire Japanese nuclear power industry 
as discussed in section 3.4.3. 

Because of low background level, KamLAND is also suitable for observation of 
supernova neutrinos. In [28], supernova neutrino detection at KamLAND was studied. 

LVD LVD (Large Volume Detector) is located in the INFN Gran Sasso National 
Laboratory, Italy. It consists of an array of 840 liquid scintillator (C n H2 n +2 with 
(n) = 9.6) counters and the active scintillator mass is 1,000 ton. The main purpose of 
the project is detection of supernova neutrinos. For details of the detector, see [2, 283]. 
Analyses of supernova neutrinos expected to be observed at LVD were performed in 
[210, 320]. 
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5. Explosion Mechanism of Core-Collapse Supernovae 

5.1. Status of Spherical Models 

Although the gross physical conditions of core-collapse supernovae are understood as 
denoted in the section 2, recent numerical simulations assuming spherical symmetry, 
however, with the current input physics (neutrino interactions and the equations of state 
of dense matter) and with/without general relativity, do not yield successful explosions 
by the neutrino heating mechanism (see Figure 74) [264, 198, 331, 201]. 

The problem apparently stems from the weak-interacting natures of neutrinos. The 
neutrino heating occurs only inefficient for the successful explosions in the spherical 
models. Whether neutrinos succeed in reviving the stalled shock wave depends on the 
efficiency of the energy transfer to the postshock layer, which in turn increases with 
the neutrino luminosity. In fact, it was pointed out that the stalled shock revives 
leading to explosions for otherwise failed explosion models by enhancing the neutrino 
luminosity artificially by a few tens percents from the original value [155] (see Figure 
75 and compare ID/2.10 and ID/2.225). It was also found by the study of the static 
configurations after the shock-stagnation that for a given accretion rate there is a critical 
luminosity for the shock revival [49] (see Figure 76). Janka et al. reported that the 
manipulation of omitting the velocity dependent term (0(v/c) ~ 10%) in the neutrino 
transport equations increased the neutrino energy deposition in the heating region and 
was sufficient to covert a failed model into a exploding one [156]. These facts suggest 
that all we have to obtain for the successful explosions, is the relatively small amount 
of boost of the neutrino luminosity and energy from the values we have obtained in the 
failed explosion models. 

In the last couple of years, both numerics and microphysics have been developed. 
The former, in particular, has seen major progress [264, 198, 331, 201, 45]. Ever 
since Wilson first proposed the neutrino heating mechanism [358], the accurate 
treatment of neutrino transport has been an important task, however mainly due to 
the computational intensity, some approximations, such as the multi-group-flux-limited 
diffusion approximation as the most familiar example [55, 56, 234], have been employed 
in the ID spherical symmetric simulations. This situation has changed completely lately 
(see [61] for a complete set of references). A couple of groups [264, 198, 331, 201], have 
published the state-of-the-art direct solutions of the Boltzmann equation for neutrinos, 
some of them extended even to 2D computations in the multi-group-flux-limited diffusion 
approximation [45, 202]. Although they have still not found successful explosions, the 
importance of the accurate treatment of neutrino transfer should be never missed. 

The microphysics such as neutrino reaction rates and equations of state have also 
been studied in detail. Recently, shell-model calculations of nuclear properties [191] 
revealed that the treatment of the electron capture rates for various nuclei should 
be changed significantly from the previous ones [55]. In almost all the supernova 
simulations, the electron captures rates on nuclei were cut off above a few ~ 10 10 g cm -3 
because only the resonant Gamow- Taylor transitions could be treated in the average 
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Figure 74. Comparison of the radial position of the shock waves as a function of time 
obtained by the two independent groups. The thin and thick lines are based on the 
simulations of a 15M© progenitor star, performed by vertex code of the Garching 
group [264] and AGILE-BOLTZTRAN code of Oak Ridge-Basel group [199], respectively. 
It is shown that spherically symmetric models with standard microphysical input fail 
to explode by the delayed, neutrino-driven mechanism. This figure is taken from 
Liebendorfer et al [201]. 

heavy nuclei, which is calculated by the employed equation of state (EOS) [55]. Despite 
of the quantitative change of the electron capture rates and hence the lepton fraction 
at the central portion of the core, subsequent shock propagations were found to show 
no significant change in comparison with the previous studies due to cancellation effects 
[136] (see Figure 77). 

As for the EOS of dense matter, we have at least two kinds of EOS now available 
based on different realistic descriptions of nuclear interactions, namely EOS by Lattimer 
& Swesty (LS EOS) [187] and EOS by Shen et al (SHEN EOS) [287]. Sumiyoshi 
et al. performed ID adiabatic hydrodynamic simulations employing the two kinds of 
equations of state, respectively, and found that there does appear the difference in the 
chemical compositions between the EOS's during the infalling phase, however, which 
disappears in the later phases, and hence no significant differences such in the remnant 
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Figure 75. Shock positions as a function of time after core bounce taken from [155]. 
ID or 2D represents the one or two dimensional models. The numbers indicate the size 
of the neutrino luminosities in unit of 10 52 erg s _1 injected from the proto-neutron star. 
Given a neutrino luminosity of 2.1 x 10 52 erg s" 1 , it can be seen that the unsuccessful 
explosion model in the ID simulation turns to lead the successful explosion in the 
2D simulation. Generally, the shock can propagate to the outer regions in the 2D 
simulations due to the convection in the hot bubble, which boosts the neutrino-heating 
efficiency. 

masses and the explosion energies are obtained [303]. Currently they performed the 
general relativistic neutrino radiation hydrodynamic core-collapse simulations assuming 
spherically symmetric and investigated long-term postbounce evolution after core- 
bounce employing the two equations of state. They found that, for both EOSs, the 
core does not explode and the shock waves stall similarly in the first 100 milliseconds 
after core bounce (see Figure 78). 

As for the neutrino reaction rates, the effects of the inter-nucleon many body effects 
near nuclear density have been elaborately studied (see Figure 79 and [52, 53, 267, 369].) 
Inter-ion correlations work for suppressing the neutrino-nucleus elastic scattering 
[141, 149], however, have been pointed out to lead to non noticeable changes of the 
dynamic during core-collapse (see Figure 80, [60, 223]). Relatively small corrections to 
the standard neutrino interaction processes [55] such as the detailed reaction kinematics 
of nucleon thermal motions, recoil, and weak magnetism contributions [142] as well 
as nucleon bremsstrahlung (NN' ^ NN'uu) [118], pair-annihilation/creations among 
different neutrino flavors [46], such as v e + v e — > v^ T + v^ T (see Figure 81), and 
the quenching of the axial vector coupling constants in dense matter [62] have been 
partly or fully incorporated and their importance has been evaluated in the recent 
ID computations [266, 331, 45]. These nuclear corrections effectively work for the 
suppression of the opacities leading to the larger neutrino luminosities, which is able to 
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Figure 76. Critical curve for the shock-revival as a function of the neutrino luminosity 
(L Ve in unit of 10 52 erg s _1 ) emitting from the central protoneutron star and the 
accretion rate (M in unit of M Q ) through the stalled shock. The "Model Trajectory" 
in this figure is taken from the results obtained in the numerical simulations by Bruenn 
(1982) [58], in which no explosion was obtained. If the evolution line of a real core 
crosses the critical curve into the hatched region, a neutrino-driven supernova should 
begin. This figure is taken from Burrows & Goshy (1993) [49]. 



act in favor of enhancing the explosion [266]. 

However even with these sophistications, the successful explosion has not yet been 
found for the present [266, 45]. We may still missing some important microphysical 
processes if we are to obtain the successful explosion assuming spherical symmetry. The 
phase transitions with various possible geometrical structures from isolated nuclei to 
uniform nuclear matter (the so-called "nuclear pasta") have been elaborately studied 
from a nuclear physics point of view [346, 347]. Due to this non-uniformity of the nuclei, 
the neutrino opacities are shown to be lowered, which are conventionally estimated by a 
single heavy nucleus approximation [144]. The inelastic neutrino-nucleus interaction has 
been pointed out to play the same important role of thermalizing neutrinos as effective 
as the neutrino scattering on electrons [57]. These microphysical ingredients, which have 
not been routinely taken into the modern supernova simulations, should deserve further 
investigations. 

5.2. Multidimensional Aspects in Core-Collapse Supernovae 

Ever since SN1987A was observed, most researchers in this field think that the dynamics 
of supernova is aspherical. The following observations of SN1987A have been attainable 
because SN 1987A is the nearest supernova for us and thus have been observed by 
modern observational instruments. For the explanation of the observed shape of the 
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Figure 77. The electron fraction, entropy, density and velocity as functions of the 
enclosed mass at the beginning of bounce for a 15 M Q model. The thin line is a 
simulation using the conventional capture rates by Brucnn parametrization [55] while 
the thick line for a simulation using the new rates by Langanke Martinez-Pincdo (LMP) 
[189, 190]. This figure is taken from Hix et al. (2003) [136]. 



lightcurve, the unexpectedly early appearance of X, 7-ray emissions, and Doppler 
features of spectral lines, the existence of a large-scale mixing between the deep 
interior and the hydrogen envelope is suggested (for a recent review, see [252]). The 
expanding debris of SN1987A is directly confirmed to be globally asymmetric by images 
of Hubble Space Telescope (HST) [261] and its axis roughly aligns with the small 
axis of the rings [345] (see the left panel of Figure 82). In recent years, it is noted 
that the same features are drawn for other core-collapse supernovae (see [138] and 
references therein). Spectroporalimetry shows that substantial asymmetry is common 
in core-collapse supernovae, indicating bi-polar explosions with axis ratios up to ~ 2 
[345, 244, 137, 343, 344]. The degree of the asymmetry tends to increase with time 
when greater depths are observed [344, 203]. Both suggests that a connection of the 
asymmetries with the central engine. For core-collapse supernovae with good time and 
wavelength coverage, the orientation of the polarization vector tends to stray constant 
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Figure 78. Radial positions of shock waves in models employing SHEN EOS and 
Lattimer & Swesty EOS are shown by thick and thin lines, respectively, as a function 
of time after bounce. The evolutions at early and late times are displayed in left and 
right panels, respectively. This figure is taken from Sumiyoshi et al (2005) [304]. 

both in time and with wavelength, which suggest that there is a global symmetry axis 
in the ejecta [344, 203]. Two oppositely directed jets [107] with the ejected material in 
a toroidal structure [351] around the center have been observed in the remnant of Cas 
A supernova. Young neutron stars are observed with high space velocities (typically 
300 - 400 km/s [213, 204], with highest values greater than 1000 km/s [17]), which are 
most likely imparted to the neutron star by a kick at the moment of the explosion. 
Interestingly, the recent X-ray observations have shown the correlation between the 
direction of pulsar motions and the spin axis of their supernovae in Vela and Crab 
pulsars [127, 260] (see the middle and right panels of Figure 82). All these observational 
evidences can be naturally interpreted as an evidence that the inner portions of the 
explosion are strongly aspherical. 
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Figure 79. Effect of inter-nucleon interactions on the neutrino-nucleon scattering 
rates taken from [369]. Left and right panels show the contours of the suppression 
and enhancement factors, which is the ratio of the scattering rate with to without 
the corrections. The top and bottom panel shows the contribution from the density 
correlation function (R\, vector currents) and the spin-density correlation {R2, axial- 
vector currents), respectively. Note the difference of the density scale in the right and 
left panels. It can be seen that the scattering rate is suppressed due to the corrections 
in the high density regime (p > 10 14 g cm -3 ) (see left panels), while the vector current 
contribution (right top panel) is enhanced in the low density region (p < 10 14 g cm" 3 ) 
and in the vicinity of the liquid-gas phase transition regions (the dark regions in the 
right panels). The above features are generally common in case of the other values of 
the proton fraction, Y p . See for details Yamada and Toki (2000) [369]. 
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Figure 80. Effect of ion-ion correlations on neutrino- nucleus elastic scattering during 
core-collapse of a 25M Q progenitor star taken from [60] . Left panel shows the optical 
depth along a radial path from the stellar center to the surface of as a function of the 
neutrino energy when the central density is 10 14 g cm~ 3 . Right panel shows the lepton 
fraction profiles at the same time with the right panel. From the left panel, the optical 
depth of the lower-energy neutrinos is shown to be lowered due to the screening effect 
of nucleus (compare the lines labeled as A and B). On the other hand, this change 
leads to a smaller change of the lepton fraction <~ 0.015 in the central region (right 
panel) . 
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Figure 81. The pair-production rates by the process v e + v e — > v^ T + v^ r (top line) 
and e+ + e~ — > i/^ )T + i^ )T (bottom line) as a function of the degeneracy parameter 
of electron- type neutrinos (r] Ve ) and of electrons (r) e ), respectively. Here T = 12 MeV 
and T = is assumed. The new rate, namely, v e + v e —* i/^ )T + ^ )T , is shown 
to dominate over the electron-positron pair reaction, which was considered only the 
sources for producing i/ M)T and v^ T . This figure is taken from Buras et al. (2003) [46]. 



So far many physical ingredients have been suggested in order to produce such 
asymmetric explosions. We review each ingredient one by one in the following. 

5. 3. Roles of Convections and Hydrodynamical Instabilities 

The convections and hydrodynamical instabilities have been long supposed to be 
responsible for producing the aspherical explosions. We give a summary of past studies 
about them, paying attention to where and why the convections and instabilities are 
likely to occur in the supernova cores and their possible roles of producing the aspherical 
explosions. 

5.3.1. convection near and below the protoneutron star It has been widely recognized 
that the convection in the protoneutron star (PNS) could play a crucial role in enhancing 
the neutrino luminosities from the PNS. In fact, Wilson and his collaborators obtained 
the exploding models [359, 360], in which a neutron-finger convection, which will be 
stated below, in a PNS was assumed, otherwise they did not see the successful explosions. 
This illuminates the importance of the convection in the PNS. 

After the shock wave stalls (t ~ 10 ms after core bounce), the outer parts of the PNS 
are convective unstable, because the deleptonization occurring in the shocked material 
outside the neutrino sphere produces a negative lepton gradient and the weakening of the 
prompt shock wave gives rise to a negative entropy gradient in the same region. Epstein 
first proposed that these two factors are an sufficient condition for the so-called prompt 
convections to occur [86]. In fact, the existence of these instabilities were demonstrated 
in most hydrodynamic simulations ([48, 50] see, however, [59]). 
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Figure 82. Observational evidences for asphericity in core-collapse supernovae. Left 
panel shows the inner debris of the supernova ejecta with an axis ratio of ~ 2 (inner 
red region) and the ring around it produced in the presupernova era (taken from [345]). 
The middle and right panel shows the nebulae of Vela and Crab pulsars, respectively 
(taken from [260] and [349]). 

After this prompt convection, the so-called "Neutron finger" instability develops 
inside the PNS [360]. This instability is driven in the presence of a positive entropy 
gradient and a negative lepton gradient, both of which are mostly satisfied inside the 
PNS. In the long-run 2D hydrodynamical simulations in the PNS, Keil et al. found 
that the neutrino luminosities increase ~ 50% due to the convection at times later than 
200 — 300 ms after core bounce (see Figure 83), which are expected to crucial for reviving 
the stalled shock wave [172]. 

It is noted that their neutrino transport was coupled only to the radial directions 
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Figure 83. Convection inside the protoneutron star. Top left and right panels show 
the absolute values of the velocity for two instants (about 0.5s (left) and 1 s (right) 
obtained in the 2D simulations of [172]. The growth of the convective region can be 
seen. Bottom panel shows the time evolutions of neutrino luminosities L u and mean 
energies of v e and v e for the model in the top panels without ("ID"; dotted) and with 
convection ( "2D" : solid) . Significant rise in the neutrino luminosities and energies can 
be clearly seen. Top panels are taken from Janka and Keil [151], and the bottom panel 
is taken from Keil et al [172]. 
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and thus suppresses the neutrino transport in the angular directions, essentially 
underestimating the stabilizing effect of the neutrino transport. On the other hand, 
Mezzacappa et al found only mild convective activity in the region near the neutrino 
sphere [235]. Since the neutrino transport in Mezzacappa et al [235] was assumed to be 
spherical, it might result in overestimating the stabilizing effect [235] (see Figure 84). 
The recent two-dimensional numerical simulation shows that the PNS convection really 
does occur, however, is not important for boosting the neutrino luminosity, because the 
convectively active layer is formed rather deep inside the PNS and is surrounded by 
a convectively stable shell [45]. However, one may argue that the conclusion may be 
still subject to change because their neutrino transport, albeit in the state-of-the-art 
manner, is not fully spatially two-dimensional, and hence could not reproduce all of the 
properties of the convective flows. 

5.3.2. convection in the hot-bubble regions In the hot bubble, in which the neutrino 
heating dominates over the neutrino cooling (the regions between R g and R s in Figure 
8), convections are expected to occur by the negative entropy gradient. In fact, a 
dynamical overturn between the hot, neutrino- heated, rising materials above gain radius 
and the cold postshock matter beneath has been demonstrated in the two-dimensional 
[128, 129, 50, 155, 269] or the three-dimensional [99] numerical simulations. As clearly 
demonstrated in [155], these convections indeed aid the shock revival and can lead to 
explosions in the case where the spherical models fail. 

In most of the preceding studies, it is noted that the neutrino luminosities from the 
protoneutron star are given by hand at some inner boundary (the so-called lightbulb 
approximation). In order to boost the neutrino-heating efficiency by the convection 
in the hot bubbles, a sufficient deposition of the neutrino luminosity from the hot 
protoneutron star is required, which unfortunately no previously published models with 
the neutrino transport in the whole regions have failed to reproduce. 

5.3.3. instability induced by the non-radial oscillation of the stalled shock waves It is 
recently pointed out that the stalled shock waves are subject to low-mode (I = 1 or 2) 
aspherical oscillations [38]. Interestingly, it is not due to the convection induced by the 
negative gradients of entropy or lepton fraction but due to the so-called vorticity-acoustic 
cycle. In the cycle, the vorticity perturbations given at the stalled shocks propagate 
inwards and reflected at some inner boundary, which is presumably the surface of the 
protoneutron star, producing the acoustic waves, and propagate outwards and create 
new vorticity perturbations when reaching the stalled shock waves. This closed cycle 
amplifies the growth of the aspherical oscillations (see Figure 85). This cycle was first 
discovered in the context of stability analysis of the accretion disk around the black holes 
[94, 95]. In the simulations by Blondin et al. [38], the oblique shock waves are found 
to feed vorticity in the postshock region and lead to growing turbulence. Although the 
model computations by them are not so realistic in the sense that no neutrino transport 
was included, which should affect the growth of the pressure waves predominantly by 
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Figure 84. Two-dimensional plots showing the entropy distribution of the 15 Mq 
model in a simulation "without" neutrino transport (top panel) and with neutrino 
transport (down panel). Taken into account the neutrino transport, albeit assuming 
spherical symmetry, it is seen that the instability in the protoneutron star is damped 
out within a short time. This figure is taken from Mezzacappa et al (1998) [235]. 

neutrino cooling, such instability seemingly appears in some recent realistic supernova 
simulations [157]. Combined with the convections and the hydrodynamical instabilities 
stated above, it is expected to have a good chance to obtain a successful explosion, 
because the neutrino heating mechanism seems very close to explosion as it is in the 
spherical collapse. However, the difficulties of multi-dimensional treatment of neutrino 
transport have hampered the definitive answer to this problem. 

Although some smoothed particle hydrodynamic (SPH) simulations have found 
explosions induced by the combination of neutrino heating and convection in the heating 
region [129, 99], there has been persistent concern with their approximate treatment of 
neutrino transfer. This situation may be changing recently. Some groups are preparing 
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Figure 85. Left panel shows the time evolution of the gas entropy (from left to 
right in this panel) , illustrating the growth of the instability induced by the aspherical 
perturbation, which is added at the stalled shock initially. Right panel shows the time 
evolution of the two-dimension simulation of the left panel in spacetime diagrams. 
The color at a given radius and time corresponds to the angle-average of the absolute 
value of the tangential velocity (right top) and absolute value of the deviation of 
pressure from the steady-state solution, AP/P (right bottom). The downflow of the 
tangential velocity is seen to be reflected at the inner boundary, producing the pressure 
perturbation propagating outwards (see the stripes in the top and down panels). Direct 
evidence of the vorticity-acoustic cycle are seen in these velocity plots. These figures 
are taken from Blondin et al [38]. 

for full spatially mult i- dimensional neutrino transport simulations [45, 202]. In these 
simulations, the dynamics of the whole core is computed with a code implemented with 
a Boltzmann solver. But still no successful explosions have been reported so far. We 
may have to look for alternatives to the convections and the hydrodynamical instabilities 
in order to get the successful explosions. 
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Figure 86. Initial angular velocity as a function of the radius obtained in the recent 
stellar evolution calculations [120, 121]. 15 (blue) and 20 (green) indicate the mass 
of the progenitor in unit of the solar mass (Mq). Much smaller angular velocities by 
models ml5b4 and m20b4 were evolved taking into account of the magnetic fields. The 
dotted red lines were constructed with the rotation law of Eq. (241) using the central 
angular velocity fig of E15 (4 rad s^ 1 ). Labeled number in the dotted red lines, such 
as 100, 500, represents the values of Rq in Eq. (241). This figure is taken from [258]. 



5.4- Roles of Rotation 

In addition to the hydrodynamical instabilities, rotation can produce the large 
asphericity in the supernova cores. It is well known that the progenitors of core-collapse 
supernovae are a rapid rotator on the main sequence [327]. Typical rotational velocity 
on the equator are on the order of 200 km/s, which is the significant fraction of their 
breakup rotational velocity [108]. Recent theoretical studies suggest a fast rotating core 
prior to the collapse [120], although this is not conclusive at all when the magnetic fields 
are taken into account [121, 122] (see Figure 86). Neutron stars, which are created in 
the aftermath of the gravitational core collapse of massive stars at the end of their lives, 
receive the rapid rotation, which are believed to be observed as pulsars. From the above 
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facts, rotation seems naturally to be taken into account in order not only to clarify 
the explosion mechanism, but also to explain the observed properties of core-collapse 
super novae. 

So far there have been some works devoted to the understanding of the effect of 
rotation upon the supernova explosion mechanism [239, 43, 274, 236, 366, 100, 101]. 
Among the studies, the systematic study of the rotational core-collapse, changing 
the initial angular momentum distributions and the strength has been done (see, for 
example, [366, 175, 258]). In the following, we show how the hydrodynamics in rotational 
core-collapse deviates from the one in spherical symmetry (section 2.3). 

5.4-1- hydrodynamics in rotational core-collapse The story of rotating core is not 
greatly different from the standard picture from the core bounce to the explosion as 
stated in section 2. In addition to the gravity and pressure gradients, one has to take 
into account the effect of centrifugal forces. In Figure 87, the time evolution of the 
shock wave in a rotational core-collapse simulation [175] is shown. The initial angular 
velocity profile of this model is based on the rotational progenitor model by [120], which 
is approximately fitted by the following way, 

«(0 = n x - J^, (241) 

where Q(r) is an angular velocity, r is a radius, and fl , R are model constants taking 
f2 = 4 rad s^ 1 and R = 1000 km (see Figure 86 and A1000). As a sideremark we 
note that r can be also interpreted as the distance from the rotational axis (the so- 
called cylindrical rotation). At first glance, the rotation profile seems more natural due 
to the Poincare and Wavre theorem [327], however, the shellular rotation profile, in 
which the angular velocity is constant with radius, is also pointed out to be natural 
because the horizontal turbulence is likely to be much stronger than the vertical one 
during stellar evolutions with rotation (for an elaborate description of the presupernova 
models including rotation, the reader is referred to [126]). Since even the-state-of-the 
art rotational progenitor models are based on the spherically symmetric models, it may 
not be decided yet which rotation profile is correct. Considering the uncertainties of the 
progenitor models, the rotating pre-collapse models so far have been made by changing 
parametrically not only the rotational profile (shellular or cylindrical), but also the 
rotational velocity and the degree of the differential rotation, and then adding them to 
the spherically symmetric progenitor models. 

Now let us return to the discuss of the model in Figure 87 again. At the central 
density of 2.3 • 10 14 g cm -3 , the model bounces at the pole first. The bounce epoch is 
slowed down of the order of 10 msec than the one in absence of rotation due to the 
centrifugal forces acting against the gravitational pull. The occurrence of core bounce 
below the nuclear density (~ 3.0 x 10 14 g cur 3 ) is a general feature in case of rotational 
core-collapse. This can be understood from the fact that rotation acts like a gas with 
an adiabatic index of 7 = 5/3. (see, for example, [327, 236, 271]) Due to this additional 
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Figure 87. Time evolution of the shock wave near core bounce for a typical rotational 
core-collapse model (we call it model A), whose initial angular momentum distribution 
is based on the rotational progenitor model by [120]. The two snapshots show color- 
coded maps of entropy (ks) per nucleon together with the velocity fields. The left 
panel is for 2 msec after the core bounce when the shock wave begins to move upwards 
at the pole, the right panel is for 4 msec. In the case of weak differential rotation, 
it is found that the hot entropy blob is formed near the rotational axis and it moves 
upwards. Note the difference of the scale of the plots. These figures are taken from 
[175]. 

pressure support, whose gradient is steepest along the rotational axis, the core bounce 
occurs at the sub nuclear density at the pole first. 

After core bounce, the shock wave begins to propagate a little bit faster along the 
rotational axis supported by the buoyancy of the hot entropy blob (S > 9.5kB/nucleon) 
(see Figure 87). At this time R, which is the aspect ratio of the shock front, is 1.1, 
which is lowered to 1.0 by the shock stagnation (see the right panel of Figure 88). The 
initially prolate configuration (R = 1.1) is stretched in the direction of the equator 
after the shock is weakened at the pole by the neutrino energy loss and ram pressure 
of the infalling material and finally stalls in the iron core. Note that a large angular 
momentum tends to push the matter parallel to the equatorial plane. 

Depending on the total angular momentum and its distribution imposed initially 
on the iron core, the effects of rotation on the core dynamics have many variety. In 
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Figure 88. Final profiles of the representative models. They show color-coded contour 
plots of entropy (fcs) per nucleon. The initial value of T/| W| is 0.5% for the models of 
the left (model A) and central panels (model B), 1.5% for the model of the right panel 
(model C). Here T 7 /] W| is the ratio of rotational to gravitational energy. The model 
of the central panel has stronger differential rotation than that of the left panel. The 
value of i?o of model B (central panel) is taken to be 100 km while the initial value of 
r/|W| is the same as that of model A with Rq — 1000 km. 

Figure 88, the entropy distributions at the shock-stall for some representative models 
[175] are presented. In fact, a variety of the final profiles is immediately seen in the 
figure. As the initial rotation rate becomes larger, the shape of the stalled shock wave 
becomes more oblate due to the stronger centrifugal force (compare the left with right 
panel in Figure 88). If the initial rotation rate is the same, the shape of the stalled 
shock is found to be elongated in the direction of the rotational axis as the differential 
rotation becomes stronger (compare the left with the central panel of Figure 88). This is 
related with the production of high entropy blobs by core bounce. As mentioned above, 
the entropy blob formed near the rotational axis floats up parallel to the axis and then 
stalls for weak differential rotations. This makes the shock prolate at first. Then the 
matter distribution returns to be spherical or oblate due to the centrifugal forces. On 
the other hand, in the case of strong differential rotations, the shock wave formed first 
near the rotational axis hardly propagates and stalls very quickly. The high entropy 
blob begins to grow near the equatorial plane in this case. This then induces the flows 
towards the rotational axis. As a result, the final configuration becomes prolate. 

When rotation is taken into account, the shock wave can generally reach further out 
than the one without rotation (-R sta ii <~ 200km). However, it was pointed that rotation 
does not good to the prompt explosion [366]. This is because the centrifugal force tends 
to halt the core collapse, which then reduces the conversion of the gravitational energy 
to the kinetic energy as clearly seen from Figure 89. 

5.4-2. anisotropic neutrino radiation in rotational core-collapse Then the problem 
is whether rotation does good or harm to the neutrino-driven mechanism. K. Sato's 
group in University of Tokyo has been paying attention to the effect of rotation on the 
neutrino heating mechanism. First of all, Shimizu et al. (2001) [272] demonstrated that 
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Figure 89. Relation between the explosion energy and the initial total angular 
momentum taken from [366]. q in the x axis means the normalized total angular 
momentum defined by q = J/(2GM/c 2 ) with J, M being the angular momentum 
and the mass of the core, respectively. The solid and dashed lines correspond to the 
uniform and differential rotation models, respectively. "d=0.836" in the figure indicates 
that in both cases the employed equation of state is the same (see [366] for details). 
As the initial total angular momentum increases, the explosion energy monotonically 
decreases. 



anisotropic neutrino radiations induced by rotation may be able to enhance local heating 
rates near the rotational axis and trigger globally asymmetric explosions. The required 
anisotropy of the neutrino luminosity appears to be not very large (~ 3%). In their 
study, the anisotropy of neutrino heating was given by hand and rotation was not taken 
into account, either. Kotake et al. (2003) [175] demonstrated how large the anisotropy 
of neutrino radiation could be, based on the two dimensional rotational core collapse 
simulations from the onset of gravitational collapse of the core through the core bounce 
to the shock-stall . In their study, a tabulated EOS based on the relativistic mean field 
theory [287] was implemented and the electron captures and neutrino transport was 
approximated by the so-called leakage scheme. They not only estimated the anisotropy 
of neutrino luminosity but also calculated local heating rates based on that. In the 
following, we state the main results shortly. 

In the left panel of Figure 90, the neutrino spheres after the shock-stall (~ 50 
msec after core bounce) for the spherical and the rotating model (corresponding to 
model A in Figure 87) are presented. For the rotating model, it is found that the 
neutrino sphere forms deeper inside at the pole than for the spherical model. This is a 
result of the fact that the density is lower on the rotational axis in the rotation models 
than in the spherical model because the matter tends to move away from the axis due 
to the centrifugal force. The neutrino temperature profile on the neutrino sphere for 
the pair models is presented in the right panel of Figure 90. Note that the neutrino 
temperature is assumed to be equal to the matter temperature. It is seen from the figure 
that the temperature varies with the polar angle for the rotating model. The neutrino 
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Figure 90. Shapes of neutrino sphere (left panel), neutrino temperature vs. polar 
angle on the neutrino sphere (right panel) for the rotating and the spherical models. 
These figures are taken from [175] 



temperature is higher at the pole for the rotating model than for the spherical model. 
This can be understood from the fact that the neutrino sphere is formed deeper inside 
for the rotational model than for the spherical model, as mentioned above. 

Based on the above results, the heating rates of the charged-current interaction: 
u e + n — > p + e~ outside the neutrino sphere were estimated. The neutrino emission from 
each point on the neutrino sphere is assumed to be isotropic and take a Fermi-Dirac 
distribution with a vanishing chemical potential. For the details about the estimation, 
we refer readers to [175]. In Figures 91 and 92, the contour plots of the heating rate for 
the above two models and the neutrino heating rate along the rotational axis with that 
on the equatorial plane for the rotating model are presented, respectively. It is clearly 
seen from the figure that the neutrino heating occurs anisotropically and is stronger 
near the rotational axis for the rotation models. This is mainly because the neutrino 
temperatures at the rotational axis are higher than on the equatorial plane. In addition, 
the radius of the neutrino sphere tends to be smaller in the vicinity of the rotational 
axis. As a result, the solid angle of the neutrino sphere is larger seen from the rotational 
axis. These two effects make the neutrino heating near the rotational axis more efficient. 
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Figure 91. Heating rates outside the neutrino sphere for the spherical model (right) 
and the rotating model (left). The color scale is the logarithm of the heating rate 
(MeV nucleon -1 s _1 ). The neutrino sphere and the stalled shock are seen as the thick 
lines separating the bright color from the dark color region. Note that the value within 
the neutrino sphere is artificially modified to dark colors and has no physical meanings. 
These figures are taken from [175]. 
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Figure 92. The comparison of the heating rate Q + v (MeV/(nucleon • s)) along the 
rotational axis with that on the equator for the rotating model. The pole-to-equator 
ratio of the heating rate outside the neutrino sphere is about 3. These figures are taken 
from [175]. 

By performing the linear analysis for the convective stability in the corresponding 
models, it was found that the convective regions appear near the rotational axis (see 
[175] and [100]). This is because the gradient of the angular momentum is rather small 
near the axis and the stabilizing effect of rotation is reduced there. The neutrino heating 
enhanced near the rotational axis might lead to even stronger convection there later on. 
Then, the outcome will be a jet-like explosion as considered in [272, 221, 222] (see Figure 
93). Here it should be mentioned that such jet like explosions play important roles in 
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Figure 93. Possible consequence of anisotropic neutrino radiation from the central 
protoneutron star, with the surface temperature of 4.5 MeV with 10 % neutrino flux 
enhancement along the rotational axis treated as the lightbulb approximation. The 
color map shows the dimensionless entropy distribution and the velocity fields (left 
t = 82 msec and right t = 244 msec after the shock stall). This figure is taken from 
Madokoro et al. (2003) [221]. 
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Figure 94. Critical luminosity for the accretion shock revival as a function of the 
accretion rate. Solid line denotes the spherical model, and the other lines correspond 
to the rotating models with the different initial angular velocities (O = 0.03 rad s" 1 ; 
dotted line and 0.1 rads^ 1 ; dashed line). It can be seen that the rotation can lower 
the critical luminosities for a give accretion rate. These figures are taken from [370] . 

reproducing the synthesized abundance patterns of SN1987A [248, 249]. 

Yamasaki and Yamada (2005) [370] recently reported the steady accretion flows onto 
the protoneutron star with a standing shock, in which they investigated how rotation 
affects the critical luminosity required for the shock revival. Note that this study is the 
extension of the study by [49], in the sense that the former newly takes into account 
the effect of rotation. For a given mass accretion rate, they found that rotation does 
lower the critical luminosity than the one in the case of the spherically symmetric mass 
accretion (see Figure 94). This result is also in favor of rotation for producing the 
successful explosions. 

A series of SPH simulations in 2D or 3D of rotational core-collapse has been done by 
[100, 101]. They referred to the deformation of the neutrino sphere induced by rotation 
(see figure 95). However, little effect of neutrino anisotropy on the explosion was found. 
This is probably their models explode in the prompt-shock timescale. 

Buras et al. (2003) [45] reported the core-collapse simulations of a slowly rotating 
model, in which the state-of-the-art neutrino reactions are included with the multigroup 
neutrino transport along the radial rays. Although the shock wave reaches further out 
in the rotating model than the one in the spherically symmetric model, the shock wave 
in the rotating model is shown to stall in the iron core (compare "sl5r" (rotating model) 
and "sl5" (spherical model) in the right-down panel of Figure 96). 
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Figure 95. Density isosurface (10 11 g cm -3 (blue), 10 12 g cm -3 (red)) implying the 
shapes of the neutrino spheres obtained in the 3D rotational core-collapse simulations 
by Fryer and Warren (2004) [101] (at 45 msec after core bounce in model SN15BB-hr). 
It is shown that their shapes are deformed to be oblate and the aspect ratios of the 
neutrino spheres are about ~ 2. 

Very recently, Walder et al (2005) [342] reported the two-dimensional rotating core- 
collapse simulations, in which the multi-energy neutrino transport with the flux-limited 
diffusion approximation was employed, and calculated the anisotropies of the neutrino 
flux in the rotating cores (see Figure 97). The degree of the anisotropy obtained in their 
simulations is almost the same with the aforementioned results [175] when the initial 
rotational velocity of the core are close with each other. In addition, they pointed out 
that the degree of the anisotropy becomes much smaller if the core rotates rather slowly 
as suggested by the recent stellar evolution calculation [122] and thus concluded that 
the rot at ion- induced neutrino heating anisotropy could not be a pivotal factor in the 
supernova explosion mechanism. 

Here it should be mentioned that the lateral neutrino heating, which is essential 
for the anisotropic neutrino radiation to work, can neither treated very appropriately 
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Figure 96. Left is the snapshots of the stellar structure for a rotating model 
(ft = 0.5 rad s _1 imposed on 15M Q progenitor model) at 198 msec after core bounce. 
The left panels show the rotational velocity (top) and the fluctuations of entropy (in 
percent) versus the enclosed mass, emphasizing the conditions inside the neutron star. 
The right panels display the rotational velocity (top) and the entropy as functions of 
radius. The arrow indicate the velocity field and the white line marks the shock front. 
Right panel shows the shock radii (R s ) vs. postbounce time (i p b)- The 2D models are 
compared to the corresponding ID simulations (thin lines). These figures are taken 
from Buras et al (2003) [45]. 



by the neutrino transport along the radial rays [45] nor by the flux-limited diffusion 
approximation [342]. Especially the diffusion approximation could lead to the 
underestimation of the anisotropy ratios. The fully spatially multidimensional radiation- 
hydrodynamic simulations seem required to give us the answer whether the anisotropic 
neutrino radiation does really help the explosion. 

5.4-3. convection in rotating protoneutron stars As in the non-rotating cases, rotation 
should affect the convective unstable regions in the protoneutron stars (PNSs). Janka 
and Keil (1997) [151] performed the two-dimensional hydrodynamic simulation in the 
PNS treating the neutrino transport in the grey flux-limited diffusion approximation 
and pointed out that convections occurs only close to the equator (see the left panel of 
Figure 98). On the other hand, a recent study by Miralles et al. (2004) [245] suggests 
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Figure 97. Integrated net gain (in unit of erg g _1 s _1 ) in the rotational core-collapse 
simulations with multigroup neutrino transport with a flux limitter. fio i n each panel 
shows the initial angular velocity of the core. The inner region 600 km on a side is 
shown. With increasing initial angular velocity, the heating rate is shown to be more 
and more concentrated along the rotational axis. This figure is taken from Walder et 
al (2005) [342]. 

that the convective unstable regions are formed along the direction of the rotational 
axis preferentially (the right panel of Figure 98). Note that the computations by [245] 
are based on the linear stability and only valid in a steady-state situation, which is not 
always satisfied in the PNS. The results between the two studies, namely, the appearance 
of the convective regions seems to change with the time evolution of the PNSs [245]. 
The effects of the heating anisotropy produced by the convections inside the rotating 
PNSs on the explosion mechanism are still an open question. 

As mentioned, the recent evolution models by [122] suggest that the transport of 
angular momentum during the quasi-static evolutionary phase of the progenitor deprives 
the core of substantial fraction of its angular momentum, particularly when the magnetic 
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Figure 98. Convection inside the rotating protoneutron stars. In the left panel, 
the contour of the absolute value of the fluid velocity evolved about 750 msec after 
core bounce evolved by a numerical simulation [151] is shown. From the panel, the 
convection is shown to be suppressed near the rotational axis (vertical) and develop 
strongly near the equatorial plane. In the right panel, the direction, to which the 
convective motions are likely to occur at 500 msec after core bounce, is shown, which 
is obtained by a stability analysis [245]. The convection is more effective in the polar 
region. The left and right panel is taken from Janka & Keil (1997) [151] and Miralles 
et al (2004) [245], respectively. 



torque is taken into account [273]. If this is really the case, the rotation will play no 
significant role in dynamics of core-collapse as shown by [45, 242, 342]. 
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5. 5. Roles of Magnetic Fields 

Another possible cause for the asphericity of supernovae may be magnetic fields. Soon 
after the discovery of pulsars, which are the magnetized rotating neutron stars, the role 
of rotation plus magnetic fields in the supernova explosion was scrutinized [257, 37]. 
Because of the magnetic flux conservations, a seed magnetic field in the stellar core can 
grow significantly during core-collapse. In addition, winding of field lines due to the 
differential rotation, which is natural after the collapse of a spinning core, can further 
amplify the toroidal field component. If the magnetic pressure becomes comparable 
to the thermal pressure, magnetohydrodynamical forces can drive an explosion [238], 
and accelerate axial jets [246, 274, 195, 12]. Also magnetic buoyant instabilities could 
produce mass motions along the rotational axis. It should be mentioned that the 
necessary condition for the working of the above mechanisms is that the core should be 
very strongly magnetized initially. The magnetic field strength is required to be more 
than ~ 10 16 G if the magnetic stress is to be comparable to the matter pressure in the 
supernova core after core bounce. Since the canonical value for the pulsar, ~ 10 12 G, is 
negligibly small in terms of the above effects on the dynamics of collapse, little attentions 
have been paid to the magnetic supernovae. 

However, it has been recently recognized these days that some neutron stars are 
indeed strongly magnetized as high as B > 10 14 G. The strong dipolar magnetic fields 
at the surface can be estimated as follows, 



^dipoic ~ 3.0 x 10 14 W [G], (242) 

which is distinct from the normal radio pulsars by their long periods (P) and high- 
period derivatives (P). Some of them are soft gamma repeaters (SGRs) and others are 
anomalous X-ray pulsars (AXPs) or high magnetic field radio pulsars (HBPs) [371, 115], 
which are collectively referred to as "magnetars" [84] (see Figure 99 and [361] for review 
of the current observations). Although they are supposed to be a minor subgroup of 
neutron stars and the field strength might be lowered to the normal pulsars [146], these 
situations do revive the study of the magnetic supernova again [368, 178, 326, 268, 179]. 



5.5.1. effect of magnetic fields on the prompt shock propagation As in the case of 
rotation [366], the first step to be investigated is the effect of magnetic fields (in 
combination with rotation) on the dynamics of the prompt propagation of a shock 
wave [368, 326]. 

The recent 2D magnetohydrodynamic simulations by Yamada & Sawai [368] and 
Takiwaki et al [326] showed that the combination of rapid rotation, > 0.1% and 

strong poloidal magnetic field, -E m /|H^| > 0.1%, leads in general to jet-like explosions. 
Here and ^ m /|W| represents the initial ratio of the rotational and magnetic 

energies to the gravitational energy, respectively. i? m /|W| > 0.1% roughly means that 
the iron core has B > 10 12 G before the onset of gravitational core-collapse. In the 
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Figure 99. B p — P diagram of radio pulsars and known magnetars taken from [371]. 
The number of the magnetars discovered so far is about 10. See for details [371]. 



following, we describe the magnetohydrodynamics features in such strongly magnetized 
and rapidly rotating models. 

The typical time evolutions are presented in Figure 100. The initial and 
^m/lWl imposed initially on the model is 1.5 % and 1.0 %, respectively. The initial 
field configuration is assumed to be parallel to the rotational axis with B = 5.8 x 10 12 G. 
After core bounce, the fast magnetohydrodynamic shock is launched in the direction of 
the rotational axis. From Figure 101, it can be seen that a seed magnetic field grow 
significantly during core-collapse due to the magnetic flux conservation (see the right 
panel of Figure 101 from A to C). Before core bounce (point C), the poloidal magnetic 
fields dominate over the toroidal ones. After core bounce, winding of field lines due 
to the differential rotation near the surface of the inner core can drastically amplify 
the toroidal components (see the right panel of Figure 101 after point C). In Figure 
102, the properties of the shock wave propagating further out from the central region 
are presented. From the top left panel, the shock wave produced after core bounce 
becomes magneto-driven, in the sense that the magnetic pressure becomes as strong as 
the matter pressure behind the shock wave (see the top right panel of Figure 102). Note 
that the magnetic pressure is much smaller than the matter pressure in the unshocked 
inner core due to its high density, however, can be much larger in the distant region 
from the center, because the matter density drops much steeper than the magnetic fields 
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which is almost constant along the rotational axis. The magneto-driven shock wave is 

B 2 

shown to be collimated along the rotational axis because the the hoop stress F hoop = 
which collimates the shock wave, is dominant over the gradient of the magnetic pressure, 

dB 2 

Fm&g — \~qx~i acting to expand the shock wave (see the bottom panel of Figure 102). 
Here X represents the distance from the rotational axis. These are the reasons for 
producing the jet-like explosion. 

Furthermore it was pointed out that the explosion energy increases with the initial 
magnetic fields strength, on the other hand, monotonically decreases with the the initial 
rotation rate [326]. This is because the collimated shock wave requires relatively lower 
energy to expel the matter near the rotational axis as the initial field strength becomes 
larger. As a result, the stronger initial magnetic field is favorable to the robust explosion. 
Interestingly, the models with the smallest magnetic fields studied in [368, 326] still 
produced the jet-like explosion although it takes longer for winding up the fields to 
launch the jet. This might suggest that even much smaller magnetic fields could be 
amplified in the collapsed core and play an important role for explosion. 

5.5.2. possible mechanisms for producing the pulsar-kicks During core-collapse of such 
strongly magnetized models, the strength of the magnetic fields substantially exceeds 
the QED critical value, -Bqed = 4.4 x 10 13 (G), above which the neutrino reactions are 
affected by the parity- violating corrections to the weak interaction rates [143, 15, 16, 10]. 
Kotake et al. (2005) [178] estimated the corrections based on the results in the 2D 
MHD simulations and discussed its role for producing the pulsar kicks, which will be 
summarized shortly below (see also [185, 186, 184] for a review). 

In the top left panel of Figure 103, the configuration of the poloidal magnetic fields 
after the shock stall (~ 50 msec after core bounce) is presented. For the initial condition 
for the model, the strong poloidal magnetic fields of 2 x 10 12 G was imposed with the 
high angular velocity of 9 rad/s with a quadratic cutoff at 100 km radius in the iron 
core. It is seen from the panel that the poloidal magnetic fields are rather straight and 
parallel to the rotational axis in the regions near the rotational axis, on the other hand, 
bent in a complex manner, in the other central regions due to the convective motions 
after core bounce. The top right panel of Figure 103 shows the ratio of the neutrino 
heating rate corrected from the parity-violating effects, AQ^ B _^ to the heating rate 
without the corrections, Q^ B=0 . It is noted that the suppression or the enhancement 
of the heating rate through parity-violating effects is determined by the signs of the 
inner product of n • B, where n, B are the unit vector in the direction of the incoming 
neutrinos and along the magnetic field. If the product is positive (negative), then the 
corrections to the heating rates results in the suppression (enhancement) of the heating 
rate. Reflecting the configurations of the magnetic fields, it is seen that the values of 
the ratio become negative in almost all the regions (see also the bottom right panel of 
Figure 103 for clarity). This means that the heating rate is reduced by the magnetic 
fields than that without. Especially, this tendency is most remarkable in the regions 
near the rotational axis and the surface of the neutrino sphere (see the regions colored by 
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Figure 100. Velocity fields (upper panels) and the poloidal magnetic fields (lower 
panels, white lines) on top of the density contours (log 10 p). Left and right is for core 
bounce and for 10 msec after core bounce. These figures are taken from Yamada and 
Sawai (2004) [368]. 
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Figure 101. Left panel is the time evolution of the central density and right is the 
time evolution of mean poloidal (dotted) and toroidal (solid) magnetic fields for a 
typical MHD model. Comparing the panels, one can see that "C" represents the time 
of core bounce. After core bounce, the toroidal component dominates over the poloidal 
one. These figures are taken from Sawai et al. (2005) [268]. 



blue in the top right panel of Figure 103). This is because the magnetic fields are almost 
aligned and parallel to the rotational axis in the regions near the rotational axis. In the 
bottom left panel of Figure 103, the contour of the ratio in the 360 latitudinal degrees 
region of a star, is prepared in order to see the global asymmetry of the heating. Since 
their simulations assumed the equatorial symmetry, the above features in the northern 
part of the star become reverse for the southern part. As a result, it was found that 
the heating rate is reduced about ~ 0.5% in the vicinity of the north pole, on the other 
hand, enhanced about ~ 0.5% in the vicinity of the south pole in the computed model 



133 




II iTi 




JC -.Ul 



Figure 102. A jet like explosion obtained in the rapidly rotating and strongly 
magnetized model studied in Takiwaki et al (2005) [326]. Top left panel shows the 
contour of the entropy per nucleon near the shockbreak out of the iron core showing 
a jet like explosion. Top right panel shows the contour of the logarithm of the plasma 
beta, f3, which is the ratio of the magnetic to matter pressure, representing that the 
shock wave is magneto-driven, since (3 is greater than 1 just behind the shock wave. 
Bottom panel shows the contour of the ratio of the hoop stress to the gradient of the 
magnetic pressure, which demonstrates that the shock wave is collimated by the hoop 
stress near the rotational axis. 
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Figure 103. Various quantities for a model with the strong magnetic fields and 
rapid rotation initially (see the context for the details of the model). Top left panel 
shows contour of the logarithm of the poloidal magnetic fields (log [.Bp (G)]) with the 
magnetic field lines. Top right panel represents the ratio of the neutrino heating rate 
corrected from the parity-violating effects, AQ+ B ^ to the heating rate without the 
corrections, Q^ B=0 - Note in the panel that the values of the color scale are expressed in 
percentage and that the central black region represents the region inside the neutrino 
sphere. Bottom right panel is the same with the top right panel except that the bottom 
right panel only shows the regions with the negative values of the ratio. Thus the white 
region shows the regions with the positive values of the ratio. Bottom left panel shows 
the contour of the ratio in the whole region, which is prepared in order to see the 
global asymmetry of the heating induced by the strong magnetic fields. These figures 
are taken from Kotake et al (2005) [178]. 

in [178]. If the north/south asymmetry of the neutrino heating persists throughout the 
later phases, it is expected that the pulsar is likely to be kicked toward the north pole. 

There is an another class of the mechanism for producing the kicks, which relies on 
the cause of the asymmetric explosion as a result of the global convective instabilities 
[51, 112, 102] caused by the pre-collapse density inhomogeneities [31] (see Figures 104, 
105) or the local convective instabilities [155] formed after the onset of core-collapse. 
More recently, Scheck et al [269]. recently pointed out that the random velocity 
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perturbation of ~ 0.1 % added artificially at several milliseconds after bounce could lead 
to the global asymmetry of supernova explosion and cause the pulsar-kick (see Figure 
106). While the orientation of pulsar- kicks is stochastic in their models, which is truly 
consistent with the observation of the normal pulsars, the asymmetry of the neutrino 
heating in the strong magnetic fields mentioned above might predict the alignment of 
the magnetic axis and the kick velocity in highly magnetized neutron stars. 

5.5.3. effect of toroidal magnetic fields While most of the magnetohydrodynamic 
(MHD) simulations in the context of core-collapse supernovae choose poloidal magnetic 
fields as initial conditions [195, 37, 246, 238, 274, 12, 368], recent stellar evolution 
calculations show that toroidal magnetic fields may be much stronger than poloidal 
ones prior to collapse [121, 122, 273]. Motivated by this situation, Kotake et al (2004) 
investigated the models with predominantly toroidal magnetic fields, changing the 
strength of rotation and the toroidal fields systematically. 

The angular velocity profile for the model with the strongest toroidal magnetic 
fields in their simulations is given in the top left panel of Figure 107. The initial values 
of T/|VF| and _E m /|W^| are 0.5,0.1%, respectively, where _E m /|W^| represents magnetic 
to gravitational energy. In addition, the initial profiles of rotation and magnetic field 
are chosen to be cylindrical with strong differential rotation for this model. 

In the top left panel, the negative gradient of angular velocity, dVt/dX < 0, can 
be found, where Q is the angular velocity and X is the distance from the rotational 
axis. Such a region is known to be unstable to non-axisymmetric perturbations [26, 4], 
as will be discussed in subsections of 5.5.4 and 5.5.5. The characteristic time scale for 
the growth of the instability called the magnetorotational instability (MRI) is given as 
tmri = 47r|<if2/<ilogX|~ 1 . The top right panel of Figure 107 shows the contour of t M ri 
for the model. The typical time scale is found to be ~ O(10) ms near the rotational 
axis. This suggests that MRI induced by non-axisymmetric perturbations can grow on 
the prompt shock time scale. The field strengths in the protoneutron star become as 
high as ~ 10 16 G (see the bottom left panel of Figure 107), and the ratio of magnetic 
stress to matter pressure gets as high as 0.9 behind the shock wave (see the bottom 
right panel of Figure 107). 

As for the anisotropic neutrino radiation observed in purely rotating case 
(subsection 5.4.2), the feature is shown to be not changed significantly by the inclusion 
of very strong toroidal magnetic fields (~ 10 16 G in the protoneutron star). Combined 
with the anisotropic neutrino radiation that heats matter near the rotational axis more 
efficiently the growth of the instability is expected to further enhance heating near the 
axis. Furthermore, the magnetic pressure behind the collimated shock wave is as strong 
as the matter pressure in the vicinity of the rotational axis. From these results, one 
might speculate that the magnetar formation is accompanied by a jet-like explosion if 
it is formed in the magnetorotational collapse described above. 
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Figure 104. Density inhomogenties (left) produced during oxygen shell burning and 
their root-mean-square azimuthal averages (right). These figures are taken from Bazan 
and Arnett (1998) [31]. 
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Figure 105. Kick velocities (x, y, z directions) of the neutron star due to the density 
inhomogenities prior to collapse. The inhomogenities were parametrically imposed in 
two ways, namely on the oxygen shell only and the core oscillation, with the density 
variations of < 25%. After the oscillatory behaviors, which were pointed out to be due 
to the neutrino emission from material accreting onto the neutron star, the motions 
of the neutron star tend to damp in time with the kick velocity less than 200 km s" 1 . 
This figure is taken from Fryer (2004) [102]. 
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Figure 106. The global asymmetry of I = 1 mode and the kick velocities of the 
protoneutron star. Left panel shows the growth of I = 1 at 1 s after core bounce. In the 
right panel, the velocities of the neutron star (top) and the explosion energies (bottom) 
as a function of time measured from core bounce are given for some representative 
models (see [269] for details). At 1 s, it can be seen that the kick velocity becomes as 
high as ~ 600 km/s in one of the models wit the explosion energy of ~ 1 x 10 51 erg. 
These figures are taken from Scheck et al (2004) [269]. 
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Figure 107. Contour plots of various quantities for the model with the strongest 
magnetic field in the computations [177]. The top left panel shows the logarithm of 
the angular velocity (s _1 ). The top right panel represents the magnetorotationally 
unstable regions to non-axisymmetric perturbations. The plot shows the growth time 
scale (s) of MRI. Note in this panel that the regions with white represent stable regions 
against MRI. The bottom left panel shows the logarithm of magnetic field strength (G). 
The bottom right panel displays the logarithm of the ratio of magnetic stress to matter 
pressure in percentage. 
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Figure 108. Schematic figure prepared for the derivation of the growth condition of 
MM. When the radial displacement of Ar of the fluid element is imposed, the MRI 
develops if the sum of the gravity and the centrifugal force becomes larger than the 
magnetic tension (see text for details). 

5.5.4- foundations of magnetorotational instability So far in this section, we have 
concentrated on the models which have large magnetic fields prior to core-collapse. 
The most mysterious is the origin of this large magnetic field. The core might have 
strong magnetic fields prior to collapse already, although the evolution models indicate 
quite the contrary [122]. We may not need the strong magnetic field initially if the 
magnetorotational instability (MRI) sufficiently develops in the supernova core. This 
MRI was first discussed in the context of accretion disks by [26]. Before we mention 
the role of MRI in the supernova, we summarize, for convenience, the basic properties 
of MRI in the context of accretion disks around the black hole, which are common 
situations in the central part of the active galactic nuclei. 

First of all, let us derive the growth condition for the MRI. See a schematic Figure 
108, describing the small displacement of the magnetic fields in the accretion disk. Before 
the displacement, the magnetic fields with strength B are assumed to be uniform and 
parallel to the rotational axis of the disk (z axis). The displacement is assumed to occur 
at the radial distance r from the central region of density p. Furthermore, let us assume 
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that the disk is in a Kepler rotation, namely, the angular velocity of the disk yields, 



where M is the mass of the central object, presumably the black hole. 

As in the Figure 108, let us displace the fluid element Ar inwards. Since the ideal 
MHD approximation can be well satisfied in the accretion disks, the magnetic fields are 
frozen-in to the matter. Thus the field line also moves Ar inwards as in Figure 108. 
Since the angular velocity of the magnetic field is constant, the centrifugal force becomes 
smaller due to the displacement. The change of the centrifugal force can be written, 

5F Iot = (r - Ar) P n 2 . (244) 

On the other hand, the gravity in the radial direction increases, 



W„ = -^(l + 2£). (245) 

As a result, the net force in combination of the centrifugal force and the gravity becomes, 

5F tot = 5F TOt + 5F grav = -3 P n 2 Ar, (246) 

which acts to move the displaced field line inwards. If this force is stronger than the 
magnetic tension, which acts to put the displacement back, the instability develops. 
The magnetic tension is B 2 /r' , here r is the curvature radius of the magnetic field (see 
Figure 108). From a simple geometric calculation, one can obtain the magnetic tension, 

B 2 2B 2 Ar , nArrS 

6F ™ = - = W (247) 

where A is the wavelength of the perturbation. Thus the condition for the growth of the 
instability becomes, 

SF mag + SF top = - 3pfi 2 ) Ar < 0. (248) 

As a result, the wavelength of the perturbation satisfying the condition A > A c = 
4^2/3v A /Q becomes unstable, where va is the Alfven velocity defined as va = :^f=- 
This is the condition for the onset of the MagnetoRotational Instability (MRI). 
Growth rate of the MRI can be roughly estimated as 

VA 

^growth ~ (249) 
A c 

To investigate the dependence of the growth rate on the wavelength of the perturbation, 
one need more detailed analysis (see for a review [26]) than the one done above. In Figure 
109, the dispersion relation between the frequency (o>) and the wavenumber (k z : along 
the unperturbed magnetic field) of perturbations is shown. The above relation was 
found by the linear analysis by adding the axisymmetric perturbations to the Keplerian 
disk with uniform magnetic fields parallel to the rotational axis and with the assumption 
that the fluid is incompressible (see for detail [27]). The oscillation, which is marginally 
stable (00 = 0) in the absent of the magnetic fields, is shown to be unstable (uj < 0) in 
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Figure 109. The growth rate of MRI again axisymmetric perturbations in the 
incompressible and magnetized disk with a Kepler rotation. In the region with u> < 0, 
the MRI develops. This figure is taken from Balbus & Hawley (1991) [27]. 



the presence of the magnetic fields. The maximum growth rate (minimum of the curve 
in Figure 109) is shown to be, 

3. 



with 



k z v A 
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(250) 



(251) 



Since this is very rapid indeed, the amplitude with the unstable mode can be significantly 
amplified during several rotation periods. The MRI can amplify the fields even if the 
initial field is very weak. It should be noted that the growth rate is independent of 
the strength of the magnetic field. Instead, the growth rate of MRI is determined by 
the angular velocity (see Eq. (250)). This is the main characteristic of the MRI, which 
is distinct from other magnetic instabilities such as Parker instability. As the MRI 
develops, the initial poloidal magnetic fields are stretched to the toroidal ones, whose 
strength grows exponentially with time. It is noted that the MRI can also develop 
nonaxisymmetrically. The evolution of the magnetic energy due to the nonaxisymmetric 
MRI in the accretion disk threaded by the purely toroidal fields is shown in Figure 110. 
The MRI and its associated angular momentum transport have been paid great attention 
in the community of the accretion disks around the black hole. 
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Figure 110. The spatially averaged magnetic energy normalized by the gas pressure as 
a function of the rotation period of the accretion disk obtained in the 3D simulations 
of accretion disk threaded by the toroidal magnetic fields. /?o represents the initial 
value of the plasma beta imposed on the disk. It can be seen that the magnetic energy 
increases exponentially regardless of the initial (3 and reaches to the quasi-steady state 

with -fgas //"mag 

10. This figure is taken from Matsumoto et al (1999) [224]. 



5. 5. 5. possibility of the growth of magnetorotational instability in supernovae The MRI 
mentioned in the last subsection was first applied to the dynamics of supernova core by 
Akiyama et al (2003) [4]. Compared with the linear growth of the toroidal magnetic 
fields in case of field wrapping, the field is expected to grow exponentially due to MRI 
also in the supernova core. If it is true, we may need the only small seed magnetic 
fields. When MRI develops efficiently in the supernova core, the field might reach to 
the saturation strength, which is estimated to be as high as 

~4xl0»( £ iV'MMf ^r) n (252) 

VI x 10 13 g cm" 3 / VlO km/ VlOOO rad s' 1 J 1 J ' v ; 

where we employ the typical values near the surface of the protoneutron star (see 
Figure 111). Here the saturation field is determined by the condition that the toroidal 
component of the Alfven speed, v\ = ;^f=, comes into rough equipartition with the 
rotation velocity, v TOt = rQ, in analogy with the case in the accretion disk [26]. 

In such a case, the magnetic stresses generated by the MRI could be the origin 
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Figure 111. The angular velocity (left panel) and the saturation magnetic fields 
(right panel) as a function of radius taken from Akiyama et al (2003) [4]. It is seen 
from the left panel that at core bounce, the angular velocity decreases sharply with 
the radius near r ~ 10 km, where is near the surface of the protoneutron star. This 
negative gradient of the angular velocity (dil/dr < 0) is the criteria for the onset of 
the MRI (see [26] for its derivation). They pointed out that the magnetic fields could 
be amplified as high as 10 16 G in the later phases as shown from the right panel. 

of the viscous energy deposition [332]. By employing an a prescription for the viscous 
dissipation, which is often used in the study of accretion disk [270], the robust explosion 
was pointed out to be obtained because the MRI heating aids the neutrino heating (see 
Figure 112). 

Although there remains persistent concern of the treatment of the rotation and 
magnetic fields in the above ID model computations [4, 332], the obtained implications 
seem to be important and should be examined by the multidimensional MHD simulations 
[332]. Currently the 3D MHD simulations have just begun to be investigated (see, for 
example [200]). Very recently nucleosynthesis in the magnetic supernovae has been 
reported [251]. Although the study of magnetic supernovae is still in its infancy with 
respect to the treatment of the neutrino physics, this field seems to be blossoming with 
the recent developments both of the transport method [45, 202, 342] and the growing 
computing power. 
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Figure 112. Effect of viscous heating due to MM on the explosion mechanism taken 
from Thompson ct al. (2005) [332]. Left panel shows the energy deposition rate as a 
function of radius for a rotating model with the initial period of 2 seconds, including 
viscous heating via the MRI with a = 0.1. In this panel, total energy deposition rate 
(<Ztot = Qv + <Zmri, solid line), viscous energy deposition (9mri, dot-dashed line), and 
neutrino deposition rate from each species {v e short-dashed line, v e dotted line) are 
shown. The order of 10 % increase of the heating rate due to the viscous heating 
above the gain radius (R ~ 120 km) can be seen. This increase was pointed out to 
be sufficient to instigate the explosion. Right panel shows the various timescales of 
the same model in the left panel. th v , t~c u , T g -„, T a( jv, tmri, and t«j T q T , represent 
the timescales of neutrino heating, neutrino cooling, the advection of the infalling 
matter, the viscous heating of MRI, the total heating including the contributions of 
MRI, as a function of radius, respectively. R„ e and R g mark the position of the v e 
neutrinosphcrc and the gain radius, respectively. In the regions above the gain radius 
from 130km < R <^ 170 km, the net heating time scale becomes shorter than the 
advection timescale of the infalling matter as a result of the viscous heating. Otherwise 
in the slowly rotating models Jo > 5 s, the explosions were not obtained in their 
computations, because there are no regions, in which the condition, Tq v < r a d v , is 
satisfied. 
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6. Gravitational Waves from Core-Collapse Supernovae 

From this section, we review the studies about gravitational waves from core-collapse 
supernovae. As mentioned in section 1, the detection of the gravitational signal is 
important not only for the direct confirmation of general relativity but also for the 
understanding of the explosion mechanism supernovae themselves. In combination with 
neutrino signals mentioned in section 4, the gravitational wave will enable us to see 
directly the innermost part of an evolved star, where the key physics related to the 
explosion mechanism such as the angular momentum distribution and the equation of 
state are veiled. 

In this section, we use the convention that Latin indices run from 1 to 3, Greek 
from to 3, where is the time component of four- vectors and that partial derivatives 
d/dx^ of tensor T ai "' an with respect to a coordinate x a are denoted by T° 



,a- 



Rfiv ^9fjLvR — rA Tjxvi (253) 



6.1. Physical foundations 

Before we discuss the gravitational waves from core-collapse supernovae, we summarize 
the physical foundations of Gravitational Wave in Einstein's theory of gravity, for later 
convenience and completeness of this review (see, also [348, 330, 232, 247]). 
As well known, the Einstein equations, 

-g, v R - — 4 

express the relation between matter distribution in spacetime, on the right hand side 
representing the matter stress-energy tensor, and the curvature of spacetime on the left 
hand side, represented by the components of the Ricci tensor given by the contraction 
of R^ u = R a ' mv of the Riemann curvature tensor, and the scalar curvature R = R a a . 
The Riemann curvature tensor is connected to the metric g^ v of spacetime through the 
Chiristoffel symbols Pg : 

pa ._ -pa -pa , -per -pa -pa -pa (OKA\ 

n /3/iv = 1 0v,n ~ 1 0n,v "+" 1 /3u L aji ~ 1 j3fi L uvi l Z04 J 

with 

1 

If one defines the Einstein tensor as 

Gpv = R^ u — -^g^uR, (256) 
the Einstein equations can be compactly written, 

G^y = ^ 4 Tfj, u . (257) 

It is noted that the Einstein equations reduce to Poisson's equation, V 2 = 4nGp in 
the Newtonian limit. To solve the Einstein equations, which consists of 10 nonlinear 
partical derivative equations, has been a challenging problem. So far only some exact 
solutions, such as Schwarzshild and Kerr solutions etc, have been found in the very 
idealized physical situations. 



r ^ = n9 a ° Xg^,p + 9ap,v - 9fa a ). (255) 
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6.1.1. weak-field approximation The gravitational waves far from the source can be 
characterized as linear metric perturbations propagating on a flat backgroud. Taking a 
first order perturbation from the Minkowskian metric: i]^ = diag(l, — 1, — 1, — 1), the 
metric g^ v of spacetime can be decomposed as, 

9 in, = ?V + V" ( 258 ) 
with the demand that 

| VI < !■ ( 259 ) 

Roughly speaking, the weak field conditions can be satisfied when GMj (c 2 R) <C 1 with 
M and R being the characteristic mass and the length scale of the system. If we take 
M = M , the approximation is well valid in the distance r 3> GM & /c 2 ~ 1.5 km. 
To first order in h, the Ricci tensor becomes, 

i^r^-r^ + o^ 2 ), (260) 

and the affine connection is 

= \v Xp (h P u,, + h mu - V,p) + 0(h 2 ). (261) 

Since we treat the linearized equations up to the order 0(h) in the following, the raising 
and rising of all indices can be done using rj^, not g^ v , that is, 

, rfc ,, A ^,^_|_.JL etc , (262 ) 

because using g^ makes the equations non-linear again. 

Introducing Eq. (261) to (260), one can obtain the first-order Ricci tensor, 

R$ = -(DV - h\ M - h\ M + h\ tllu ), (263) 

here □ = r]^ u d J? dxV is the d'Alembert operator. Then the Einstein field equations (Eq. 
(253)) read 



16ttG, 



□V " h »M ~ h \M + = ~ T ^- ( 264 ) 

Since the above equations are invariant under arbitrary coordinate transformations, the 
solutions cannot be determined uniquely. Let us introduce the most general coordinate 
transformation in the following form, 

xf t ^x< t = xf i + e{x! i ). (265) 

Here it is noted that de^/dx" should be at most of the same order of magnitude as 
h^y not to violate the weak-field condition. Since the metric of spacetime in the new 
coordinate is given by 

/„„ dx'^dx" ,„ 

then 

fcV = h iw _ e n xr) x» _ e v ^ = h iw _ e w _ (267) 
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One can easily check that the new K^ v satisfy the linearized Einstein equations by 
introducing Eq. (267) to Eq. (264). 

The above property, namely gauge invariance of the field equations, is a nuisance 
when one actually solves the field equations. To circumvent this problem, one has only 
to choose the coordinate system. The most familiar and convenient choice is to work in 
a harmonic coordinate, such that, 

<TV = 0. (268) 

Using Eq .(261), one can obtain equivalently, 

If h^ u does not satisfy Eq (269), one can find the new h'^ u by performing the coordinate 
transformation (Eq. (265)) with e u satisfying the condition, 

□e, = KS* ~ \ W (270) 

It should be noted that there still remains the freedom of the coordinate transformation. 
For example, perform the coordinate transformation [x^ — > x' 11 = x^ + e' M (x)) with e v 
satisfying the following condition, 

Ue v = 0. (271) 

Then the condition in Eq. (269) is indeed satisfied. We return to this problem soon in 
subsection 6.1.2. 

Using the harmonic gauge condition (Eq. (269)) in Eq. (264), the field equations 
now read, 

UhT = ^—7^- (272) 

Finally one can find the physical formal solution in a form of the time-retarded Green 
function, 

V(t, x) = *£[ dV T ^ t ;J^ X ' ) . (273) 

6.1.2. wave solutions in vacuum Now we move on to consider the solution of the 
linearized Einstein equations in vacuum (T^ v = 0). Then Eq. (272) reduces to 

□V = °- (274) 

The well-known plane wave solution is 

hfu, = e MI/ exp(iA; A x A ). (275) 

Introducing the solution into Eq. (274) yields k\k x = 0, which means that the 
gravitational waves travel along the null geodesies at the speed of light. Since Eq. 
(275) should satisfy the harmonic condition (Eq. (269)), 

- \ k »< = °- (276) 
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Let us express the remained freedom of the coordinate transformation (Eq. (271)) in 
the following form, 

e»(x) =ic ( *exp(ik x x x ), (277) 
where c M is the constant vector. Introducing Eq. (277) to Eq. (267), one obtains 

= &iiv "I - k^c u + k u c^. (278) 
For instance, let us consider a wave traveling in the z direction with wave vector, 
k 1 = k 2 = 0, k° = k 3 = k > 0. (279) 
Using Eq. (278), we determine c M satisfying e' 00 ,e' 0i = 0, namely 

eoo = eoo + 2/c c = (280) 
e 'oi = e oi + koQ + kiC = 0. (281) 

Noting ko = —k, the component of c M becomes, 

eoo e i e 02 e 03 + kc , , 

Co = W Cl = T' C2 = T' C3 = "Hfe— (282) 

Since Eq. (276) is invariant under the gauge transformation (Eq. (278)), we omit ' in 
the following. Using Eq. (280 and 281), Eq. (276) becomes, 

ke 3^ - ^M e " + e 22 + e 33) = 0. (283) 
When we see the zero component of the above equation {p, = 0), 

en + e 22 + e 33 = 0, (284) 

because 

e 3 o = e 03 = 0. (285) 
Using Eq. (284), Eq. (283) means, 

e 3 z = 0. (286) 
Taking I = 3, Eq (284) becomes 

en + e 22 = 0. (28 7) 

Finally, the can be written in the following form, 

/ \ 

en ei2 
ei2 -en 
\ J 

So far, the spatial component of k^ is taken to be z axis, Eq. (288) indicates that 
with arbitrary k^ satisfies, 

eij 5 ij = (Tranceless), (289) 
tijk? = (Transeverce). (290) 



e^ 



(288) 



149 



y 

z 

/_ x 




Phase %S2 k 3xJ2 2k 



Figure 113. Schematic figure describing the changes of the proper distance due to 
the plus h + and cross h x modes of the gravitational wave, which propagates along 
the +z axis. Note that horizontal and vertical directions express the x and y axis, 
respectively. 

It is noted that the above choice of the gauge is the so-called transverse-traceless (TT) 
gauge. We will use it in the following and denote it by the superscript TT. Since e 0fl 
was originately correspond to the freedom of the coordinate choice, they have nothing 
to do with the physical freedom. After all, the nonvanishing two components of en and 
e22 remain. We call the true physical freedom as the gravitational wave. 



6.1.3. polarization of gravitational waves First of all, we pay attention to the two 
spatial components of the gravitational wave propagating +z direction, (hJJ in Eq. 
(288)), 

(K 

h x -h+ I . (291) 



TT 



\ 








Here the corresponding metric, which expresses the superposition of the two plane waves 
can be written in a form, 

ds 2 = -alt 2 + (1 + h + )alx 2 + (1 - h + )dy 2 + 2h x dxdy + dz 2 , (292) 

where h+ = h + (ct — z), h x = h x (ct — z). As clearly seen from the metric, the 
gravitational wave changes the proper distance in the plane (x — y plane) perpendicular 
to the propagating direction (z axis). When the gravitational wave with h+ (plus mode) 
propagates, the proper distance in the x axis becomes longer (shorter) and the one in 
the y axis becomes shorter (longer) (see the top panel of Figure 113). 
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Next, let us rotate the coordinate 7r/4 along the +z axis, 

X '.) = ( C ° S7r/ t Sin?r/ ' )( X ) = ( & X + V \ ) (293) 
1/ / \ -smvr/4 costt/4 J \ V J \ -^(-x + y) J 

Then the metric in Eq. (292) reads, 

ds 2 = -dt 2 + (1 + h x )dx' 2 + (1 - h x )dy' 2 + 2h+dxdy + dz 2 . (294) 

Thus one can see that the cross mode of the gravitational wave h x represents the change 
in the proper distance tilted 7r/4 with respect to the change formed by the plus mode (see 
the lower panel of Figure 113). In this way, gravitational waves act tidally, stretching 
and squeezing space in a quadrupole manner and thus object that they pass through. 



6.1.4- meaning of TT gauge In this subsection, we discuss why TT gauge is important 
and how we can extract the TT components from the arbitrary tensor hij. 

First of all, we should be cautious that the solution of Eq. (274) is not only the 
gravitational wave. For example, let us consider the static solution around the mass 
point M in the Newtonian limit, 

2GM\ 



9oo 



= ~ i 1 " » 9oj = 0, 9ij = 5 tJ . (295) 



For the distant observer r 3> 2GM/c , the 00 component of the metric can be written, 

<7oo = -l + e-, (296) 
r 

with e being 2GM/c 2 r . Comparing Eq. (296) with 

9iiv = Vnv + h^u, (297) 

one can notice that the solution includes a static contribution. Furthermore, the 
freedom of the arbitrary coordinate transformation could be remained. We should omit 
these components in order to take the components of gravitational wave, which we are 
interested in. In other words, it is necessary to take TT part from the arbitrary h^. As 
a sideremark, let us summarize the procedure to do so in the following [247]. 

First of all, we perform the infinitesimal coordinate transformation (Eq. (267)) in 
order to satisfy the condition : h'°^ = 0. The condition is equivalently written, 

de 

= 4, = /loo -2-^, (298) 

and 

= K = hm - lj L- 1 JL. (299) 

From Eq (298), e can be determined and using it, ^ can be determined from Eq (299). 
As a result, one can confirm that ti = can be really satisfied. Then the problem is 
how one can find the TT components given the arbitrary tensor h^. There is a general 
formula to do so, in which the arbitrary tensor can be decomposed as follows, 

hij = hi? + (Wij + Wj,i - ^SijW'j) + S -fh\, (300) 
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S lm ^ = 0, 5«h™ = 0, (301) 

with Wi is an arbitrary longitudinal vector. Define h^ = hij — 5ijh l i/3, then Eq. (300) 
becomes 

W id + Wj,i - \^W\, = hij - hf?. (302) 



AW i + -(W l j) i = hi J , (303) 



By taking divergence of Eq. (302), Eq. (302) becomes 

1 

3 1 

because one can drop the TT components. By taking divergence Eq. (303) again, Eq. 
(303) reads, 

A(W^) = \h%. (304) 

To summarize, given hij, one can find W\ by solving the Poisson equation of Eq. (304). 
Then introducing the solution of W\ into Eq. (303), Eq. (303) becomes, 

AWi = h{, - \{W\{) t . (305) 

Solving the Poisson equation of Eq. (305), one can obtain Wi. Then with Eq. (302), 
TT component of h^ can be determined. 

The above procedure for extracting TT components from the arbitrary tensor is a 
little bit complicated, but is far easier for the plane- wave solution in Eq. (275). This is 
because the spatial derivative can be replaced by ik. The projection tensor to change 
the arbitrary vector to transverse is given by, 

Pij = Sij - TiiUj, (306) 

with 

n, = % (307) 
k 

Thus the transverse component of hy can be given by 

h[j = P t l Pj m h lm . (308) 

Noting P/Pij = Pij and P\ = 2, the procedure to make the tensor traceless can be 
given by, 

2 

When the source of gravitational wave is far distant from us, the plane-wave 
approximation in Eq. (275) is well satisfied. Thus we have only to do the procedure in 
Eq. (309), in order to extract the TT components. 



h TT = plpm him _ -jl(plm hlm) _ (309) 
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6.1.5. quadrupole formula In order to extract the components of gravitational wave, 
we should take TT components of Eq. (273). Let us write Eq. (273) again with using 
the geometrical unit (G — c — 1) for simplicity, 



h^(t,x) = 4 / d 3 x 



T^(t - 




X — X 


,x) 




X 


f 

— X 





(310) 

Integral with respect to x is performed in the source volume. Even for a supernova in 
our galactic center (r = 10 kpc), the distance to us (r) is far larger than the scale of the 
source, 

r = \x\ > \x'\ ~ L ~ O(10 km). (311) 

Note that we set L to be typical size of the inner core, because the gravitational wave 
is emitted most strongly at the epoch of core bounce as will be mentioned later. 

Neglecting the terms higher than 0([L/r} 2 ), the denominator of Eq. (310) becomes, 

(312) 



with 



\x — x \ ~r — n ■ x , 



x 

n = — . 

r 



(313) 



Taking TT of Eq. (310), Eq. (310) up to the lowest order of L/r becomes, 

= d 3 x T^ T (t -r + n-x,x), (314) 

here we remained the spatial components of h^ v in order to extract the gravitational 
wave. Let us expand Eq. (314) as follows, 



d 



T S T (t -r + n.x',x')=J2 d^ T S T (t ~ r > x ') 



(n ■ x'y 



m=0 



ml 



(315) 



, d 



= T< j 1 (t-r,x) + n-x-T i l j 1 (t-r : x) + - (316) 

This expansion is allowed only when the motion of the source is much slower than the 
speed of light. Writing the second term in Eq. (316) as follows, 



n-x'^ T (t-r,x) 



< 



d_, 
dt' 



n'(l-r.x) ~-7;f(*-r,a0,(317) 



where v is the typical velocity of the source, one can understand the reason clearly. This 
approximation is often referred as the slow-motion approximation. We employ this in 
the following. 

With Eq. (315), Eq. (314) simply reads, 



h T J = ~ £ n kl n k2 ■ ■ ■ n km H]? ^^ {t - r), 



m=0 



here 



2 / (9 \ m r 
13 m\ \dt) J K 



x kl x h2 ■ ■ ■ x km d 3 x. 



(318) 



(319) 
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Now let us take the lowest order (m = 0) in Eq. (319), 

Hij = 2 J Tij d 3 x. (320) 

Then we introduce the following identity, 

(T^xW)^ = (T^x» + T^x"),? = 2T, U , (321) 
which can be readily proved using the energy-momentum conservation, 

T^\ v = 0. (322) 
Using the identity, the right hand side of Eq. (320) becomes, 

J T ij d 3 x=^J d 3 x(T 00 x^) fi0 + 2 J d 3 x(rtV), M + J d 3 x(T H x'^) iH ) 
= l -(J d 3 x(T 00 x i x j ) fi0 + 2 J dS k (T ko x i x j ), + J rfS fc (T u xV),,) 
= \ J r/ :1 .r(V"".r',-M, )ll 

= \%pj l >y -' J l >r ( 323 ) 

here we used the Gauss's theorem from the second to the third column and we assumed 
the Newtonian perfect fluid, T 00 ~ p. Note that in the final column, represents the 
mass quadrupole moment. Taking the TT part of iy using Eq. (309), one can find the 
so-called reduced mass quadrupole moment, 

JJ T = J d 3 x ,,(,■',-' -^r). (324) 

Finally, we can reach to the quadrupole formula for the gravitational waves, 

2G-- TT/ t\ /• d 
— I ' 

c 4 r ■ 13 K 

where we have recovered G and c. 



W(t.*) = £W(t-l). C s s)- (325) 



6.1.6. angular dependence of quadrupole formula In this subsection, we discuss the 
angular dependence of the gravitational wave seen from the distant observer. 

For the purpose, it is convenient to find the non-zero components of hj^ in the 
spherical coordinate, 

x = r sin 9 cos 0, (326) 

y = r sin 9 sin 0, (327) 

z = rcos9, (328) 

not in the Cartesian coordinate. We assume that the gravitational wave propagates 
along the r direction. Due to the TT nature of the gravitational wave, h rT) h r $, and 



154 

h rt p vanish. By performing a simple coordinate transformation, for example for the 66 
component, 

TT dx l dxi 

h ee = % m m , (329) 
one can find the following nonzero components, 



hut) = r 



cos 2 6 + 1 
T 

= -h ge sin 2 6, (331) 



COS t7 ~r 1 1 

^ sin 2 _ ^« sin cos ^ cos - /i^, sin 6 cos # sin (330) 



^ cos # sin + h® cos # cos 20 



where 



r 2 sin# 2 ^ 

+ fc~ z sin # sin0 — hy Z sin # cos0, (332) 

h% = , (333) 

represents the second time derivative of the reduced quadrupole moments in the 
Cartesian coordinate. 

Now let us define the two independent components as follows, 

, hee , he,/, ,~~ A \ 

h+ = —, h x = — •~n' 334 

Then we perform a simple estimation for the gravitational wave emitted from the 
rotating star in axisymmetry. 

Given the density distribution of the rotating star of p(R, Z) in the cylindrical 
coordinate, the quadrupole moments are 

I xx = J p(R, Z) R 2 cos 2 4> R dR dz d(j) = n J p(R, Z) R 3 dR, (335) 
Iyy = J p(R, Z) R 2 sin 2 R dR dz d(p = I xx , (336) 

I xy = J p(R, Z) i? 2 sin0 cos0 R dR dz d(p = 0, (337) 

furthermore assuming the equatorial symmetry of the rotating star, the remaining 
components are, 

I zz = J p(R, Z) z 2 R dR dz dtp = 2tt J p(R, Z) z 2 R dR, (338) 

hz = Iyz = 0. (339) 

As easily understood, if the star rotates stationally, no gravitational waves are emitted 
because = 0. The gravitational waves are emitted when the rotating star contracts 
or expands dynamically, because the time derivatives of the quadrupole moments have 
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non-zero values. It should be noted that the gravitational waves can be emitted from 
the "axisymmetrically" rotating stars, when the motion is dynamically changing. 
From Eq. (330) and (332) with Eq. (340), the gravitational waves are 



h+ — (I X x 
r 



I zz ) sin 2 0, 



0. 



(340) 
(341) 



The gravitational waves are most strongly emitted in the direction perpendicular to the 
rotational axis (Eq. (340)). Intuitively, it is natural because the dynamical behavior 
of the rotating star can be seen most drastically for the observer in the direction 
perpendicular to the pole. On the contrary, gravitational waves from the merging 
neutron stars are most strongly emitted in the direction of the rotational axis. 



6.1.7. quadrupole formula for supernovae To end this section, we introduce the 
quadrupole formula in a useful form, which is often used for the computation of 
gravitational wave from core-collapse supernovae. 
First of all, let us define the tensor ft™, 

( \ 



Jij ~ 



\ 



Wlm Xi m 

X, m - sin 2 9W, 



lm 



J 



with 



lm ~ dcf>\de 
d 2 



w, 



lm 



d 1 d 2 \ 



1 d 2 

06 2 COt d0 sin 2 9 0(f) 2 



(342) 

(343) 
(344) 



where Y[ m is the spherical harmonics. 



After tedious calculations, one can check that hf^ can be expressed by using ft™, 



(OA)) 



16n 1 
IT 24' 



{h% + h% - 2fc£)Re(/j?(M)) 



^\h%Mf%{6A)) - sJ~h%Mf%{0A)) 



More compactly, the above equation may be written, 

m=2 



(345) 



(346) 



Although we have so far considered the gravitational wave up to the quadrupole, the 
gravitational wave in all order is shown to be expressed by the multiple expansions in 
the following form, 

°° ' 1 til rll 

ftf = EE; - + y im (t - r)<%) , (347) 
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where 



rim 



16tt /(/ + !)(/ + 2) \ 1/2 



T 00 Y lm *r l d 3 x, 



(2Z + 1)!!V 2(1-1)1 
represents the mass quadrupole, 

s i m= 327T ^(/ + 2)(2/ + l)V/2 



(348) 



(2/ + l)!!V2(/-l)(/ + i; 



J ejnX^-TojYj-^r 1 - 1 d 3 x(349) 



represents the mass-current quadrupole with Yj~ 1,lm * being the pure orbital spherical 
harmonics, and f\™, d l ™ are the pure-spin harmonics (see [330] for a complete 
description). Note that n\\ — n ■ (n — 2) • • • 1. 

In case of axisymmetry (m = 0), the gravitational wave up to the quadrupole 
(/ = 2) becomes, 

' d2 ;I 2 °(t-r)f™. (350) 



h 



TT 



v rdt 2 ' 

Noting in Eq. (348) that Y 20 = y^(3cos 2 # - 1) and T 00 = p, Eq. (348) reads 
G 32vr 3 / 2 



j20 _ 



nl POO g 

Jo ^ Jo drp<K 2^ ~ 2' )r4 ' 



(351) 



where \i = cos 9. And becomes, 
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fij 8V7T 




/o 







V 





sin 2 # 
sin 2 9 

Finally one can obtain the nonvanishing component in the following, 

1 



/ 



(352) 



h. = h, 



TT 



-h 



TT 



... i l l^e^l. 

8V 7T r 



where 



d 2 

aE2 — u T 20 
A ™ ~ d? 1 ■ 



(353) 



(354) 



It is noted that the numerical treatment of the second time derivatives in Eq. (354) are 
formidable. Using the hydrodynamic equations of perfect fluid, the time derivatives can 
be eliminated. For example, we take the first time derivative of Af 2 , 
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Using the equation of mass conservation expressed in the spherical coordinates, 
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In this way taking one more time derivative using the Euler equations, the well 
known form of A^q in the literature can be obtained, 

c V15 Jo Jo 

- 6v r v e fiy/1 - ll 2 - rd r $(3n 2 - 1) + 3<9 e $ fiy/l -fx 2 }, (358) 

where d r = d/dr, d e = d/d6. With Eq. (353) and Eq. (358), one can extract the 
gravitational waves from the numerical simulations assuming axisymmetry. 

Now that we have mentioned the physical foundations, we move now on to 
the discussion of the gravitational waves in core-collapse supernovae from the next 
subsection. The readers, who are interested in the detection techniques of gravitational 
waves by interferometric detectors, please see, for example [145] for a review. 
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6.2. Gravitational waves at core bounce 

If the gravitational collapse of the supernova core proceeds spherically, no gravitational 
waves can be emitted. The gravitational core-collapse should proceed aspherically and 
dynamically for the emissions of the gravitational waves. 

As mentioned in subsection 5.4, stars are generally rotating. This stellar rotation 
has been long supposed to play an important role in the gravitational waves from core 
collapse supernovae. The large-scale asphericities at core bounce induced by rotation 
can convert the part of the gravitational energy into the form of the gravitational waves. 
In this section, we review the gravitational waves emitted at core bounce in rotating 
supernovae (see also [250] for a review). 



6.2.1. characteristic properties First of all, we make an order-of-magnitude estimate 
of the amplitude and frequency of the gravitational waves emitted at core bounce for 
later convenience. 

A characteristic amplitude of gravitational waves at core bounce can be 
approximately estimated with the help of the standard quadrupole formula (see Eq. 
(325), for example, [271]) as follows, 
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where D is the distance to the source, is the second time derivative of the quadrupole 
moment of J^, M, R, and T dyn represents the typical mass and radius of the inner core 
and the dynamical timescale at core bounce, respectively. We assume that the supernova 
occurs at our galactic center at the distance of 10 kpc. e is a parameter, representing 
the degree of the nonsphericity as well as the degree of compaction, which may be 
optimistically estimated to be the order of 10 % in rapidly rotating supernova cores. 

In addition, the typical frequency of the gravitational waves, which can be 
approximately estimated by the inverse of the dynamical timescale, is expected to be in 
the following range, 

z/ GW ~ — *— ~ 0(100 Hz) ~ 1 kHz. (360) 

J- dyn 

We will later see that these values have the right order of magnitude. In the following, 
we review the study of gravitational wave at core bounce in rotating core-collapse 
supernovae. 



6.2.2. waveforms in rotating core-collapse supernovae Realistic progenitor models 
[120, 122, 363] (in section 2.2), which current researchers in this area can obtain, had 
been hard to access for those in the early studies. Some were forced to assume the 
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collapse of oblate spheroids with pressureless dust [328], others tried to include the 
effect of the internal pressure additionally, by which closer situations at core bounce 
were examined [284, 285, 286]. Although the obtained waveforms and the amplitudes 
of the GWs are quantitatively different from the ones in the current numerical studies, 
these pioneering studies were valuable in the sense that they constructed the formalism 
of the standard quadrupole formula still used by current studies [335], and that they 
obtained the qualitative understanding of the effect of initial angular momentum of the 
stars and the stiffness of equations of state on the waveforms near core bounce [286]. 
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Figure 114. Typical waveforms of gravitational waves (bottom) with the time 
evolution of the central densities (top) panels taken from [79]. Bottom left, middle, 
right panels correspond to type I, II, and III waveforms, respectively. The vertical 
dotted lines represent the epoch of core bounce. The peak spike and the subsequent 
spikes represent the gravitational waves emitted at core bounce and at the oscillation 
of the inner core produced by its inertia after core bounce. 

In the early 1980's, E. Miiller in the Max-Planck Institute for Astrophysics (MPA) 
performed two-dimensional (2D) axisymmetric core collapse simulations with better iron 
core models and calculated the quadrupole GW emission [240] . As for the microphysics, 
he employed a finite-temperature equation of state, however excluded the treatment of 
neutrino energy loss for simplicity. Due to the poor computational intensity at that 
time, only a small set of models could be investigated. However, it was found that 
differential rotation enhanced the efficiency of the GW emission. Afterward, developing 
the 2D hydrodynamic code employed in [240], Monchmeyer et al accounted for electron 
capture and treated neutrino transfer by making use of a leakage scheme for simplicity 
[237] . By computing four models changing the initial angular momentum parametrically, 
they categorized the obtained shapes of the waveforms into two distinct classes. The 
waveform categorized as Type I is distinguished by a large amplitude peak at core bounce 
and subsequent damping ring-down oscillations. This waveform is obtained when the 
initial angular momentum is small, which leads to core bounce near at nuclear density. 
A typical waveform for a Type I is shown in the left panel of Figure 114. They also 
found the waveform identified as Type II, which shows a several distinct peaks caused 
by multiple bounce (see the middle panel of Figure 114 for an example of a Type II 
waveform). From their study, it was found that the gravitational signals at core bounce 
are largest with amplitudes less than ~ 10~ 20 for a source at the distance of 10 kpc 
in the frequency range of 5 x 10 2 — 10 3 Hz. One may find these values are roughly 
in agreement with the one obtained in the simple order-of-the-magnitude estimation in 
Eq. (359). Note here that the waveforms shown in Figure 114 are from the study of 
Dimmelmeier et al [79] discussed later. 

Since a clear criteria determining the types of the waveforms might not be obtained 
by the study of Monchmeyer et al. due to their limited models, Zwerger et al. simulated 
the collapse of a large number of the initial models (78 models) with varying amounts of 
the initial rotation rates, the degree of differential rotation, and the stiffness of equation 
of state [372]. In order to make this large survey possible, they employed a simplified 
equation of state and did not take into account electron capture and neutrino transport. 
Their initial models were constructed in a rotational equilibrium by the method of [88]. 
In contrast, all the preceding studies, with the exception of [44], constructed the initial 
models just by adding the angular momentum to the spherically symmetric progenitor 
models by hand. In the study by [372], the rotational equilibrium was produced by 
a polytropic equation of state with the initial adiabatic index of T r = 4/3, and then, 
the core-collapse was initiated by dropping the adiabatic index down to 4/3, varing 
the values from T r =1.28 to 1.325. In addition to the cold part, the equation of state 
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consists of the thermal part and the stiff part in order to take into account the shock 
heating and the repulsive action of nuclear forces, respectively. 

With these computations, they found that the type of the waveforms was mainly 
determined by the stiffness of the cold part of equation of state, T r . The type I and II 
was obtained for the models with relatively softer (T r <~ 1.31) and stiffer equations 
of state (1.32 <~ T r ), respectively. In addition, they found a smooth transition from 
type II to type I while fixing other parameters, such as the initial rotation rate and the 
degree of differential rotation. They explained the cause of the transition as follows. As 
the value of T r becomes smaller, the core-collapse is enhanced. This results in the core 
bounce at the higher density. Since the typical interval between the multiple bounces 
should be an order of the dynamical timescale tdyn ~ W^/°> the higher density at 
bounce results in a shorter interval between the subsequent bounces. This makes the 
transition to type II to type I. As for the degree of the differential rotation, they did 
not find a large effect on the transition of the waveforms. In addition to the above 
waveforms, it is noted that they observed an another class of the waveform, the so- 
called type III (the right panel of Figure 114) for their models with the extremely lower 
values of T r = 1.28. Between the initial models constructed in rotational equilibrium and 
those not in rotational equilibrium, they found no significant changes in the waveforms. 
Employing the extensive sets of the initial models, they pointed out that the maximum 
amplitudes of the GWs were in the range of 4 x 10~ 22 <~ h <~ 4 x 10~ 20 for a source 
at the distance of 10 kpc with the typical frequencies between 500 to 1000 Hz. 

More recently, Kotake et al (2003) calculated the waveforms by performing 2D 
rotational core-collapse simulations, in which they employed a realistic equation of state 
(EOS) and took into account electron captures and neutrino transport by the so-called 
leakage scheme [174, 179]. Employing the-state-of-the-art equation of state, it was found 
that the typical frequencies and the amplitudes of the gravitational waves are consistent 
with the previous studies qualitatively and quantitatively. Furthermore they pointed 
out the importance of detecting the second peaks of the gravitational waves, because 
they will give us the information as to the angular momentum distribution of evolved 
massive stars. This is explained below. 

The waveform for a Heger's 15M rotational progenitor model studied in [174] is 
given in the left panel of Figure 115, which can be categorized into the type I waveform. 
While in the right panel of Figure 115, the waveform (type II) is given for the model, 
which has a cylindrical rotation law with the strong differential rotation, in which the 
initial angular velocity is assumed to yield to a quadratic cutoff at 100 km radius. 
Comparing the panels, one can see that the sign of the second peak is negative in the 
left panel, while it is positive for the model in the right panel. Note that the second 
peak is defined to be the place where the absolute amplitude is the second largest. They 
found that this characteristics that the sign of the second peak is negative is common 
to the models with the strong differential rotation and the cylindrical rotation law. 
The absolute amplitudes of the peak and the second peak studied are shown in Figure 
116. In addition to the first peak, the second peaks are also shown to be within the 
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Figure 115. Time evolutions of the amplitude of gravitational wave for the 
representative models taken from [174]. Note that the distance of the source is assumed 
to be located at the distance of 10 kpc. 



detection limit of the first L1GO for a source within 10 kpc. It seems quite possible for 
the detectors of next generation such as the advanced LIGO and LCGT to detect the 
difference of the sign. Therefore, it may be possible to obtain in this way the otherwise 
inaccessible information about the angular momentum distribution of evolved massive 
stars. 

They also discussed the relation between the maximum amplitudes of gravitational 
wave and the initial which is the ratio of the rotational to the gravitational 

energy. From Figure 117, it was found that the largest amplitude is obtained for 
the moderate initial rotation rate (i.e., T/|W^| init = 0.5%) when one fixes the initial 
rotation law and the degree of differential rotation. This is understood as follows. 
The amplitude of gravitational wave is roughly proportional to the inverse square of 
the typical dynamical scale, tdyn- Since tdyn is proportional to the inverse square root 
of the central density p, the amplitude is proportional to the density. As a result, the 
amplitude becomes smaller as the initial rotation rates become larger because the density 
decreases then. On the other hand, the amplitude is proportional to the quadrupole 
moment, which becomes larger in turn as the total angular momentum increases. This is 
because stronger centrifugal forces not only make the mass of the inner core larger, but 
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Figure 116. Detection limits of TAMA [8], first LIGO [333], advanced LIGO [357], 
and LCGT [188] with the amplitudes from numerical simulations [174]. The open 
squares represent the maximum amplitudes for all the models, while the pluses and 
the closed squares stand for the amplitudes of the second peaks for the models with 
strong differential rotation and cylindrical rotation law (negative sign of the peak) and 
for the other models (positive sign of the peak), respectively. Note that the source is 
assumed to be located at the distance of 10 kpc. 



also deform it. The amplitude of gravitational wave is determined by the competition 
of these factors. As a result, the amplitudes is found to become maximal for moderate 
initial rotation rates. This is also noticed by the earlier work by Yamada and Sato (1995) 
[367]. Combining the result by Finn [91], who pointed out by a perturbation technique 
that the amplitudes of the GWs are proportional to the square of the angular momentum 
h ~ J 2 in slowly rotating cases, the relation between the initial angular momentum and 
the peak amplitudes up to the rapidly rotating cases could be understood. 



6.2.3. effects of magnetic fields In addition to rotation discussed so far, Kotake et al 
(2004) investigated the effect of magnetic fields on the gravitational signals [179]. 

They extended the quadrupole formulae (Eq. (325)) in a form including the 
contributions from the electromagnetic fields. To begin with, I 20 in Eq. (354) should 
be replaced by M^ 2 as follows, 

^o 2 = |;<o 2 , (361) 
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Figure 117. Relations between T/|VF|i n i t and the peak amplitude |/i TT | m ax for all 
the models [174]. In the figure, "CS, SS, SL, CL seq" represent the model sequences 
whose initial rotation profiles differ (see [174] for details). Note that the distance of 
the source is assumed to be located at the distance of 10 kpc. 



where the mass quadrupole formula is given as 



M, E2 = 



G 32tt 3 / 2 



J dfij drp^^-^r 4 , (362) 



20 ~ VI5 

where p* represents the total energy density including the contribution from the 
magnetic field, 

By a straightforward, however tedious, calculation to replace the time derivatives by 
the spatial derivatives applying the continuity equation, the equation of motion, and 
the induction equation, 
dB 



P* = P+7, — o- (363) 



, )f Vx(vxB), (364) 

noting the divergence-free constraint (V • B = 0), Af 2 can be transformed into the 
following form, 

^20 = ^20 , quad + ^20 , Mag' (3Q5) 

where A^q quad is the contribution from the matter: 

A 20 , quad = T^ " 7ff ( l ^ / ^ dr pW{^ ~ 1) + V(2 - 3p 2 ) - v£ - 6v r V 9 
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^20 Mag = ^i? 7xB + ^20 Pm ^ s the contribution from the magnetic field: 
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^20 jxB' ^20 p m represent the contribution from j x B part and from the time derivatives 
of the energy density of electro-magnetic fields, respectively. Only the first time 
derivative of the magnetic fields is remained, because this is the leading order and 
the numerical treatments of the second time derivatives are formidable. Then the total 
gravitational amplitude can be written as follows, 



h tt = . TT 
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TT 



(369) 



where the quantities of the right hand of the equation are defined by Eqs. (353), (366), 
(367), and (368). Note that q r and qe in Eq. (366) represents the gravitational waves 
contributed from the artificial viscosity (see e.g., [237]). When one uses an artificial 
viscosity of von Neumann and Richtmyer, which is a most popular one, the concrete 
form of qi is, 



q t = V, [I 2 p (V 
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where % — r,9 with I defining the dissipation length. Using above quadrupole formula 
including contributions from the electromagnetic fields [368, 179], they calculated the 
waveforms by performing the 2D magnetohydrodynamic core-collapse simulations [179]. 

With these computations, they found that the amplitude is affected in the strongly 
magnetized models whose initial Ein/I W| is greater than 0.1 %, where E^/l W| represents 
the magnetic to the gravitational energy. This is natural because the amplitude 
contributed from the electromagnetic fields should be an order of 
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p c c* vseveral x 10 17 G/ 

with B c , p c being the central magnetic field and the central density near core bounce, 
respectively. Thus, strongly magnetized models, whose central magnetic fields at core 
bounce become as high as ~ 10 17 G, can affect the amplitude. 

It was furthermore found that the contribution of the electromagnetic fields changes 
in the opposite phase to the matter contribution (see Figure 118). Together with a 
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Figure 118. Waveforms for a model with the strongest magnetic field obtained in 
the magnetohydrodynamic core-collapse simulation [179]. In the figure, "quad","j x 
£?" represent the contributions from the mass quadrupole moment and from the 
electromagnetic field, respectively (see Eq. (369)). The total amplitude is denoted 
as "tot" . Note that the source is assumed to be located at the distance of 10 kpc. 



slight offset of the electromagnetic part, the negative part of the amplitude becomes 
less negative, while the positive part becomes more positive (see Figure 118). As a 
result, the peak amplitude at core bounce is found to be lowered by ~ 10 % They 
confirmed that the amplitudes of second peaks and the difference of its sign, from which 
one may know the information of the angular momentum of the core as mentioned, are 
still within the detection limit of the first LIGO for the galactic supernova, although 
the characteristics of second peaks are reduced by the incursion of the strong magnetic 
fields. 

6.2.4- e ff ec t s of realistic equations of state We turn to the effect of the equation 
of state (EOS) on the gravitational signals. Needless to say, EOS is an important 
microphysical ingredient for determining the dynamics of core collapse and, eventually, 
the gravitational wave amplitude. As a realistic EOS, Lattimer-Swesty (LS) EOS [187] 
has been used in recent papers discussing gravitational radiations from the rotational 
core collapse [258, 242]. It has been difficult to investigate the effect of EOS's on 
the gravitational signals because available EOS's based on different nuclear models are 
limited. Recently, a new complete EOS for supernova simulations has become available 
[287, 303]. The EOS is based on the relativistic mean field (RMF) theory combined 
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Figure 119. Waveforms (left panel) for the models with the relativistic EOS (the 
solid line labeled as MSL4) and with LS EOS (the dashed line labeled as MSL4-LS) 
and the relation between the central density and the effective adiabatic index 7 near 
core bounce (right panel). These figures are taken from [179]. 



with the Thomas- Fermi approach. 

By implementing these two realistic EOS's, Kotake et al. (2004) looked into the 
difference of the gravitational wave signals [179]. The left panel of Figure 119 shows the 
waveforms for the models with the relativistic EOS (model MSL4) or the LS EOS (model 
MSL4-LS). The maximum amplitudes for the two models do not differ significantly (see 
the left panel of Figure 119). The important difference of the two EOS's is the stiffness. 
As seen from the right panel of Figure 119, LS EOS is softer than the relativistic EOS, 
which makes the central density larger at core bounce and thus results in the shorter 
time interval between the subsequent bounces. On the other hands, softer EOS results 
in the smaller lepton fraction in the inner core, which reduces the mass quadrupole 
moments at core bounce. By the competition of these factors (see Eq. (359)), the 
maximum amplitude remains almost the same between the two realistic EOS's, while 
the typical frequencies of the gravitational wave become slightly higher for the softer 
equation of state. Furthermore it was found that the aforementioned type III waveform 
observed in a very soft EOS polytropic equation of state [372] does not appear when 
the realistic equations of state are employed. 
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6.2.5. non-axisymmetric simulations In addition to the above 2D simulations, several 
3D simulations have been computed. The first 3D hydrodynamic core-collapse 
simulations well beyond the core bounce was performed by [263]. The initial condition 
for their study was based on the configuration at several milliseconds before core 
bounce in the rapidly rotating 2D models of Zwerger et al [372]. In addition to the 
configuration, they imposed low mode (m = 3) density perturbation and followed the 
growth of the instability, where m denotes the azimuthal quantum number. They 
observed the three clumps merged into a bar-like structure due to the growth of the 
non-axisymmetric instability (see the left panel of Figure 120). In fact, their models are 
rapid rotator whose value of T/|W| at core bounce exceeds the critical value, beyond 
which MacLaurin spheroids become dynamically unstable against tri-axial perturbations 
(T/|W| > T/|Vr|dyn — 27.0%). However, they found that the maximum amplitudes of 
the gravitational waves were only ~ 2% different from the 2D cases by Zwerger et al. [372] 
(see the right panel of Figure 120). 

3D rotational core-collapse simulations of Fryer and his collaborators [101, 102] seem 
in favor of the above result. They have investigated whether the core fragmentation, and 
thus, the significant deviations of the gravitational radiation from 2D studies happen 
or not. As for the numerical computations, they have performed 3D smoothed particle 
hydrodynamic (SPH) simulations, with a realistic equation of state and the flux-limited 
diffusion approximation method for neutrino transfer. As for the initial model, they 
employed an rapidly rotating model with the initial value of T 1 / 1 | initia i of ~ 3 %. 
As a result, they also found that no fragmentation, although near core bounce the 
value of approaches to a critical value of ~14 %, beyond which the secular 

instability due to the non-axisymmetric perturbations sets in. Since the maximum 
value of T/ 1 | i n i t i a i predicted by the recent evolution models are <~ 0.5% [120], they 
concluded that the fragmentation or dynamical bar instabilities are unlikely to occur 
with any of the currently-produced supernova progenitors. 

While the past studies terminated the 3D simulations at several tens milliseconds 
after core bounce, Ott et al (2005) investigated the growth of the non-axisymmetric 
structure until the rather later phases (> 100 msec) after core bounce [259]. They 
found that the growth of the m — 1 mode, the so-called one-armed instability 
[65, 288, 289, 276], precedes the growth of the bar-mode (m = 2) instability (see 
Figure 121), where m denotes the azimuthal quantum number. Since the criterion 
of the growth of the one-armed instability is lower than that of the bar-mode instability 
[65, 288, 289, 276], they pointed out that the initial rotation rate required for the 
sufficient gravitational radiation from a galactic supernova enough to be detected by 
the future detectors can be as small as T j\W\ init = 0.2%, which is much smaller than 
the one previously assumed for igniting the growth of the bar-mode instability. 

6.2.6. general relativistic studies All the above computations employ the Newton 
gravity. In the following, we give a brief description of general relativistic (GR) studies. 
The study by [79] may be one of the representative GR studies for computing the 
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Figure 120. Gravitational radiation in a 3D model taken from [263]. Left panels 
shows the contour of density in the equatorial plane after about 3 msec after core- 
bounce showing the growth of the three arms. From the left panel, it can be seen that 
the difference of gravitational amplitudes between the 3D models (dashed and dashed- 
dotted lines) and the axisymmetric models (solid line) at core-bounce (30 msec) are 
too small to see by eye. Cross mode of the gravitational waves h x , which is of genuine 
3D origin begins to grow after core bounce (t > 33 msec), however, does not grow later 
on. (see [263] for details). 

GWs in stellar-collapse in the sense that they compared the properties of the GWs 
obtained by Newtonian and GR simulations systematically. As in the work of [372], they 
characterized the model difference by the degree of differential rotation, initial rotation 
rates, and adiabatic indices of equation of state and computed 26 initial models. The 
conformally flat (CF) metric was used to approximate the space time geometry in their 
GR hydrodynamic simulations. As well known, the CF approximation gives the exact 
solution of Einstein's equation in spherical symmetry. The approximation may not be 
so bad unless the configurations are extremely deviated from the spherical symmetry. 
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Figure 121. Time evolution of various quantities in a 3D model with T/| W|i n it = 0.2% 
calculated by Ott et al. (2005) [259]. Time is measured from the epoch of core bounce 
tb. Top panel shows that the amplitude of the m = 1 mode precedes that of the m = 2. 
Middle panel shows the time evolution of T/| W| and the core's maximum density. It 
is shown from the panel that after the epoch of t — t& ~ 100 ms, when the m = 2 
mode begins to be amplified, the transfer of the angular momentum becomes active 
which results in the increase of the maximum density and the decrease of the T/| W|. 
The bottom panel shows the gravitational strain at the distance to the source r as 
viewed down the rotational axis (solid curve) and as viewed along the equatorial plane 
(dotted curve) . One can see that the waveform traces the time evolution of the m = 2 
mode. Note in the panel that rh = 100 cm corresponds to h ~ 3 x 10~ 21 for a galactic 
supernova. 
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Figure 122. Prospects of detection of the gravitational wave signal from 
axisymmetric rotational supernova core collapse in relativistic (black filled circles) and 
Newtonian (red unfilled circles) gravity studied in Dimmelmeier et al [79] . The figure 
gives the (dimensionless) gravitational wave amplitude /i TT and the frequency range 
for all 26 models. For a source at a distance of 10 kpc the signals of all models are 
above the burst sensitivity of the LIGO I detector (except for some low amplitude, 
high frequency models), and well above that of the LIGO II interferometer. It can 
be seen that the typical frequencies of the gravitational waves are blue-shifted when 
relativistic effects are taken into account. This figure is taken from Dimmelmeier et al 
[79]. 
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Figure 123. Comparison of the waveforms between the fully general relativistic (solid 
line) and conformally flat approximation (dotted line) calculations. This figure is taken 
from Shibata and Sekiguchi (2004) [290] . 



However, they were forced to employ the quadrupole formula for extracting the GWs 
because CF approximation eliminates the GW emission from the spacetime. With these 
computations, it was found that the qualitative features of the GWs obtained in the 
Newtonian studies of [372] are almost true for their studies. Quantitatively, however, it 
was pointed out that relativistic effects make the central density at core bounce much 
higher than that in the Newtonian gravity. This is simply due to the enhancement of the 
gravity due to the GR effect. As for the peak amplitudes at core bounce, no significant 
differences between the GR and Newtonian case were found, while the typical frequencies 
at the peak amplitudes are blue-shifted for the GR models. This may be because the 
higher central density makes the timescale at core bounce shorter, which leads to the 
higher frequency (see Figure 122). 

Fully general relativistic collapse simulations from core-collapse to the formation of 
a neutron star have been performed by the group of M. Shibata [290] (see the reference, 
therein). It is mentioned that not only in their studies but also in other GR studies, the 
polytropic equations of state are employed in order to reduce the computational costs 
required for including a realistic equation of state. With these computations, it was 
found that not only the evolution of central density during core-collapse, bounce, and the 
formation of PNS, but also the waveforms are qualitatively in good agreement with those 
in the study of [79], except for a factor of ~ 2 difference of the amplitudes in the ring- 
down phase (see Figure 123). As a result, they concluded that the approximated method 
by [79] is appropriate for following the axisymmetric stellar core-collapse associated 
with the formation of the neutron stars and for estimating the emitted gravitational 
waves. Needless to say, the full GR calculations are indispensable for estimating the 
gravitational waves from the core-collapse of very massive stars associated with the 
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Figure 124. Results of the first 3D full GR calculations taken from Shibata and 
Sckiguchi (2005) [291]. Left panel shows the gravitational waveform of the model with 
T/|W|i n it ~ 1-8% with the strong differential rotation (model M7c2 in [291]). R + , x is 
the amplitude computed by the gauge-invariant extraction method, and ^4+, x is the 
one by the quadrupole formula (see [291] for details). Right panel shows the snapshot 
of the contour of the density on the equatorial plane at t = 90.7 ms with the vector 
fields (arrows), showing the bar- mode instability does develop, leading to the significant 
change in the amplitudes near the corresponding time. The maximum amplitude in 
the right panel can be translated into a dimensionlcss strain of h ~ 10 -18 for a galactic 
supernova with frequencies about 1 kHz. 

formation of the black hole. 

The numerics of the general relativistic studies have seen major progress recently. 
Dimmelmeier et al. recently succeeded in extending their 2D GR code to 3D with the 
CF approximation [80]. Shibata and Sekiguchi (2005) have performed the fully general 
relativistic 3D simulations and pointed out that the amplitudes of the gravitational 
waves can be enhanced by a factor of 10 than the axisymmetric collapse due to the 
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growth of the dynamical bar-mode instabilities, when the core initially rotates very 
differentially with rapid rotation of 1% < T/\W\ init < 2% [291] (see Figure 124). They 
discussed that the enhancement of the self-gravity due to the general relativistic effects 
results in a more efficient spin-up of the core, and thus, leading to the growth of the 
instability, which would be underestimated in the study of Rampp et al mentioned 
above. Very recently, fully general relativistic and magnetohydrodynamic simulations 
have also been reported to be practicable [81, 83]. The vary wide varieties of relativistic 
astrophysical events are expected to be investigated by the new-coming 3D GR studies. 
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6.3. Gravitational waves from convection and anisotropic neutrino radiation 

All the studies, which we reviewed so far, paid attention to the gravitational signals 
produced near core bounce due to the large-scale aspherical motions of matter induced 
by core's rotation without/with magnetic fields. In addition, two other sources of the 
GW emissions have been considered to be important in the later phases after core 
bounce, namely convective motions and anisotropic neutrino radiation, both of which 
can contribute to the non-spherical part of the energy momentum tensor of the Einstein 
equations. While gravitational waves from convective motions are originated from the 
aspherical motions of matter as well as the ones at core bounce, gravitational waves 
from neutrinos have some different features, which will be explained in the next section. 

6.3.1. foundation of gravitational waves from neutrinos As mentioned, the 
gravitational-wave signals at core bounce are emitted as bursts in which the wave 
amplitude rises from zero at core bounce, oscillates for several cycles and then approaches 
to zero (see for example Figure 114) due to the hydrodynamic motions of the central core. 
In addition, there is another class of gravitational wave, that is, bursts with memory, in 
which the wave amplitude rises from zero and then after the neutrino bursts settles into 
a non-zero final value [42] . The gravitational waves from anisotropic neutrino radiation 
from core-collapse supernovae is categorized to this class, which has been originately 
pointed out in late 1970's by [87, 335]. The detectability of such effect was discussed by 
[42] through the ground-based laser interferometers. According to [87], we summarize 
the formulation of gravitational waves from neutrinos in the following, which will be 
useful in the later discussions. 

At first, one should define a concrete form of an energy- momentum tensor of the 
neutrino radiation field to compute the gravitational waves. Then the following form is 
naturally assumed, 



radiation fields of neutrinos being released at the speed of light from the point x = to 
an observer at a distance of r. The functions L u (t) and f(Q,t) are the rate of energy 
loss and the angular distribution of neutrino radiation, respectively, at time t. Given the 
energy-momentum tensor, one can calculate the transverse-traceless (TT) gravitational 
field from the corresponding source as follows, 



in the observer's frame (x, y, z) (see Figure 125). For convenience, we write Eq. (372) 
as follows, 






(374) 
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Figure 125. Source coordinate system (x , y , z ) and observer coordinate system 
(x,y,z). The observer resides at the distant point on the z-axis. The viewing angle 
is denoted by £ which is the angle between z and z axis. The z axis coincides with 
the symmetry axis of the source, presumably the rotational axis. Central red region 
indicates the anisotropic neutrino radiation from a supernova. 



Introducing Eq. (372) to (373) and performing the integration with respect to r and 
we may obtain, 

In addition, we make use of the approximation that the gravitational wave signal 
measured by an observer at time t is caused by radiation emitted at time t = t — r. 
Hence we take t — t' = const = r, which means that only a neutrino pulse itself is 
assumed to cause a gravitational wave signal. This procedure is equivalent to eliminate 
the non-zero value of h v in case of isotropic neutrino radiation appeared in Eq. (375). 
We verify this in the later section. 

In the geometrical setup shown in Figure 125, we assume that the z-axis lies on the 
(x',z) plane for convenience. In this case, the two polarization states of gravitational 
waves satisfying the transverse-traceless conditions become 

<tt = K X = ~h™, (376) 

and 

K,tt = K y = K x , (377) 

in the observer coordinates. It is noted that the sum of the squared amplitudes 
\h + \ 2 + \h x \ 2 is invariant under the rotation about the z-axis. Using the following 
relations between the two coordinates, 

sin^ cos0 = sin # cos cos £ + cos 6 1 sin ^, (378) 
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sin 9' sin 0' = sin 9 sin 0, (379) 
cos O' = — sin 9 cos sin £ + cos cos £, (380) 
one can obtain the following, 



hl T + = - i " dt dQ'(l + cos0) cos20 L v (i)f(ti,i), (381) 

^ J — oo J An 

while the counter part of the amplitude, h x , is obtained just replacing cos 20 with sin 20, 
which immediately becomes zero by integrating over the angle due to the axisymmetric 
source we consider here. 

In the above equation, it should be noted that 9 and are required to be expressed 
in terms of the angles , 0' with respect to the source coordinate valuables, and the 
viewing angle of £. In the following, we consider two cases, in which the observer is 
situated parallel to the z axis (£ = 0) or perpendicular to the z axis (£ = tx/2). In the 
former case, one easily obtains, 

h T v l = - I ° dt dfi'(l + costf') cos2<// L v (t)f(Q',t), (382) 

r J — oo J An 

which becomes zero in case of the axisymmetric radiation source. Here the subscript p 
suggests that the observer is situated in the polar axis relative to the source coordinate 
frame. In the latter case, the observer is positioned perpendicular to the source's z axis 
(seen from the equator), and the field becomes, 

t— — 

= - I C dt j <m'*(&',<p') L u (t)f(Q',t), (383) 

r J — oo J An 

where, 

cos 2 $ — sin 2 ■&' sin 2 tp' 
cos 2 + sin 2 sin 2 ip' 

Since the amplitudes from neutrinos become largest seen from the equatorial plane of 
the source, we regard the Eq. (383) as the base formula to compute the amplitudes 
from neutrinos, we use the gravitational waves from neutrinos are anti-beaming 

Next, we make an order-of- magnitude estimate of the amplitude of the gravitational 
waves for neutrinos [87, 241]. From Eq. (383), one can estimate the amplitudes, h u , as 
follows, 

2,G rt~R/ c 

h u = — dt L v {t)-a{t'), (385) 



*(7?,^) = (l + sui7? cosy) — ^ ; - ; , 2 ^ f . (384) 



c 4 D _ 

where D is the distance to the source and a(t') is the time-dependent anisotropy 
parameter, 

a{t) = I dn'*(#',(p) f(ti,t), (386) 

J An 
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representing the degree of the deviation of the neutrino emission from spherical 
symmetry. Inserting typical values into the above equation, one can find the typical 
amplitudes, 

K - 1.6 x 10 2 cm I (A) r L " \ (¥) , (387) 
£> VO.l/ V10 52 erg s^ 1 / VI sJ v ; 

where we take an emission time of At = 1 s assuming constant radiation and the 

optimistic values of a ~ 0.1 suggested from the numerical results [241]. Thus it is 

expected that the gravitational- wave amplitude from neutrinos can be larger than the 

one emitted at core bounce in rotational core-collapse (see Eq. (359)). The typical 

frequency of the gravitational waves from neutrinos is expected to be lower by an order 

of magnitude than the one at the core bounce in rotational-core collapse, because the 

dynamical scale is not determined at the central core p ~ 10 14 g cm -3 but at the 

neutrinosphere p ~ 10 12 g cm -3 . 

It is noted that the gravitational memory stems from the change in the transverse- 
traceless part of the Coulomb-type (oc 1/r), and thus, can appear in other astrophysical 
events. Recently, the memory effect generated by a point particle whose velocity changes 
via gravitational interactions with other objects is studied [277]. Further, the memory in 
the context of jets in gamma-ray bursts is studied, which predicts that such gravitational 
waves are likely to be detected by the space-based laser interferometers such as LISA 
and DECIGO/BBO [275, 135]. 

Now we return to mention the gravitational waves from neutrinos and convections 
from core-collapse supernovae in the following sections. 



6.3.2. gravitational waves from convections and neutrinos in non-rotating stars As 
mentioned in section 5.3, the convections are likely to occur in the protoneutron stars 
and in the hot bubbles regions. In Figure 126, a typical GW waveform due to the 
convective motions and the associated anisotropic neutrino radiation inside the non- 
rotating protoneutron star (PNS) is presented. From the left panel, the time interval of 
the each GW signal from convections (thick line) is found to be very short with an order 
of milliseconds, while the waveform associated with the neutrinos (thin line) shows much 
less time structure. The short interval of the GWs from matter reflects the timescale 
of the convective motion inside the PNS, which may be roughly estimated as follows, 
tconv <~ -RpNs/^conv ~ O(ms) (i? PNS /20km)/(t; conv /l x 10 9 cm s _1 ) with R PNS and v com 
being the size of the PNS and the typical velocity of the convective motions. Since the 
amplitudes both from matter and neutrinos result from the small-scale motions induced 
mainly by the negative gradient of the lepton fraction, the amplitudes become much 
smaller than the ones at core bounce in rotational core-collapse (typically < 1/10). 
Due to the smaller amplitudes and the higher frequencies of the emitted gravitational 
waves, both unlike the ones obtained at core bounce in rotational core-collapse, they are 
marginally within the detection limits for the laser interferometer in the next generation 
(see the right panel of Figure 126). 
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Figure 126. Computations of the gravitational wave in the non-rotating protoneutron 
star (PNS) taken from [242]. The left panel shows the waveform due to the convective 
motion (thick line labeled as "flow" ) and the anisotropic neutrino radiation (thin line 
labeled as "z/s"). The insert shows an enlargement of the signal until ~ 100 ms from 
the start of the simulation. The right panel shows the total gravitational spectrum 
with the sensitivity curves for some detectors. The source is located at a distance of 
10 kpc. These figures are taken from Miiller et al (2003) [242]. 

Not only inside the PNS, but also outside the PNS, the convections are likely 
to be induced, as stated earlier (see section 5.3). Burrows & Hayes (1996) [51] 
performed 2D simulations, in which the density of the precollapse core was artificially 
reduced 15 % within 20 degree of the pole, and demonstrated how the initial density 
inhomogeneity affects the gravitational waves both from the convective motions and 
the anisotropic neutrino radiation. The large density inhomogeneity assumed in their 
work was predicted by the stellar evolution calculations, pointing out that they could be 
formed and amplified during silicon and oxygen burning stages [31, 112]. The obtained 
properties of the waveform is presented in Figure 127. They discussed that the total 
amplitude could be detected by the advanced LIGO, with a signal-to-noise ratio of 10, 
for a supernovae at a distance of 10 kpc. 

With almost the same motivation for the investigation, Fryer et al performed 
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Figure 127. The gravitational wave strain, K^J times the distance to the supernova, 
D, versus time. Core bounce is at 0,215 seconds. Each line shows the corresponding 
contributions to the gravitational waves. The contributions of matter motion and 
neutrinos to the GW amplitude can be seen of opposite sign at core bounce. Only for 
the first 20 ms after core bounce, the gravitational waves from neutrinos are shown to 
dominate over the mass motions. This figure is taken from Burrows & Hayes [51]. 



3D SPH simulations of the inhomogeneous core-collapse and discussed the waveforms 
[101, 102]. The computed signatures of the gravitational waves are consistent with the 
study of [51]. They discussed that such gravitational waves are within the detection 
limits for the advanced LIGO for the galactic supernova (see Figure 128). As shown 
by the simple order-of- magnitude estimates (Eq. (387)), it is seen from the right panel 
that the peak amplitude form neutrinos becomes as high as the one at core bounce 
[h ~ 10 -20 ) with the relatively lower peak frequency. 

Both in the above studies of [51] and [101, 102], a large density inhomogeneity 
(~ 0(10)% fluctuations in the density) prior to core-collapse is assumed in their initial 
conditions. On the other hand, it is noted that a recent study pointed out by the linear 
stability analysis in the cores of supernova progenitor stars that the timescale for the 
growth of the nuclear burning (the so-called e mechanism) is much longer than the time 
until the commencement of core-collapse, hence such a large inhomogeneity may not 
develop [233]. 

6.3.3. gravitational waves from anisotropic neutrino radiation in rotating stars 
Another possibility to induce the anisotropy of neutrino emissions is the stellar 
rotation. Recently, Miiller et al (2004) [242] performed the rotational core collapse 
simulations employing the elaborate neutrino transport with the detailed microphysics 
and calculated the gravitational waves from neutrinos. They found that the gravitational 
waves from the neutrinos grows due to convections and dominate over those of the matter 
at core bounce (see Figure 129). 

In their study, the initial models were limited to the rather slower rotating cases 
based on a recent stellar evolution models, in which the magnetic braking is taken 
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Figure 128. Gravitational from anisotropic neutrino radiation as a function of time 
and observer location (left) and its spectrum (right), obtained in the 3D simulations 
[103]. In the model, a 25 % core oscillation perturbation is assumed in the iron core 
of 15M Q progenitor star for producing the asphericity of the neutrino radiation field. 
From the right panel, it can be seen that the gravitational wave from the neutrinos 
peaks at lower frequencies and within the detection limits for the Advanced LIGO for 
the galactic supernova at the distance of 10 kpc. These figures are taken from Fryer 
et al [103]. 



into account [121, 122]. Recently, Kotake et al (2005) [176] calculated the waveforms 
from neutrinos employing a series of more rapidly rotating models with changing the 
rotational profiles and the degree of differential rotation parametrically in order to see 
the effects of anisotropy of the neutrino radiation induced dominantly by rapid rotation. 

In the left panel of Figure 130, the waveform for a typical model studied in [176] 
near core bounce is given. Note that the model is based on the recent stellar evolution 
calculation while excluding the effects of the magnetic fields [120]. About 2.1 ms after 
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Figure 129. Computations of the GW in rotational core collapse obtained by the 
state-of-the-art simulations (O = 0.5 rad s _1 is imposed on 15M Q progenitor model) 
[242]. In the left panel the waveforms contributed from neutrinos and the matter 
flows are shown. It is shown that the gravitational waves from neutrinos dominate 
over the ones from the flows almost always. The middle panel shows the spectrum 
of gravitational waves contributed from neutrinos. The right panel shows the total 
gravitational spectrum with the sensitivity curves for some detectors. Comparing the 
middle with the right panel, one can see that the gravitational waves in the lower 
frequency (< 100 Hz) are dominated by the neutrinos. These figures are taken from 
Miiller et al. (2003) [242]. 



core bounce at point A in the panel, the amplitude begins to rise more steeply than 
before (see point B in the figure). This epoch corresponds to the so-called neutronization, 
which occurs when the shock wave goes over the neutrino sphere. Here it should be noted 
that the neutronization occurs anisotropically in their rotating models. Due to the non- 
sphericity of the shapes of the shock wave and the neutrino sphere, the neutronization 
can occur more than once while it occurs only once for spherical models. In which 
direction the neutronization occurs is determined by the shapes of the deformed neutrino 
sphere and the anisotropically propagating shock wave. The first neutronization (at 
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Figure 130. Left panel is the waveform due to the anisotropic neutrino radiation 
for model S10. In the panel, points A, B, C, and D represents the epoch of core 
bounce (t = 242.7 ms), the onset of the first neutronization (t = 244.8 ms), the peak 
of the neutronization, which corresponds to the onset of the second neutronization 
(t = 245.3 ms), and the offset of the second neutronization (t = 246.1 ms), respectively. 
Right panel represents the angular dependence of <& in Eq. (388). Note that the angle 
is measured from the rotational axis. These figures are taken from [176]. 



point B in the panel) occurs near the pole at the radius of ~ 20 km along the rotational 
axis. The shape of the shock wave formed by core bounce is prolate initially because the 
bounce occurs near the rotational axis. This prolate shock wave crosses the neutrino 
sphere, whose shape is deformed to be oblate due to rotation. Since the area of surface, 
where the first neutronization occurs, is small, the neutrino luminosity is relatively low 
at this epoch (L v ~ 10 52 erg s -1 ). About 0.5 ms after the first neutronization at point 
B in the left panel of Figure 130, the amplitude shows sudden fall from t = 245.3 ms 
(point C) to t — 246.1 ms (point D). After the first neutronization occurs at the pole, 
the shock wave in the polar region is weakened by the neutrino energy loss and ram 
pressure of the infalling material. Then, the subsequent shock wave is formed in the 
vicinity of the equatorial plane and begins to move rather parallel to the plane. Since 
the first shock wave almost stalls along the rotational axis during a several millisecond 
after the first neutronization, the shape of the shock wave becomes rather oblate due to 
the propagation of the second shock wave rather parallel to the equatorial plane. As this 
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Figure 131. Waveforms due to the anisotropic neutrino radiations for some 
representative models with the shell-type (left panel) and the cylindrical rotation (right 
panel) profiles. Note that the source is assumed to be located at the distance of 10 
kpc. These figures are taken from [176]. 



oblate shock wave propagates along the equatorial plane, it crosses the oblate neutrino 
sphere. At this moment, the second neutronization occurs, where its luminosity becomes 
maximum, whose value reaches as high as < 10 54 erg s -1 . 

They explained the signatures of the gravitational waves at the two epochs of 
the neutronization as follows. At the first neutronization, the neutrino emissions are 
concentrated near the rotational axis. For convenience to understand the relation 
between the direction of neutrino emissions and the resultant properties of the 
waveforms, we analytically integrate over ip' in Eq. (383) and obtain the formula in 
the closed form as follows, 

where $(#) is the latitudinal angle dependent function, 

. /. / / f 2lT , / 1 + sin cos <£>' \ 

(see the right panel of Figure 130) and dL u (9,t)/dfl is the latitudinal angle dependent 
neutrino luminosity. From the formula, one can readily see that no gravitational waves 
can be emitted if the neutrino emissions are spherical. 
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From the feature of $ seen in the panel, the enhancement of the neutrino radiation 
near the rotational axis does not work sufficiently to change the amplitude although 
the amplitude shows a rise at the first neutronization due to the non-zero contributions 
of $ and the neutrino luminosity near the rotational axis (see from point B to C in 
Figure 130). On the other hand, at the second neutronization, the neutrino emissions 
are concentrated in the vicinity of the equatorial plane. The value of the function $ 
is negative (9 = n/2) and its absolute value is large compared to the value near the 
rotational axis (9 ~ 0) (see the right panel of Figure 130). In addition, the neutrino 
luminosity is much larger at the second neutronization. This is because the area of the 
surface at the second neutronization is quite larger than that at the first neutronization 
near the pole, since the oblate shock wave crosses the oblate neutrino sphere. As a 
result, the amplitude of the gravitational waves from the neutrinos shows a steep fall at 
the second neutronization. 

Moreover they demonstrated that the differential rotation mainly determines the 
waveforms. In Figure 131, the waveforms h u from near core bounce up to the final time 
of their simulation are presented for the representative models with the shell-type (left 
panel) and cylindrical rotation profiles (right panel), respectively. It can be seen from the 
figures that the amplitude due to the neutrinos in the later times becomes much larger 
for the models with the stronger differential rotation (compare S10 (weakest differential 
rotation with initial angular velocity cut of 1000 km) with SI (strongest differential 
rotation with initial angular velocity cut of 100 km), CS10 (weakest one) with CS1 
(strongest one), respectively). Simultaneously, it is found that this feature is regardless 
of the rotational profiles (compare the left (shell-type rotation) with the right panels 
(cylindrical rotation) in Figure 131). As the differential rotation becomes stronger, the 
shape of the shock wave becomes more prolate, which makes the neutrino emission more 
stronger in the vicinity of the rotational axis. This makes the gravitational amplitudes 
from the neutrinos more larger (see [176] for the more detailed explanations). 

In the left panel of Figures 132, the root mean square (rms) sensitivity curves with 
the rms gravitational wave spectra for the standard model is given. From the panel, 
it can be seen that the gravitational waves due to the anisotropic neutrino radiation 
dominate over those due to the mass motions at the frequency lower than several 10 
Hz. In the 21 models computed, it is found that the values of z/ eq ranges from 8 to 
58.1 Hz with the rms amplitude ranging from 7.28 x 10 -23 to 1.28 x 10~ 21 and that 
the values of /i^ eq generally become larger for the stronger differential rotation models. 
Here z/ eq represents a characteristic frequency, below which the gravitational waves from 
the neutrinos dominate over those from the mass motions, and the corresponding rms 
gravitational wave amplitude is named as h Ufiq . In the right panel of Figure 132, the 
values of h Ufiq with z/ eq for all the models are plotted with the sensitivity curves of laser 
interferometers. It can be seen that the gravitational waves from the neutrinos are the 
detection limits of the detectors in the next generation such as LCGT and advanced 
LIGO. Thus, it is suggested that the detection of the gravitational wave at the low 
frequency range (< 100 Hz) becomes more hopeful due to the contributions from the 
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Figure 132. Detection limits of TAMA, first LIGO, advanced LIGO, and Large- 
scale Cryogenic Gravitational wave Telescope (LCGT) with the expected gravitational 
wave spectrum obtained from the numerical simulations. The left panel shows the 
gravitational-wave spectrum contributed from neutrinos (solid) and from the matter 
(dashed) in a rotating model with fl = 4 rad s _1 imposed initially on a 15 Mq 
progenitor model. In the right panel, the open circles and the pluses represent the 
amplitudes of h Ujeq with the characteristic frequencies of v c(l for the models with the 
cylindrical and the shell-type rotation profiles, respectively. Under the frequency of 
f eq , the gravitational waves from the neutrinos dominate over those from the matter 
contributions. From the panel, it is seen that the gravitational waves from neutrinos 
dominate over the ones from the matter in a lower frequency (/ < 100 Hz). Note that 
the source is assumed to be located at the distance of 10 kpc. These figures are taken 
from [176]. 



anisotropic neutrino radiation. 
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7. Concluding Remarks 

The aim of writing this article was to provide an overview of what we currently know 
about the explosion mechanism, neutrinos, and gravitational waves in core-collapse 
supernovae. As we have discussed, much progress has been made in each topics. 

Recently, multidimensional studies and simulations of core-collapse supernovae have 
come into blossom again since 1990's when the direct observations of global asymmetry 
in SN 1987A were reported. The current trend might be ascribed to the fact that 
spherically symmetric supernova simulations have not yet produced explosions, albeit 
with the probably "ultimate" sophistication of the neutrino-transport method including 
the state-of-the-art microphysics. A step beyond the spherical models is in steady 
progress. It seems most natural and convincing to investigate the effect of the asphericity 
on the neutrino heating mechanism. Many ingredients to produce asymmetry in the 
supernova core have been considered, e.g., stellar rotation, magnetic fields, convection, 
and the standing accretion shock instability. In order to see the real outcome of them, 
we should be able to perform at least two-dimensional fully angle-dependent radiation- 
hydrodynamic calculations. In fact, several groups are really pursuing it with the use 
of advanced numerical techniques. 

Understanding the explosion mechanism of core-collapse supernovae is important 
not only for itself but also for the theoretical understanding for the other (astro)physical 
relevance, such as neutrino and gravitational-wave emissions. Conversely, we are 
now being able to understand them from the observations. In fact, neutrino and 
gravitational- wave astronomy are now becoming reality. 

Neutrino is a powerful tool to probe deep inside of the supernova while we can see 
just its surface by electromagnetic waves. In fact, supernova gave the first stage for the 
Neutrino Astronomy when we observed neutrinos from SN1987A. Neutrinos reflect the 
physical state of the core and even the density structure of the mantle if we consider 
neutrino oscillation. Actually neutrino oscillation is a necessary physics when we want to 
extract information on supernova from observed neutrinos because neutrino oscillation 
changes the neutrino spectra. However, this cannot be done so easily because there 
are some ambiguities in one of neutrino oscillation parameters, #i 3 , and mass hierarchy, 
although we have other parameters with high accuracies. Conversely, this means that 
supernova can be a unique laboratory which probes the unknown neutrino parameters 
which can not be studied by other experiments such as solar, atmospheric, accelerator 
and reactor experiments. 

Gravitational wave is also a powerful tool, which gives us the information deep inside 
the supernova core. If the supernova core rotates rapidly, the burst-like gravitational 
waves are emitted at core bounce, and their amplitudes are as strong as the currently 
running interferometric observatories could detect for a galactic supernova. While this 
is not the case, gravitational waves emitted in the later phase due to the convective 
motions and the anisotropic neutrino emissions will be most promising for the detections. 
Especially, the gravitational signal from neutrinos have different features from the 
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other, in the sense that they have memory, the detection may need somewhat different 
technique. If we could detect the gravitational waves from anisotropic neutrino emissions 
simultaneously with neutrinos themselves, it will bring about the great progress in 
the understanding of the explosion mechanism, because the anisotropy should play 
a important role. Moreover, the mutual understanding of the explosion mechanism, 
the supernova neutrinos, and the gravitational waves, which we reviewed somewhat 
separately in this article, will be greatly progressed. 

Moreover, supernova study is indispensable for the understanding of some 
hot astrophysical issues, such as the central engine of gamma-ray bursts and the 
nucleosynthesis in the explosion of population III stars. Further study will have a 
great impact not only on disclosing the mechanism of such astrophysical phenomena, 
but also on unveiling the fundamental properties of particle physics. 
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